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I. BACKGROUND AND PURPOSE 



This report and associated computer software applications 
are submitted in fulfillment of the thesis requirements for 
the degree of Master of Science, Aeronautical Engineering from 
the Naval Postgraduate School in Monterey, CA. The thesis 
text explains the fundamental theory and equations used in 
writing a set of software applications which predict an 
antenna's radiation parameters. Each application predicts the 
radiation characteristics of a given type of antenna over a 
planar surface which serves as a model of either earth or 
seawater. The thesis requirement was generated by a request 
from the Naval Maritime Intelligence Center ( NAVMAR INTCEN ) . 
At the request of NAVMARINTCEN, Mathcad mathematical software, 
DOS version 3.1, was used to write the applications. The 
NAVMARINTCEN request was extensive enough to be tasked as two 
separate thesis topics. This thesis fulfills the first half 
of the NAVMARINTCEN requirements. Mathcad applications are 
complete for the following types of antennas: 

1. Vertical Monopole Antenna 

2 . Elevated Vertical Dipole Antenna 

3. Elevated Horizontal Dipole Antenna 

4 . Arbitrarily Oriented Dipole Antenna 

5. Inverted L Antenna 

6. Long-Wire Antenna 

7. Terminated Sloping V Antenna 

8. Side- Loaded Vertical Half Rhombic Antenna 

9. Terminated Sloping or Horizontal Rhombic Antenna 

10. Sloping Double Rhomboid Antenna 

11. Vertically Polarized Log-Periodic Dipole Array 

12. Horizontally Polarized Log- Periodic Dipole Array 

13. Horizontal Yagi-Uda Array 
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II. INTRODUCTION 



New American Heritage Dictionary defines antenna as 'an 
apparatus for sending and receiving electromagnetic waves'. 
An electromagnetic wave is a time-varying, self -propagating, 
interrelated combination of electric and magnetic fields. 
Antennas radiate electromagnetic fields as a function of their 
time-varying surface charge density and surface current 
density distributions. There are many different types of 
antennas, each with its own set of radiation characteristics 
or radiation parameters . These radiation characteristics are 
related to an antenna's radiated electromagnetic field 
distribution, and they determine the useful applications for 
a particular antenna design. 

Current personal computers with high-speed processors can 
quickly and accurately predict antenna radiation parameters 
from electromagnetic theory. The Mathcad routines described 
in this report predict an antenna's radiation parameters based 
solely on its physical dimensions and electromagnetic theory. 
Accurately predicted radiation characteristics provide 
intelligence analysts with a reliable estimate of an antenna's 
capabilities, limitations, and vulnerabilities . Chapters III 
and IV are a review of applicable electromagnetic theory, and 
the remaining chapters describe the calculations of the 
associated Mathcad computer code, and are intended to be a 
user's guide for the associated Mathcad applications. 
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III. FUNDAMENTAL ELECTROMAGNETIC THEORY 



A. MAXWELL'S EQUATIONS AND HELMHOLTZ'S EQUATIONS 

Maxwell's equations are the basis for electromagnetic 
field calculations. Their differential form is given by [Ref 
1: pp. 321-325] 

V X E = (3.1a) V • E = -£ (3.1c) 

dt e 

VXH = J + (3.1 b) V • H = 0 (3.id) 

o t 

Equation 3.1a is Faraday's Law, 3.1b is Ampere's Circuital 
Law, 3.1c is Gauss's Law, and 3 . Id postulates the inexistence 
of magnetic charge. 

Maxwell's equations dictate the relationship between the 
electric and magnetic field intensities (E and H) and the 
charge/current density source distributions (p and J) (i.e., 
the antenna). Maxwell's equations are often applied in a 
source- free, current -free region, and an e )<Jt time dependence 
is assumed. The time -harmonic, free space Maxwell's equations 
are [Ref 1: p. 340] 

V X E = -ju\i 0 H (3.2 a) V * 2? = 0 (3.2c) 

V X H = jue 0 B (3.2b) V • H = 0 (3.2 d) 

where 

e Q =free space permittivity = — — • 10' 9 / farads \ 

0 36n \ meter ) 

li 0 =free space permeability = • 10" 7 | ~ met ex ) 
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Rather than solve Maxwell's equations directly, the 
method of vector potentials is often used where the scalar 
potential (V) and magnetic vector potential (A) can be 



obtained from the nonhomogeneous Helmholtz equations 



V 2 ^ + $ 2 V = (3.3) V 2 * + p 2 A = -|i 0 «7 (3.4) 



0 



where H=(l//x 0 ) (V X A) [Ref 1: pp. 338-340] . 

The solutions to Helmholtz's equations are the retarded 
scalar potential and the retarded vector potential. The value 
of V and A at some distance (R) from the source depends on the 
source's charge/current density at an earlier time (t-R/c), 
where c=l / v'/Xq e 0 is equal to the speed of light in free space. 
The delay (R/c) is the time required for electromagnetic waves 
to propagate through a distance, R, in free space from source 
to observation point. The solutions to the nonhomogeneous 
Helmholtz equations are [Ref 1: pp. 338-340] 



observation vector, and R' is the source vector. The wave- 
number is most often expressed as co/c = 2n/\, where X is the 
wavelength. The Mathcad applications use the expression 2i r/X 
to calculate / 3 . 




(3.5) 




(3.6) 



where the wavenumber (ft) is ciVV 0 e 0 in free space, R is the 
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The Mathcad applications described herein assume that the 
antenna under analysis is radiating in free space above a 
defined ground plane. The free space Maxwell's equations and 
nonhomogeneous Helmholtz equations (using the free space 
wavenumber) provide the theoretical basis for deriving the 
equations for an antenna's radiated electromagnetic fields. 
The radiated E and H fields can be calculated directly from 
their source's charge and current densities using Maxwell's 
equations. However, as previously mentioned, it is often 
simpler to find the magnetic vector potential using 
Helmholtz's equations. Once the vector potential is 
determined, the electric and magnetic fields are calculated 
from H=(l//x 0 )(V X A) and E-(l/jwe 0 )(V X H) [Ref 1: pp. 338- 
341] . Existing theoretical equations for the radiated 
electric fields of many antenna types have been derived in 
this manner. These equations are used in the Mathcad 
applications to predict an antenna's radiated electric field 
distribution. Other radiation parameters are then calculated 
on the basis of the predicted electric field distribution. 

B. ELECTROMAGNETIC FIELD REGIONS 

The space surrounding an antenna is divided into three 
regions: the reactive near- field, radiating near- field, and 
far-field [Ref 2: pp. 22-24]. The reactive near-field 
occupies the space immediately surrounding the antenna out to 
a radius of about 0.62- (DVX) 1 ' 4 , where X is the wavelength and 
D is the maximum dimension of the antenna. Reactive 
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electromagnetic fields dominate in this region. The radiating 
near- field occupies the space from the boundary of the 
reactive field out to a radius of about 2* (D 2 /X) . Radiation 
fields dominate in this region, and radiation patterns are 
often a function of both radial and angular coordinates. The 
far- field region occupies all space outside the radiating 
near- field which meets two additional far- field requirements: 
27TR/X >> 1 and R > 5D . Radiation fields also dominate in the 
far- field, but field components are primarily transverse, and 
radiation patterns are normally a function of directional 
variables only. The Mathcad antenna applications are valid 
only in the far- field region of a given antenna. 

C . RADIATION PATTERNS 

A radiation pattern is a three-dimensional representation 
of an antenna's radiated electromagnetic field distribution or 
power distribution at a fixed distance from the antenna. 
Because it is difficult to depict three-dimensional images, 
the patterns are most often plotted in a defined vertical or 
horizontal plane. The Mathcad antenna applications plot two- 
dimensional radiation patterns in polar coordinates depicting 
the far-field electric field distribution. Patterns are 
plotted in the 0=0 and 0 = tt/ 2 vertical planes, and also in a 
horizontal plane at an elevation angle set by an index from 
the applicable 'Elevation Angle Index Table' (Table 3.1 or 
3.2). Field magnitudes are normalized with respect to the 
maximum radiated electric field intensity. The magnitude to 
which each pattern is normalized is displayed below its plot. 
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TABLE 3.1 

ELEVATION ANGLE INDEX TABLE FOR NON-ARRAY ANTENNAS 



ELEVATION 


INDEX 




ELEVATION 


INDEX 


0.285° 


630 




46° 


470 


2° 


624 




48° 


463 


4° 


617 




50° 


456 


6° 


610 




52° 


449 


8° 


603 




54° 


442 


10° 


596 




56° 


435 


12° 


589 




58° 


428 


14° 


582 




60° 


421 


16° 


575 




62° 


414 


18° 


568 




64° 


407 


20° 


561 




66° 


400 


22° 


554 




68° 


393 


24° 


547 




70° 


386 


26° 


540 




72° 


379 


to 

00 

0 


533 




74° 


372 


30° 


526 




76° 


365 


32° 


519 




78° 


358 


34° 


512 




80° 


351 


36° 


505 




82° 


344 


O 

CO 

CO 


498 




84° 


337 


o 

o 


491 




86° 


330 


42° 


484 




88° 


323 


44° 


477 




88.857° 


316 
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TABLE 3.2 

ELEVATION ANGLE INDEX TABLE 
FOR LOG- PERIODIC AND YAGI-UDA ARRAYS 



ELEVATION 


INDEX 




ELEVATION 


INDEX 


0.57° 


314 




45.71° 


235 


2.28° 


311 




48 . 00° 


231 


4.57° 


307 




50.28° 


227 


6.86° 


303 




52.57° 


223 


9.14° 


299 




54.86° 


219 


11.43° 


295 




57.14° 


215 


13.71° 


291 




59.43° 


211 


16.00° 


287 




61.71° 


207 


18.28° 


283 




64.00° 


203 


20.57° 


279 




66.29° 


199 


22.86° 


275 




68.57° 


195 


25.14° 


271 




70.86° 


191 


27.43° 


267 




73 . 14° 


187 


29 . 71° 


263 




75.43° 


183 


32.00° 


259 




77.71° 


179 


34.29° 


255 




80.00° 


175 


36 . 57° 


251 




82.29° 


171 


38 . 86° 


247 




84 . 57° 


167 


41.14° 


243 




86.86° 


163 


43.43° 


239 




89.14° 


159 
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The Mathcad applications plot the far- field radiation 
patterns for an antenna's space wave, surface wave, and total 
radiated electric field. The space wave is comprised of the 
electromagnetic waves which propagate directly from antenna to 
the field point and electromagnetic waves reflected from the 
surface below the antenna. The surface wave is composed of 
electromagnetic waves ducted along the surface, much like 
waveguide propagation. The total field is composed of the 
space wave and surface wave combined. 

Space wave patterns typically vary with directional 
spherical coordinates, 6 and <p, but they are predominately 
invariant with distance from the antenna. The surface wave 
attenuates exponentially with distance, and far- field surface 
wave patterns are a function of distance from the antenna. 
Total field patterns are still predominately invariant with 
distance, because the maximum surface wave intensity is 
usually many orders of magnitude less than the maximum space 
wave intensity. Even at distances which just meet far-field 
requirements, the surface wave contribution to the total field 
is barely significant for parameter values of interest. 

D. POWER DENSITY, POYNTING VECTOR AND RADIATION RESISTANCE 

The instantaneous real power flux density of an electro- 
magnetic wave is called the Poynting vector ((P) , 

Qix.y.z) = E X H = 9te{E(x,y, z) e J&>c } X 9te{H{x,y, z) ©*»*} (3.7) 

= — 3le{rxir*} + — Sle{E X He 2jut ) 

2 2 
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The time invariant term on the right hand side of equation 3.7 
is the time average power flux density or average Poynting 
vector. Far- field E and H fields are related by the simple 
expression | E | = (l/r? 0 ) • | H | , where i? 0 =v Vo/ € o =1207r is the 
intrinsic impedance of free space. The average Poynting 
vector solely in terms of the radiated electric field 
intensity is given by [Ref 1: pp. 382-385], 

& uvv {x. yi z)=2-me{E XH*) = ±$le{E X iElUliJE-l! (3.8) 

The total power flux density is composed of both the 
Poynting vector and the imaginary part of (E X H*) which 
represents the reactive power of the radiated electromagnetic 
fields. Far- field power flux density is predominately real 
and is approximately equal to the average Poynting vector. An 
antenna's total average radiated power is calculated by 
integrating the poynting vector over a Gaussian surface which 
surrounds the antenna in the far-field [Ref 1: pp. 382-386]. 

P iad = P av9 = j s da = \jj*e{E XH'Y ds (3.9) 

The Mathcad applications use equations 3.8 and 3.9 for power 
calculations by integrating over the hemisphere which encloses 
the antenna with the ground plane as the lower boundary. 

Radiation resistance relates an antenna's total average 
radiated power to its peak input current. If a DC current 
equal to the rms input current were applied to a resistance 
equal to the radiation resistance, the dissipated power would 
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equal the antenna's total average radiated power [Ref 3: pp. 
47-48]. For the Mathcad antenna applications the radiation 
resistance is defined in terms of the total average radiated 
power and input current by 



■^rad 



2P 

K 



rad 



or 



P rad ~ " 



I J. |2 



rad 



(3.10) 



A sinusoidal input current (or voltage response across the 
input terminals for Yagi-Uda Arrays) with a maximum of unity 
is assumed in the Mathcad applications. 

E. DIRECTIONAL GAIN, DIRECTIVITY AND POWER GAIN 

An isotropic radiator is a hypothetical antenna which 
radiates electromagnetic energy equally in all directions. 
Its radiation pattern in any plane passing through the source 
is a circle. All real antennas are directional radiators 
which radiate or receive electromagnetic energy more 
effectively in some directions compared to others. Some 
antennas radiate an omnidirectional pattern which is 
directional in elevation, but non- directional in azimuth. 

Directive gain (D g ) is the ratio of an antenna's Poynting 
vector in a given direction, divided by the power flux density 
of an isotropic source with equal total average radiated 
power. Directivity (D 0 ) is the maximum value of the directive 
gain, and it is the primary measure of an antenna's 
directional properties. The power flux density of an 
isotropic radiator is equal to the total average radiated 
power divided by the surface area of a sphere at a given 
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distance from the antenna (47rR 2 ) . The expression for 
directive gain is [Ref 2: pp. 29-31] 

^ iso ^ iad 

The Mathcad applications calculate only directivity using 
9 (0 ,<b ) 4ni? 2 P (0 ,d) ) 

jj _ v v max • ^max ' _ * v max ' ^max 1 ^ ]_2) 

® iso P iad 

Power gain accounts for an antenna's radiation efficiency 
in addition to its directional properties. Power gain (G d ) is 
defined as 47rR 2 times the ratio of the Poynting vector 
magnitude in a given direction to the total input power to the 
antenna. Total input power includes losses due to ohmic power 
dissipation and impedance mismatches between transmission line 
and antenna. The ratio of total average radiated power to the 
total input power is the radiation efficiency , p eff . The 
relationship between power gain, directive gain, and radiation 
efficiency is given by [Ref 2: pp. 42-45] 

P eff=^ (3.13) 

^ in 

G d (Q , <]>) =p eff Dg(B , <J>) (3.14) 

= Peff z> ? _(e,4» = Peff D 0 (3.i5) 

Since radiation efficiency and power gain account for 
ohmic losses and impedance mismatches, these parameters are 
difficult to predict based on antenna dimensions alone. For 
these reasons, the Mathcad applications do not predict power 
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gain and radiation efficiency. Gain calculations are limited 
to directivity since it depends only on total average radiated 
power and the radiated electromagnetic field distribution. 
The maximum Effective Isotropic Radiated Power (EIRP) is also 
calculated from EIRP=P rad -D 0 [Ref 3: p. 62]. 

F. HALF -POWER BEAMWIDTH 

Radiation lobes are portions of the radiation pattern 
where a local maximum field intensity is bounded on either 
side by relative minimum intensities. The major lobe, or 
main beam, is the lobe in the direction of maximum radiated 
power flux density. Split beam antennas have more than one 
main beam. A minor lobe is any lobe except the main beam. 
Beamwidth refers to the half -power beamwidth of an antenna's 
main beam (or beams) , and it can be measured in the vertical 
or horizontal. For antennas with a defined beamwidth, the 
power intensity drops off on either side of the point of 
maximum radiated average power density until it is one half of 
the maximum value. The angle between the radials which 
intersect the half power points is defined as the half-power 
beamwidth. Main beams are not necessarily symmetrical since 
the radial which intersects the point of maximum power density 
does not necessarily bisect the half -power radials. 

The Mathcad applications do not calculate half -power 
beamwidth directly, but radiation patterns provide a good 
estimate. Since the radiation patterns depict the electric 
field distribution normalized to unity, a magnitude of 
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approximately 0.7 represents the half power-points of the 
pattern. If two radials are drawn from the origin through the 
half -power points, the angle subtended by the two radials is 
the antenna's half -power beamwidth. 

G. POLARIZATION AND BANDWIDTH 

Polarization describes the time varying properties of an 
antenna's radiated electric field at a set spatial coordinate 
as electromagnetic waves propagate outward from the antenna. 
Signals are usually transmitted with linear, elliptical, or 
circular polarization. Linear polarization is where the field 
vector at the spatial coordinate is directed along a fixed 
linear path. An elliptically polarized signal is one in which 
the tip of the electric field vector traces an ellipse around 
the fixed coordinate as the wave propagates outward. Circular 
polarization is a special case of elliptical polarization 
where the magnitude of the electric field remains constant as 
it rotates about the spatial coordinate. Clockwise rotation 
of the electric field vector is right-hand polarization, and 
counter-clockwise rotation is left-hand polarization. The 
Mathcad applications do not quantitatively predict the 
polarization at a selected spatial coordinate because of the 
computational intensity of such a calculation. Only the 
predominant type of polarization expected is addressed briefly 
in the introductory remarks of each application. 

The bandwidth of an antenna is the range of frequencies 
over which its radiation characteristics meet or exceed a 
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specified acceptable performance standard. Because radiation 
characteristics of dissimilar antennas vary differently with 
frequency, there are no set rules or guidelines which define 
bandwidth for all antennas in general. Bandwidth standards 
vary depending on antenna design and intended application. 
Therefore, the Mathcad applications do not predict operational 
bandwidth. Instead, multiple frequencies can be input for a 
given configuration to determine the bandwidth for whatever 
standards are used to define acceptable performance. 

H. EFFECTIVE LENGTH AND EFFECTIVE AREA 

The effective area, or effective aperture, gives the 
power delivered to the antenna load when multiplied by the 
incident average power flux density. It is defined by [Ref 2: 
pp 59-63] 



The maximum effective area is related to directivity for all 
antennas by [Ref 4: pp. 46-47] 



This is a theoretical maximum effective area since directivity 
does not account for radiation efficiency. Like directive 
gain and power gain, actual effective area and theoretical 
effective area are related by radiation efficiency 




(3.16) 



Anax 4^° 



(3.17) 




(3.18) 
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The maximum theoretical effective area is a good approximation 
of actual effective area for frequencies within the bandwidth 
of a well designed antenna with a high radiation efficiency. 

Effective length is defined as the ratio of the voltage 
induced across the antenna terminals by the incident electric 
field and is given by [Ref 4: pp. 40-42] 



1 ~ Vin 
±eff ~T* ~ 



Inc I 



(3.19) 



Effective length is related to effective area by 



■^max 



^racflmax 



(3.20) 



Substituting equation 3.17 into equation 3.20, we get the 
maximum theoretical effective length in terms of wavelength, 
directivity, radiation resistance and intrinsic impedance of 
free space as 




Finally, the actual 
theoretical effective 



■rad^ 2 A) 

4Jtr|o 

effective length 
length by 



(3.21) 

is related to the 



eff~ P ef f-^max (3.22) 

Effective area and effective length are actually a 
function of directional coordinates with respect to the 
source. The equations given here are for maximum values of 
effective area and effective length along the radial of 
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maximum directive gain (directivity) . While effective area is 
a parameter more generally applied to all antennas, the 
Mathcad applications use equations 3.17 and 3.21 to calculate 
both a theoretical maximum effective area and effective 
length. If an antenna's radiation efficiency for a given 
frequency is known, the actual effective area and length can 
be calculated by multiplying the theoretical values by the 
radiation efficiency. 

The next chapter on the dipole antenna explains in detail 
the derivation of the radiated electric field equations for a 
vertical dipole antenna. The elevated vertical dipole Mathcad 
application uses this equation to obtain the electric field 
distribution and total average radiated power (P rad ) based on 
an assumed sinusoidal current distribution on the antenna with 
a maximum of unity. Total average radiated power is then used 
to predict the radiation parameters discussed in the preceding 
paragraphs. In follow-on chapters, only the final electric 
field equations are presented because the derivation is 
similar for all types of antennas. They are primarily summary 
chapters which describe the calculations of each individual 
application, but any new or important concepts not presented 
in previous chapters will be discussed as required. 
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IV. THE VERTICAL DIPOLE ANTENNA 



A. THE ELEMENTAL DIPOLE 

The current element, or elemental dipole, is the 
building block for calculating antenna radiation parameters. 
The elemental dipole is a hypothetical differential length 
thin-wire antenna with a sinusoidally oscillating surface 
current. The antenna length (1) is many orders of magnitude 
smaller than the wavelength (X) , so the instantaneous current 
is considered uniform over the length of the dipole. Figure 
4 . 1 depicts the elemental dipole oriented along the z axis in 
a rectangular coordinate system with its center at the origin. 




FIGURE 4.1: Elemental dipole spatial orientation. 
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The surface current is given by 

i ( t) =1 cos (wt) =ffte{Ie j< * t } (4.1) 

and the magnetic vector potential, A, is found from equation 
3.6 to be 



A=— f 

4 nJc' R 4tz 




(4.2) 



Primed coordinates refer to source points, R is the distance 
from the origin to the field point, /3=o*/ /x 0 e 0 is the free space 
wavenumber, and the magnitude of the current density (J) is 
given by I [Ref 2: pp. 100-102]. The magnetic vector 
potential in spherical coordinates is 



A r =A z cos 



( 6 ) = 



Ho Idl 
471 



e -ipJ? 

R 



COS (0) 



(4.3a) 



A e =-A z sin (0) 



HoXdl 


— — sin (0) 


(4.3b) 


4k 


R ) 




o 

ii 


(4.3c) 



The magnetic field intensity, H, 



H =— V X A=a*— ^ 
Ho * Ho* 



3 (on) - 9Ar 

~dR (RAq) ~dd 



i^=-^P 2 sin (0) 
* 4n 



1 1 
+ - 



j$R (Jpi?) 



is found from 
(4.4) 

e-J'P* (4.5) 
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and the electric field intensity, E, is found from Ampere's 
Law by 



E= , 1 V X H =— 1 



jwe 0 

Idl 



j«e c 



^=-4^ T ioP 2cos (0) 



R 271 

Idl 



r i?sin(0) d0 
1 . 1 



9 [H^sin{Q)]-a e l-^(RH^) 



(jpi ?) 2 (j pi?) 



-jpR 



E e = ~— — T] 0 P 2 sin (0) 

471 



1 ■ +— L. +. 1 



. U&R) (Jpi ?) 2 (Jpi?) 



, -jp^ 



(4.6) 

(4.7a) 

(4.7b) 



H r -Hq-E^=0 (4.7 c) 

At far-field distances, the only significant components are 
the first terms inside the parenthesis on the right hand side 
of equations 4.5 and 4.7b. The far-field E and H fields are 
functions of a single angular coordinate (6) for a vertically 
orientated element [Ref 1: pp. 602-607] . 

Retaining only those terms that vary as l/R, we get the 
far-field electric and magnetic field intensities for the 
current element as 







e -^)psin(0) 



(4.8) 

(4.9) 



Note that the far-field relationship |h|=) 7 0 |e| holds true for 
the elemental dipole's radiated fields in the far-field 
region. 

The elemental dipole's radiated average Poynting vector 
(power density) is [Ref 1: 607-612] 
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«•)=*, ||B„| |/f ( ,| = i|^li no p2sin 2 (6) 



(4.10) 



The Poynting vector integrated over a Gaussian surface in the 
far- field yields the total average radiated power 

P rad =<f 9 ava ds= ( Idl \ 2 TloP 2 f 'T [sin 2 (6) ] sin (6) d6d4> (4.11) 

J s ° y 32n 2 J o Jo 

0 32n 2 \ 3 / 



P„ w =27t (Idi j 2 r|oP 2 



rad 32* 2 ' ,0 ' 



COS ' (6) -cos (6) 



.871 Udi) 2 B 2_ Udl) 2 fl2 
rad 3 327I 2 * oP 1271 110,5 



The directive gain and directivity are found from the 
power density and total radiated power as 



D, (6,4.) = <,8i % z 4nsin 2 (8) ,3 sln2 



rad 



8n 

3 



( 0 ) 



(4.12) 



D 0 = D g ^(Q,4>) = D J^j,4j=1.5 or 10 log 10 1.5=1.76 dB (4.13) 



Radiation resistance is found directly from P rad and equation 
2.10 as [Ref 1: pp. 607-612]. 



_ _ (Jdl ) 2 „ 2 

* zad 2 Kiad 1271 



(4.14) 
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l A. } J 
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(4.15) 
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B. THE ELECTRICALLY SHORT DIPOLE 



The derivation of the electric field for ah electrically 
short dipole is the next step in deriving the field equations 
for a practical finite length dipole. The electrically short 
dipole is a center fed thin-wire antenna oriented as in Figure 
4.1 with length greater than 0.02X, but less than 0.10X. The 
current distribution is estimated as triangular where 



I{x' , y' , z') 
where I m 





o <.z'<.— 
2 




(4.16 




<.z'<. 0 

2 



is the current maximum at the center feed. 



The dotted line in Figure 
4.2 is a more exact 
representation of the 
actual current distri- 
bution, but the triangular 
estimate is a reasonable 
approximation for short 
dipoles . 

The electrically 
short dipole has a magnetic vector potential given by 




A{x,y, z) =a. 



A=a. 
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(4.17) 
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Because dipole length is small in comparison to the far- field 
distance, R is considered constant throughout the integration 
over the length of the dipole. Note that the magnitude of the 
magnetic vector potential for the electrically short dipole is 
half that of the elemental dipole for the same current input. 

The equations for the radiated electromagnetic fields are 
derived in the same way as for the elemental dipole, and the 
final expressions for the electric and magnetic field 
distributions of the electrically short dipole are given by 



Directive gain and directivity for the electrically short 
dipole are the same as those of the elemental dipole. The 
magnitudes of the radiated power density, total radiated 
power, and radiation resistance for the short dipole are given 
by [Ref 2: pp. 109-112] 




(4.18) 



(4.19) 



E r =E^=H z =H e = 0 




(4.20) 




(4.21) 



R iad =20n 2 




(4.22) 
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C. THE FINITE LENGTH DIPOLE 



Finite length dipoles are center fed thin : wire antennas 
with negligible diameter (diameter << X and radius << length) . 
Dipoles with a length greater than 0.10X can no longer use the 
assumption of triangular current distribution with acceptable 
accuracy. The finite length dipole must account for the 
sinusoidal current variation along the length of the antenna. 
The sinusoidal surface current distribution for a finite 
length dipole oriented as depicted in Figure 4.1 is 



I{x' ,y', z') 



I(z') 



a t I 0 sin[p (h-z)], 

a x I 0 sin[P(A+zO]/ 



0 <.z'<,— 

2 

<.z'<. 0 

2 



(4.23) 



where I Q is the peak current input and h=— . 



Considering the finite length dipole as an aggregate of 
elemental dipoles, we get the expression for the electric far- 
field from equation 4.7b as [Ref 2: pp. 118-120] 




oP 

4ltR 




I(x' ,y ’ , z') e j ^ z dz 



/ 



(4.24) 



For purposes of the radiated electric field's phase 
angle, the radial field distance (R) cannot be considered 
constant while integrating over the length of the dipole. The 
variable r (which stands for |R-R'|) in the exponential 
(phase) term is approximated as [R- z ' cos (6 ) ] , where R is the 
distance from the origin to the field point, and z' is the 
source coordinate of the vertical dipole along the z axis 
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[Ref 2: pp. 112-114]. The 
approximation is depicted 
in Figure 4.3, and is 
accurate in the far- field 
for r > 21 2 /A. For 

purposes of calculating 
electric field magnitude, 
the variation of r over 
the length of the dipole 
is negligible, and the distance from the origin to the field 
point (spherical coordinate R) is used in the denominator of 
equation 4.24. Substituting equation 4.23 for I(x',y , ,z') 
into equation 4.24 and [R-z'cos(0)] for r in the exponential 
term of equation 4.24 and integrating, we get 




FIGURE 4.3: Phase approximation 

for distance variable, r. 



E e = J 






AtzR 



•sm 



(0)[|°X o sin[p (h+z , )]e + ^ z ' cos{6) dz / (4.25) 

f h J 0 sin[p {h-z 1 ) ] e *^ z>c ° 3 !0) dz' 

J 0 

sin (0) • 



=J 



Tl 0 P e ‘ jP * 



4 nR 



^J 0 sin[P [h+z 1 ) ][cos [pz'cos (0) ] +jsin [Pz'cos (0) ] ]dz ' 

+ /* h J 0 sin[P ( h-z ') ][cos [Pz^os (0) ] +j'sin [Pz^os (0) ] ]dz' 
J 0 



Since the current distribution is an even function, the 
odd term of j • sin (/? • z ' • cos (0) ) integrates to zero over the 
interval from -h to h. Since the cos (/? • z ' • cos ( 8 ) ) term is 
even, the integral can be written as 
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(0) f h sin[P ( h-z ') ]cos [Pz'cos (0) ] dz' (4.26) 



o 



After integrating, the final expressions for radiated electric 
and magnetic field distributions of a finite vertical dipole 
are given by [Ref 2: pp. 118-120] 



The coefficient outside the brackets is called the element 
factor, and the expression inside the brackets is the space 
factor. 

The electric field radiation pattern is a function of 
elevation ( 6 ) and wavelength (X) . For a fixed wavelength, the 
radiation pattern for any vertical plane through the antenna 
axis is obtained by plotting the magnitude of equation 4.27 in 
polar coordinates as 6 varies from zero to 2i\ . The pattern is 
normalized by multiplying the electric field magnitude by a 
factor such that the maximum result is one. Since electric 
field distribution is independent of <f>, the pattern for -7r<;0*;O 
is symmetrical with that for Os0s7T. Figure 4.4 depicts the 
normalized radiation patterns of four finite dipoles with 
successive lengths of X/4, X/2, 3X/4, and X. These lengths 
simplify calculation of the electric field patterns since /?h 
is 7r/4 , 7r/2, 37 t/ 4, and 7r for each successive dipole. If 
dipole length is greater than one wave-length, the patterns 




cosfPhcos (0) ] - cos(Ph) 
sin (0) 



(4.27) 




cosfphcos (0) 1-cosph 
sin (0) 



(4.28) 
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become more complex and the 
number of lobes increases. 




FIGURE 4.4: Vertical plane 
patterns for dipoles of length 
l-X/4 , l-X/2 , 



l=3X/4 and 1=X 



The radiated power 
density and total radiated 
power are found just as 
they were for the elemental 
dipole and the electrically 
short dipole. The final 
expressions for radiated 
power density and total 
radiated power are [Ref 2 : 
pp. 120-124] 



9 



T 2 
^0 



avg 






cos [phcos (9) ] -cos (Ph) i 2 
sin (0) 



p rad = {* n f* P avg R 2 s±n(Q)dQd4> 

^qIq 2 f* [cosfPhcos (6) 1-cos ( p h) l 2 
4 it Jo sin(0) 



(4.29) 



(4.30) 



The Mathcad applications evaluate equation 4.30 to obtain 
total radiated power for a sinusoidal input current with a 
maximum of unity. An expression for radiation resistance is 
also derived by Balanis [Ref 2: p. 124] . However, the Mathcad 
applications calculate radiation resistance from total average 
radiated power and equation 2.10. 

When the dipole is elevated above the origin along the z 
axis, an additional phase shift term of exp [- jj3H 0 cos (0) ] must 
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be included in equations 4.27 and 4.28, where H 0 is the 
elevation of the feed at the center of the dipole; hence, 




cos[phcos (8) 1 - cos(ph) ' 
sin (0) 



(4.31) 




cosfPhcos (6)1~cosPh 
sin (0) 



(4.32) 



D. GROUND PLANE EFFECTS ON THE RADIATED SPACE WAVE 

The electric field equations derived in Chapters 3A-3C 
are free space equations which do not consider the presence of 
a ground plane. Since antenna orientation is described with 
respect to the ground plane, orientation descriptions are not 
really valid until a ground plane is introduced. The free 
space electric field distribution for a dipole oriented along 
the x axis or y axis is identical to the distribution for the 
same dipole along the z axis. The radiation pattern is 
displaced to correspond to the new position of the dipole, but 
the radiated electric field distribution relative to the 
antenna is unaffected. 

The free space electric field equations no longer provide 
an accurate description of an antenna's radiated electric 
field distribution above a ground plane. The free space 
equations only provide the direct wave component of the space 
wave which propagates directly from the antenna to the field 
points above ground (-7r/2 < 6 < 7r/2) . They are not valid for 
|0| > 7 r /2 (points below the ground plane). 
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A complete description of the space wave above a ground 
plane must include the reflected wave component. The 
reflected component consists of electromagnetic waves radiated 
at elevation angles corresponding to |0| > 7r/2 in the free 
space equations. These waves strike the ground plane and are 
absorbed, partially reflected, or completely reflected before 
reaching a field point above the ground plane. The total 
space wave is obtained by adding the direct wave and reflected 
wave components for field points defined by -7r/2 < 6 < 7r/2. 
The distribution of the reflected waves vary relative to the 
antenna as its spatial orientation changes, so the electric 
field distribution (radiation pattern) is also affected when 
the antenna is repositioned with respect to the ground. 

The reflected wave component is determined using 
reflection coefficients and image theory. Reflection 
coefficients determine how much of an incident electric field 
is reflected by a given ground plane, and they are a function 
of 6 and the complex index of refraction. The complex index 
of refraction (n) in terms of frequency, conductivity (a) , and 
relative dielectric constant (e r ) of the ground plane is 









v/^o e o 



\ 



e,+ 



X J<*>6 0 \ 



e r~J 



. 18000a 



(4.33) 



-MHz 



The vertical reflection coefficient (r v ) 



n 2 c os (0 r ) -^n 2 -sin 2 (0 r ) 
n 2 c os (0 r ) +^n 2 - sin 2 (0 r ) 



(4.34) 
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is used to find the reflected wave contribution to the 6 



component of the space wave. The horizontal reflection 
coefficient (TJ 

r _ cos (6 r ) -^/n 2 -sin 2 (6 r ) ( 4 . 35 ) 

cos (0 r ) +^/n 2 - sin 2 (0 r ) 

is used to find the reflected wave contribution to the <£ 
component of the space wave [Ref 5: pp. 629-633]. 

Antenna image theory accounts for phase differences 
between the direct and reflected wave components. Phase 
differences are caused by the different distances travelled by 
the two components to a given field point. The image is a 
'mirror image 1 located below the ground plane at a distance 
equal to the antenna height above ground. The distance from 
the image to a given field point is the distance travelled by 
the reflected wave to that field point. This distance is 
approximated by R+H o cos(0), similar in derivation to the 
approximation illustrated in Figure 4.4 for the direct wave 
distance of R-H o cos(0). The distances travelled by the two 
components differ by 2H o cos(0), but the direct wave expression 
already includes the term exp [+j/3H 0 cos (6 ) ] , so the reflected 
wave expression includes the term exp [ - jj8H 0 cos (6) ] to account 
for any resulting phase differences. 

In equations 4.34 and 4.35, the variable 0 r is the angle 
between the z axis and the radial from antenna image to field 
point. If the height of the antenna is small in comparison to 
the far- field distance, 0 r is approximated as 6 (the angle 
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between the z axis and the radial from origin to field point) . 
The Mathcad applications use this approximation to find the 
reflection coefficients for 632 discrete values of 6 (316 for 
the array applications) in the interval -it/ 2 < 6 < 7 r /2 at a 
specified distance (R) from the antenna. 

The direct wave component of the space wave is calculated 
for discrete values of 6 with the free space electric field 
equation and the phase terms derived from the antenna height. 
The reflected wave components are then calculated for each 6 
by multiplying the direct wave components by their 
corresponding reflection coefficients and applicable phase 
terms. The reflected wave components are then added to the 
direct wave components to obtain the total radiated space wave 
at each 6 . The discrete results for the space wave are 
normalized with respect to the value with maximum magnitude, 
and the space wave radiation pattern is depicted by plotting 
the normalized field magnitudes for -tt/2 < 6 < tt/2 . 

The Mathcad applications plot radiation patterns assuming 
a planar surface below the antenna. The planar earth model 
provides accurate results for (1) antenna heights within 
several wavelengths of the surface and (2) wavelengths much 
smaller than the radius of the earth. The antennas addressed 
in this report are usually mounted well within several 
wavelengths of the ground, so the planar assumption in not a 
problem for the majority of calculations. If a configuration 
is encountered where the antenna is more than four or five 
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wavelengths above the ground, the Mathcad applications still 
provide a fair estimate of the radiation parameters. However, 
if highly accurate results are required, a computer program 
which accounts for the curvature of the earth should be used. 

The final expression for the radiated space wave of a 
vertical dipole antenna is [Ref 6: p. 166] 



Eq=J 



, 1)0*0® 
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cos[Phcos (8) 1 - cos(pfr) 
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j" e JPH 0 COS (6) + p e -JPH 0 COS (0) j 



=j’6 0 J 0 - 



,-jP* 



R 



cos [PAcos (6) ] -cos (PA) 
sin (0) 



^ jpH 0 cos (0) 



+r v e 



-jPH o cos(0)1 



(4.36) 



This is the equation used in the Elevated Vertical Dipole 
application to calculate the space wave radiation pattern and 
radiated electric field distribution. 



E. THE SURFACE WAVE 

The radiated electric field of an antenna above a ground 
plane includes a surface wave component in addition to the 
space wave. The surface wave radiation pattern is a function 
of the distance from the antenna and the relative dielectric 
constant and conductivity of the ground plane. Since the 
surface wave attenuates exponentially with increasing range, 
the surface wave radiation pattern varies with distance from 
the antenna. While the space wave can be calculated at 
different distances from the antenna, its radiation pattern is 
predominately invariant with distance from the antenna. 
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Equations for the radiated electric fields of a vertical 
antenna above planar earth were derived by Sommerfeld in 1909 
[Ref 7] , but the final form of these equations is complicated 
and difficult to evaluate. Norton simplified these equations 
to a form suitable for engineering applications [Ref 8] . 
Norton's expression for the surface wave of an elemental 
vertical dipole above a ground plane is [Ref 5: p. 644] 



j30$Idl(l-r v )F 



-jfp (J?+tf 0 cos (6) ] \ 



R 



(4.37) 



£(l-u 2 ) +£*sin (6) 1+ sin ^ I u v/l-u 2 sin 2 (0) 

k 2 j 



Its derivation is not explained here, but references 6 and 7 
give a detailed derivation. In equation 4.37, R is distance 
to the field point, r and k are unit vectors perpendicular and 
parallel to the antenna axis respectively, u is the inverse of 
the complex index of refraction (n) , and F is the surface wave 
attenuation factor. The function F introduces a surface wave 
attenuation dependent upon distance, frequency, and the 
properties of the earth over which the wave is travelling. 
The attenuation factor is given by [Ref 5: pp. 644-651] 



F e =l-jjTtP~ e e ~ P ‘ [erfc (j/P~ e ) ] 



p _ -j'P-Ru 2 [l-u 2 sin 2 (8) ] 
e 2 
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cos (0) 



u ^/l-u 2 sin 2 (0) 



or P e = 



-j'P [R+H 0 cos (0) ] 
2sin 2 (0 r ) 



v /n 2 -sin 2 (0 r ) 
cos (0.) + j — 

n 2 



(4.38) 

(4.39) 



(4.40) 
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where P c is the complex numerical distance for vertical 
polarization and the erfc term is the complementary error 
function defined by 

eifc(jJP e )=-^-{ e~ x2 dx (4.41) 

v 



Mathcad is incapable of evaluating the complementary error 
function for complex numbers, so an asymptotic expansion and 
infinite series approximation [Ref 9] was written into the 
Mathcad code to evaluate the erfc terms. 

When the magnitude of j •v / p 7 is less than 2.18, the 
Mathcad applications use 



erf(x+jy) = erf(x) +• 



2nx 



(l-cos(2 xy)) + j sin(2xy) 



(4.42) 
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n 2 + 4x 2 



f n (x,y) + j g n (x, y) 
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where 



f n (x,y) =2x-2x cosh (ny) cos (2 xy) +n sinh (ny) sin (2xy) 
g n {x,y) -2x cosh (ny) sin (2xy) +n sinh (ny) cos (2xy) 

|e {x, y) | = 10' 16 \erf {x+jy) \ 



an infinite series approximation which evaluates erf(jv / p7), 
where 'erf 1 is the standard error function. [Ref 9 p. 299]. 
The solution to erfc(jV^7) is then found by 

erfc{x+jy)=l-erf{x+jy) (4.43) 
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When the magnitude of j •V , p 7 is greater than 2.18, the 



Mathcad code uses 



fKze z2 er£c{z) *1+Y (- 1 )™ 1 ’ 3 ' 5 ' (2in 1} (4.44) 
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an asymptotic expansion which solves for j ♦V'ttP e • erf c ( jv^F7 ) 



[Ref 9 p. 298] . Either of these solutions can be substituted 
directly into equation 4.38 to find the surface wave 
attenuation factor. 

A magnitude of 2.18 is used as the transition point in the 
Mathcad applications because it provides the smoothest 
transition in both magnitude and phase between the series 
approximation and asymptotic expansion results. 

For antennas which radiate an electric field with a 4> 
component, the surface wave attenuation factor for horizontal 
polarization (F m ) must be found using 



where P m is the complex numerical distance for horizontal 
polarization. The variable 8 t in equations 4.40 and 4.46 for 
P e and P m , respectively is approximated as 6 (just as it was 
for the reflection coefficients in equations 4.34 and 4.35). 



F m =l-jf£P~ m e- p "erfc( jJFj 



(4 .45) 



and P = 

m 




(4.46) 
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The numerical distances (P c and P m ) and the surface wave 
attenuation factors (F e and F m ) are evaluated fbr O o s0s9O° to 
determine surface wave radiation patterns. Their value is of 
most interest, however, directly at the surface (0=90°) where 
the magnitude of F is called the ground wave attenuation 
factor. Numerical distance is proportional to distance from 
the antenna and the square of the frequency, while it is 
approximately inversely proportional to ground conductivity. 
When the magnitude of the numerical distance is less than 
about 4.5, empirical formulas show that the ground wave 
attenuation factor varies exponentially with numerical 
distance. For numerical distance magnitudes greater than 4.5, 
empirical formulas show that the ground wave attenuation 
factor is inversely proportional to the numerical distance. 
Numerical distance magnitudes are typically much greater than 
4.5 in the far- field. This implies that the ground wave 
(0=90°) field strength varies inversely with the square of the 
distance from the antenna since numerical distance is directly 
proportional to distance from the antenna and equation 4.37 
shows an explicit 1/R dependence in addition to the 
attenuation factor. 

The Mathcad applications determine the complex numerical 
distances for the same discrete values of 0 described in 
Chapter 4.D. After determining the complementary error 
function for each numerical distance result, the surface wave 
attenuation factor is calculated for each discrete 0. The 
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surface wave is evaluated for each 0 using the equations for 
the electric field distribution, and the discrete results of 
the surface wave are normalized with respect to the value with 
maximum magnitude. The surface wave pattern is depicted by 
plotting the normalized magnitudes for -7r/2 < 0 < ir/2 . The 
numerical distances (P e and P m ) at the surface (0=90°) are also 
listed with the predicted radiation parameters. 

With all of the terms defined, equation 4.37 is 
integrated over the length of a generic finite length vertical 
dipole to obtain the radiated surface wave distribution of an 
elevated vertical dipole [Ref 6 p.166] 



E e =j60I o- 



-j'P [R+H 0 cos (6) ] 



R 



d-rjF, 



' cos [ft.fr cos (6) ] -cos ( $h) 
sin 



( 0 ) 



(4.47) 



sin* ,0) -j^EH^Mcos (6) 



I2‘ 



The total radiated far- field electric field distribution 
of an elevated vertical dipole above a planar earth can now be 
expressed by combining equations 4.36 and 4.37 to write 



Eq- j C 0 1 o ■ 



e -jP[/;-.ff 0 cos(e)] |- CQs cos (Q) j _ cos (p^) i. 



R 



sin (0) 



(4.48) 



1 +f v e 



-i2pH 0 cos(6) 



+ ( 1 -r v ) F e e 



-j'2p// 0 cos (0) 



sin*(6)--£t^ sin2(9> cos(e) 



n- 



The first two terms inside the brackets represent the space 
wave and the third term represents the surface wave. 

Chapter IV has defined the necessary terms and provided 
a detailed derivation of the total electric field distribution 
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for an elevated vertical dipole antenna. The result is a 
radiated electric field distribution with a 6 component only. 
Most antenna types radiate an electric field with both 6 and 
4> components. While a detailed derivation of an electric 
field cf> component was not presented here, the necessary 
horizontal polarization terms associated with deriving $ 
components were discussed. When these horizontal terms are 
substituted for their vertical counterparts, the method used 
to derive the electric field equation for the <p component of 
a radiated space wave or surface wave is identical to that 
used in deriving the 6 component for the elevated vertical 
dipole . 

The remaining chapters of this thesis explain the antenna 
configuration and the calculations of each individual Mathcad 
application. Only the final expressions for the radiated far- 
field electric field distributions are provided since the 
derivation of each follows the same procedure as explained in 
this chapter. There is a lot of text repetition among the 
remaining chapters because each is written as a stand alone 
users guide for its associated Mathcad application. Any 
important concepts not presented in previous chapters are 
discussed as required. 
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V. THE ELEVATED VERTICAL DIPOLE 



The orientation of the elevated vertical dipole is 
depicted in Figure 5.1, where 1 is the length of the dipole, 
H 0 is the height of the feed above ground, R is the radial 
coordinate, 0 is the elevation coordinate, and <p is the 
azimuth coordinate. 




FIGURE 5.1: Spatial orientation of the elevated vertical 

dipole antenna for its corresponding Mathcad application. 



39 



The Mathcad application for the elevated vertical dipole 



requires the following inputs: 



1 length of the dipole (meters) 

f 5 operational frequency (Hertz) 

H 0 antenna feed height above ground plane (meters) 

R distance from antenna (meters) 

e r relative dielectric constant of ground plane 

a conductivity of ground plane 



The length of the dipole is used to determine the frequencies 
which correspond to a X/4, X/2, 3X/4, or X length dipole. The 
user inputs the frequency and the radiation parameters 
discussed in Chapter 3 are calculated. Feed height (H 0 ) must 
be at least equal to the half-length of the dipole (h=l/2), 
and the distance from the antenna (R) must meet the far- field 
requirements of Chapter 3.C. 

The equation for the radiated electric field distribution 
of the elevated vertical dipole was derived in Chapter 4 as 



Eq-J 6 0I Q ■ 



-j P [tf-tf 0 cos(9)] 
R 



/ cos [ph cos (8) ] -cos (ph) \ 
\ sin (0) / 



(5.1) 



1+ r v e- J ‘ PH ° cos(d) + (l-r v ) F e e - J ' 2PH » COS(0 jsin 2 (0) -i^Zz^iMcos (0) j 



The first two terms inside the brackets represent the space 
wave and the third term represents the surface wave. I 0 is 
one since a sinusoidal current input with a maximum of unity 
is assumed. 

The requested inputs are used to calculate variables at 
the top of page 41 for 632 discrete values of 6 which are 
equally incremented from -tt/ 2 to 7r/2. 
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h half-length of the dipole (meters) 

\ 5 wavelength of the operational frequency (meters) 

j8 free space wavenumber for operational frequency 

R d distance from antenna to field point (meters) 

R r distance from antenna image to field point (meters) 

n index of refraction 

P e complex numerical distance (vertical polarization) 

r v vertical reflection coefficient 

F e vertical surface wave attenuation factor 



The calculated variables are used to evaluate the 
radiated far- field space wave and surface wave for each 
discrete 6. The space wave and surface wave results are 
combined for corresponding values of 6 to obtain the total 
radiated electric field distribution. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each value of 6 to 
depict the radiation patterns. The patterns are valid for any 
vertical plane containing the antenna axis, because the 
radiated electric field does not vary with changes in <j>. 

Equations 3.8 and 3.9 are used to integrate equation 5.1 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna to find total average 
radiated power (P rad ) . With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation characteristics 
from the equations in Chapter 3: 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 6= 90°) 



Angle^ .elevation angle of maximum directive gain (degrees) 
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As an example, the Mathcad 

was executed with the following 

dipole length 
frequency 

height of antenna feed 
distance from the antenna 
relative dielectric constant 
ground conductivity 

Figures 5.2 and 5.3 depict the 

space wave and surface wave for 



vertical dipole application 
inputs : 

3 meters 
100 • 10 6 Hertz 
1 . 5 meters 
3000 meters 
15 

5 • 10‘ 3 

radiation patterns for the 
this example . 





The following radiation parameters were predicted by the 
Mathcad application for a sinusoidal current input of one Amp: 



Total power radiated 
Radiation resistance 
Directivity 

Effective isotropic radiated power 
Maximum theoretical effective area 
Maximum theoretical effective length 
Numerical distance (vertical) 
Elevation of directivity 



31.233 Watts 
62.467 Ohms 
6.034 

188.470 Watts 
4.322 sq meters 
1.692 meters 
195.178 

12.266 degrees 
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These results are consistent with expectations for this 
particular configuration. Appendix A contains computer 
hardcopies of additional example calculations for the vertical 
dipole and compares predicted radiation parameters to those 
expected based on previous calculations. 
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VI. THE VERTICAL MONOPOLE 



The orientation of the vertical monopole is depicted in 
Figure 6.1, where h is the length of the monopole, the antenna 
feed is at ground level, R is the radial coordinate, 6 is the 




FIGURE 6.1s Spatial orientation of the vertical monopole 
antenna for its corresponding Mathcad application. 
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The Mathcad application for the vertical monopole 



requires the following inputs: 



h length of the monopole (meters) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

€ r relative dielectric constant of ground plane 

a conductivity of ground plane 



The length of the monopole is used to determine frequencies 
which correspond to a X/8, X/4, 3X/8, or X/2 long monopole. 
The user inputs the frequency for which the radiation 
parameters discussed in Chapter 3 are calculated. Distance 
from the antenna (R) must meet the far- field requirements of 
Chapter 3.C. 

The monopole's electric field distribution is derived in 
the same way as the vertical dipole equation was derived in 
Chapter 4 and is given by [Ref 6: p.152] 



E Q =j3OI 0 - 



R 



A+jB -p A-jB 



sin (0) 



sin (0) 



( 6 . 1 ) 



. (l-r v ) (sin* (8) -J^! -sin z (6) cos (e) 

sin (0) ^ n 2 



-j2F e sin(Ph) 



sin 2 (6) - VglL ^(e) C0 3 ( e)|, v / "!- g j-i!L<§.) . 

) n 2 sin(0) 



n‘ 



A= cos [Ph cos (0) ] -cos (Ph) 
B=sin [ph cos (0) ] -cos (0) sin (Ph) 



The first two terms inside the brackets represent the space 
wave and the second two terms represent the surface wave . The 
vertical monopole's far- field electric field has a 6 component 
only , just as the vertical dipole. I 0 is one since a 
sinusoidal current input with a maximum of unity is assumed. 
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The requested inputs are used to calculate the following 
variables for 632 discrete values of 0 which are equally 
incremented from -7r/2 to 7r/2 : 



X 5 wavelength of the operational frequency (meters) 

/3 free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

r v vertical reflection coefficient 

F e vertical surface wave attenuation factor 



The calculated variables are used to evaluate the 
radiated far- field space wave and surface wave for each 
discrete 0. The space wave and surface wave results are 
combined for corresponding values of 0 to obtain the total 
radiated electric field distribution. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each value of 0 to 
depict the radiation patterns. The patterns are valid for any 
vertical plane containing the antenna axis, because the 
radiated electric field does not vary with changes in <£. 

Equations 3.8 and 3.9 are used to integrate equation 6.1 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna to find the total 
average radiated power (P rad ) . With the discrete values of the 
electric field and total average radiated power determined, 
the Mathcad application predicts the radiation characteristics 
listed at the top of the next page from the equations in 
Chapter 3 . 
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Radiation characteristics determined by the vertical 

monopole application are: 

Rn, d radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

maximum theoretical effective area (square meters) 

1,^ maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

Angles .elevation angle of maximum directive gain (degrees) 

As an example, the Mathcad monopole application was 
executed with the following inputs: 



monopole length 
frequency 

distance from the antenna 
relative dielectric constant 
ground conductivity 



0 . 5 meters 
150 -10 6 Hertz 
3000 meters 
15 

5 • 10' 3 



Figures 6.2 and 6.3 depict the radiation patterns for the 



space wave and surface wave for this example. 





FIGURE 6.3: Monopole surface 

wave radiation pattern. 



The following radiation parameters were predicted by the 
Mathcad application for a sinusoidal current input of one Amp: 
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Total power radiated 
Radiation resistance 
Directivity 

Effective isotropic radiated power 
Maximum theoretical effective area 
Maximum theoretical effective length 
Numerical distance (vertical) 
Elevation of directivity 



5.607 Watts 
11.215 Ohms 
31246 

18.203 Watts 
1.033 sq meters 
0.351 meters 



293.016 

26.815 degrees 



These results are consistent with expectations for this 



particular configuration. 



Appendix B contains computer 



hardcopies of additional example calculations for the vertical 



monopole and compares predicted radiation parameters to those 
expected based on previous calculations. 
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VII. THE HORIZONTAL DIPOLE 



The orientation of the horizontal dipole is depicted in 
Figure 7.1, where 1 is length of the dipole, H 0 is the height 
of the feed above ground, R is the radial coordinate, 8 is the 
elevation coordinate, and </> is the azimuth coordinate. 




x axi s 



Horizontal Dipole Orientation 



FIGURE 7.1s Spatial orientation of the horizontal dipole 
antenna for its corresponding Mathcad application. 
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The Mathcad application for the horizontal dipole requires 



the following inputs: 



1 length of the dipole (meters) 

H 0 antenna feed height above ground plane (meters) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

£ r relative dielectric constant of ground plane 

a conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The length of the dipole is used to determine frequencies 
which correspond to a X/4, X/2, 3X/4, or X length dipole. The 
user inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are computed. Feed height, H 0 , must be 
greater than or equal to zero, and distance from the antenna 
(R) must meet the far- field requirements of Chapter 3.C. The 
elevation angle index (d) sets the elevation coordinate, 6, 
for which a horizontal radiation pattern is determined. Table 
3.1 lists possible indices and their corresponding elevation 
angles from .285° to 88.857° in increments of about 2°. 
Indices between those listed can be used to interpolate a 
better approximation of a desired elevation. 

The horizontal dipole's radiated electric field is a 
combination of d and $ components. In addition, unlike the 
vertical dipole, its radiation pattern varies with changes in 
4>. The total radiated electric field is the vector sum of the 
6 and $ components. The electric field for the horizontal 
dipole is obtained in a manner analogous to that used for the 
vertical dipole discussed in Chapter 4 . The horizontal 
dipole's 6 and $ field components are given by [Ref 6: p. 144] 



50 



(7 . 1 ) 



E Q =-j60I 0 cos (4>) 

- j’P [R-H q cos ( 0 ) ] 



(8) cos (4>) ] -cos (PA) 
-sin 2 (0) cos 2 (<J>) 



( cos [pA sin i 
[ 1-si 

cos(8)(l-r v e'^ pH ° c0S(e) ) 

(1-r )Fe -Wcos(e) y/n 2 -sin 2 (8) f sin 2 (6) .jyfof-si n 2 (8) CQS (Q) \ 

n 2 l n 2 

g <h =j60J Q sin(4>)( cog tp^ sin (8) cos ($) -cos ( ph) N 
* \ 1-sin 2 (0) cos 2 (<{)) ; 



(7.2) 



,- j'P [ J ?- H 0 COS ( 0 ) ] 
R 



^ + P e -i-2p« 0 cos (6) + ^ ^ p e -J2PH 0 COS ( 0 ) j 



The first two terms inside the brackets of each equation are 
the space wave components and the third terms are the surface 
wave components. I 0 is taken to be one, since a sinusoidal 
current input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant </> of 4 >= 0 and $=7r/2 for 632 discrete 
values of 6 which are equally incremented from -ir/2 to tt/ 2 : 



h half-length of the dipole (meters) 

X 5 wavelength of the operational frequency (meters) 

13 free space wavenumber for operational frequency 

R d distance from antenna to field point (meters) 

R r distance from antenna image to field point (meters) 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
6 for <f)=0 and </>= 7 r/ 2 . The total space wave and surface wave 
distributions are determined from vector addition of 
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corresponding 6 and $ components . The space wave and surface 
wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the 0= 0 and 0= ir/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 6 to depict the radiation patterns. The Mathcad 
horizontal dipole application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and 0=7r/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 d ) as 0 varies from 0 to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 
components to the total average radiated power (P rad ) , equations 
3.8 and 3.9 are used to integrate equations 7.1 and 7.2 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. With the discrete values 
of the electric field and total average radiated power, the 
Mathcad application predicts the following radiation 
characteristics from the equations in Chapter 3 : 



52 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

A^ maximum theoretical effective area (square meters) 

l inax maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

P m numerical distance (horizontal polarization, 0=90°) 

Angle^ .elevation angle of maximum directive gain (degrees) 

The directivity (D 0 ) , effective isotropic radiated power 

(EIRP) , maximum theoretical effective area (A^) , maximum 

theoretical effective length (l max ) , and the elevation angle of 

maximum directive gain (Angle^J are all determined for both 

the 0= 0 and 0 = 7 r /2 vertical planes. 

As an example, the Mathcad horizontal dipole application 

was executed with the following inputs: 



dipole length 
frequency 

height of antenna feed 
distance from the antenna 
relative dielectric constant 
ground conductivity 
elevation angle index 



1 . 0 meter 
150 -10 6 Hertz 
1 . 0 meter 
3000 meters 
15 
10' 2 

535 (-27°) 



The following radiation parameters were predicted by the 



Mathcad application for a sinusoidal current input of one Amp: 



Total power radiated (Watts) 




26 . 116 


Radiation resistance (Ohms) 




52.233 


Numerical distance 


(vertical ) 




292.419 


Numerical distance 


(horizontal) 


66215.4 






0=0 


0 = 7T / 2 


Directivity 




1.310 


7.215 


EIRP (Watts) 




34.221 


188.435 


Max eff area (sq meters) 


0.417 


2.297 


Max eff length (meters) 


0.481 


1.128 


Angie.™ (degrees) 




45.071 


28.241 



Figures 7.2 through 7.7 are the space wave and surface 
wave radiation patterns in the 0=0 and 0 = 7 r /2 vertical planes 
and the designated horizontal plane for this example. 
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pattern for 0=0 vertical 
plane . 




pattern for 0=0 vertical 
plane . 




pattern for 0=7r/2 vertical 
plane . 




FIGURE 7.5: Horizontal 
dipole surface wave radiation 
pattern for 0=7r/2 vertical 
plane . 
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FIGURE 7.6: Horizontal 
dipole space wave radiation 
pattern for horizontal plane. 




pattern for horizontal plane. 



These results are consistent with expectations for this 
particular configuration. Appendix C contains computer 
hardcopies of additional example calculations for the 
horizontal dipole and compares predicted radiation parameters 
to those expected based on previous calculations. 
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VIII. THE ARBITRARILY ORIENTED DIPOLE 



The orientation of the arbitrarily oriented dipole is 
depicted in Figure 8.1, where 1 is length of the dipole, H 0 is 
the feed height above ground, a is the angle of the dipole 
axis with the horizontal, R is the radial coordinate, 6 is the 




FIGURE 8.1: Spatial orientation of the arbitrarily 

oriented dipole antenna for its Mathcad application. 
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The Mathcad application for the arbitrarily oriented 



dipole requires the following inputs : 



1 length of the dipole (meters) 

H 0 antenna feed height above ground plane (meters) 

a angle of antenna axis with the horizontal (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e r relative dielectric constant of ground plane 

o conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The length of the dipole is used to determine frequencies 
which correspond to a X/4, X/2, 3X/4, or X length dipole. The 
user inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are computed. Feed height, H 0 , must be 
greater than or equal to (1/2 ) * sin (a) , and the distance from 
the antenna (R) must meet the far- field requirements of 
Chapter 3.C. The elevation angle index (d) sets the elevation 
coordinate, 0, for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The radiated electric field for the arbitrarily oriented 
dipole is a combination of 0 and $ components, and its 
radiation pattern varies with changes in <p. The total 
electric field is the vector sum of the 0 and $ components. 
The arbitrarily oriented dipole's radiated electric field is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4, and is given by 
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_ -J'P* 

E 0 =-j6OI o ^-e 



( 8 . 1 ) 






1-cos 2 (v|f ) 
A 



(cos (a) sin (<{>) cos (0) -sin (a) sin (0) ) 



-r — (cos (a) sin (d>) cos (0) +sin(a)sin(0) ) e 

v 1-cos 2 (V) 



-j2$H 0 cos (0) 



+ (l-rj F e -J' 2 P H oCos(e) 



1-sin 2 (t|t ) 



[sin 2 (0) - yn a -sin 2 (8) CQS (9) 



n 4 



|cos(a) sin(4>) — — - -sin (a) sin(0) j 



E^j60I 0 ^-^- e J ‘ pHoCOS(0) cos (a) cos (<J>) 
R 



( 8 . 2 ) 






1-cos 2 (i|f) 1-cos 2 (iJrO 

where A^cos [$h cos (\|r) ] -cos ( ph) 
cos (i|/) =cos (0) sin (a) +sin (0) cos (a ) sin (4>) 
A 2 =cos [ph cos (ijf 7 ) ] -cos (ph) 



cos (t|r ) =sin (0) cos (a) sin (<|>) -cos (0) sin (a) 



Equation 8.1 is the d component (vertically polarized), and 
equation 8.2 is the (horizontally polarized) component [Ref 
6: p. 174-175]. The first two terms inside the brackets of 
each equation represent the space wave components and the 
third terms represent the surface wave components. I 0 is 
taken to be one, since a sinusoidal current input with a 
maximum of unity is assumed. 
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The requested inputs are used to calculate the following 

variables using a constant <f> of <f>=0 and 4>=tt/2 for 632 discrete 

values of 9 which are equally incremented from - 7 r /2 to 7 r/ 2 : 

h half-length of the dipole (meters) 

X 5 wavelength of the operational frequency (meters) 

j 8 free space wavenumber for operational frequency 

R d distance from antenna to field point (meters) 

R r distance from antenna image to field point (meters) 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 

The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
9 for <fi=0 and 0 = 7 r/ 2 . The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 9 and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 9 to obtain the total radiated electric field 
distribution for the <p=0 and 4>= 7 r /2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 9 to depict the radiation patterns. The arbitrarily 
oriented dipole Mathcad application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the </>=0 and 0 = 7 r /2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
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field components for the space wave and surface wave at the 
fixed elevation angle ( 0 d ) as <f> varies from O' to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those in the vertical planes. 

To find the contributions of the 9 and $ electric field 

components to the total average radiated power, P rad , equations 

3.8 and 3.9 are used to integrate equations 8.1 and 8.2 over 

the hemispherical Gaussian surface above the ground plane at 

a fixed radius (R) from the antenna. The sum of the two 

integrals is the antenna's P rad . With the discrete values of 

the electric field and total average radiated power 

determined, the Mathcad application predicts the following 

radiation characteristics from the equations in Chapter 3: 

R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

P m numerical distance (horizontal polarization, 0=90°) 

Angle max .elevation angle of maximum directive gain (degrees) 

The directivity (D 0 ) , effective isotropic radiated power 

(EIRP) , maximum theoretical effective area (Am**) , maximum 

theoretical effective length (l max ) , and the elevation angle of 

maximum directive gain (Angle max ) are all determined for both 

the 4 >= 0 and 0=ir/2 vertical planes. 
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IX. THE INVERTED L ANTENNA 



The orientation of the inverted L antenna is depicted in 
Figure 9.1, where H is length of the vertical antenna segment, 
1 is the length of the horizontal antenna segment, the feed is 
at ground level, R is the radial coordinate, 0 is the 
elevation coordinate, and 4 > is the azimuth coordinate. 




x axis 



Inverted L Antenna Orientation 



FIGURE 9.1: Spatial orientation of the inverted L antenna 

for its corresponding Mathcad application. 
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The Mathcad application for the inverted L antenna 



requires the following inputs: 



H length of vertical antenna segment (meters) 

1 length of horizontal antenna segment (meters) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e, relative dielectric constant of ground plane 

cr conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The combined length of the vertical and horizontal antenna 
segments (H+l) is used to determine frequencies which 
correspond to a A/8, A/4, 3A/8, or A/2 long inverted L 
antenna. The user inputs the frequency for which the 
radiation parameters discussed in Chapter 3 are computed. 
Distance from the antenna (R) must meet the far- field 
requirements of Chapter 3.C. The elevation angle index (d) 
sets the elevation coordinate, 9 , for which a horizontal 
radiation pattern is determined. Table 3.1 lists possible 
indices and their corresponding elevation angles from .285° to 
88.857° in increments of about 2°. Indices between those 
listed can be used to interpolate a better approximation of a 
desired elevation. 

The inverted L's radiated electric field is a varying 
combination of 9 and $> components, and its radiation pattern 
varies with changes in 0. The total radiated electric field 
is the vector sum of the 9 and $ components . The radiated 
electric field for the inverted L antenna is obtained in a 
manner analogous to that used for the vertical dipole 
discussed in Chapter 4. The radiated fields of each antenna 
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segment are considered separately, and superposition is used 
to determine the total radiated electric field.' The vertical 
segment radiates only a 0 field component, but the horizontal 
segment radiates both 0 and $ components. The electric field 
equations for the inverted L are [Ref 6: p. 169-170] 



■^ 0 , V~j 3 0I Q 



. -Jp* 



+ d-rjF, 



R 



K+jB^ A v -jB v 
sin(0) v sin(0) 



(9.1) 



fsitf <e> -JEEH*W cos (e , ) 

sm(0)( n 2 ) 

-[sin [p ( H+l ) ] -sin (pi) cos (pB cos (0) ] 
+jsin (P 1) sin [p/f cos (0) ] ] 

■]2F. [ sin 2 (6) - (9) cos (9) \ (?) \ 
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Equation 9.1 is the 6 component for the vertical segment, 
equation 9.2 is the 6 (vertically polarized) component for the 
horizontal segment, and equation 9.3 is the $ (horizontally 
polarized) component for the horizontal segment. The first 
two terms inside the brackets of each equation are the space 
wave components and any additional terms are the surface wave 
components. I 0 is taken to be one, since a sinusoidal current 
input with maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant 0 of 0= 0 and 0=7r/2 for 632 discrete 
values of 6 which are equally incremented from - tt /2 to 7r/2 : 



h half-length of the dipole (meters) 

X 5 wavelength of the operational frequency (meters) 

/? free space wavenumber for operational frequency 

R d distance from antenna to field point (meters) 

R r distance from antenna image to field point (meters) 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F c vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
6 for 0=0 and 0=7r/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 6 and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the 0=0 and 0=7r/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
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normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 0 to depict the radiation patterns. The inverted L 
Mathcad application computes the space wave, surface wave, and 
total electric field radiation patterns and radiation 
parameters in the 0=0 and 0=7r/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 d ) as 0 varies from 0 to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 0 and $ electric field 
components to the total average radiated power, P rad , equations 
3.8 and 3.9 are used to integrate equations 9.1, 9.2, and 9.3 
over the hemispherical Gaussian surface above the ground plane 
at a fixed radius (R) from the antenna. The sum of the 
integrals is the antenna's P rad . With the discrete values of 
the electric field and total average radiated power 
determined, the Mathcad application predicts the following 
radiation parameters using the equations from Chapter 3: 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0 = 90°) 

P m numerical distance (horizontal polarization, 0=90°) 



Angle max .elevation angle of maximum directive gain (degrees) 
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The directivity (D 0 ) , effective isotropic radiated power 



(EIRP) , maximum theoretical effective area (A,™,) , 
theoretical effective length (l max ) , and the elevation 
maximum directive gain (Angle^J are all determined 
the 0=0 and 0=7r/2 vertical planes. 



maximum 
angle of 
for both 
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X. THE SLOPING LONG -WIRE ANTENNA 



The orientation of the sloping long-wire antenna is 
depicted in Figure 10.1, where 1 is length of the antenna, the 
feed is at ground level, oi is the angle of the antenna with 
the horizontal, R is the radial coordinate, 6 is the elevation 
coordinate, and <f> is the azimuth coordinate. 



z axis 




FIGURE 10.1: Spatial orientation of the sloping long-wire 

antenna for its corresponding Mathcad application. 



The Mathcad application for the sloping long -wire antenna 



requires the following inputs: 



1 length of the antenna (meters) 

a angle of antenna with horizontal (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

6 r relative dielectric constant of ground plane 

a conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The user inputs the frequency (f 5 ) for which the radiation 
parameters discussed in Chapter 3 are computed. The distance 
from the antenna (R) must meet the far- field requirements of 
Chapter 3.C. The elevation angle index (d) sets the elevation 
coordinate, 6, for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The electric field equations for the sloping long-wire 
are given by [Ref 6 : p. 180-181] 



E Q =-j3 0I 0 - 



, -J'P* 



l-COS 2 (l]l) 



[cos (a) sin (<j>) cos ( 0 ) -sin (a) sin(<|))] ( 10 . 1 ) 
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1 -cos 2 (i|r ) (, n 2 
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n 2 
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p -J pj? 

E^=j30I 0 - — cos (a) cos (0) (10.2) 

■ ' A 1 + J B 1 +r A 2 + J^2 , (1 r )F A 2+j 5 2 

1-cos 2 (0) A 1-cos 2 (V) h m 1-cos 2 ( . 

where A 1 =cos[pj cos (0) ] -cos (pi) 

S 1 =sin[pi cos (0) ] -sin ( pi) cos (0 ) 
cos (i|r) =cos (0) sin (a) +sin(0) cos (a) sin (<|>) 

^ 2 =cos [pi cos (tyO ] -cos (Pi) 

I? 2 =sin [pi cos (ijr 7 ) ] -sin(Pi) cos (00 
cos (00 =sin(0) cos (a) sin(4>) -cos (0) sin(a) 

The radiated electric field for the sloping long-wire is 22 
obtained in a manner analogous to the vertical dipole equation 
discussed in Chapter 4. The long- wire's radiated electric 
field is a varying combination of 6 and $ components, and its 
radiation patterns vary with changes in <t>. The total radiated 
electric field is the vector sum of the 6 and 0 components. 
The first two terms inside the brackets of each equation are 
the space wave components and the additional terms are the 
surface wave components. I 0 is taken to be one since a 
sinusoidal current input with a maximum of unity is assumed. 

Requested inputs are used to compute the variables on the 
next page using a constant <p of 0=0 and 0=7r/2 for 632 discrete 
values of 6 which are equally incremented from - 7 r /2 to 7 r/ 2 . 



X 5 wavelength of the operational frequency (meters) 

(8 free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

T h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



69 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
8 for 0=0 and 0=?r/2 . The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 8 and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the 0=0 and 0=7r/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 6 to depict the radiation patterns. The sloping 
long- wire Mathcad application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and 0=tt/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle ( 6 d ) as 0 varies from 0 to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 8 and $ electric field 
components to the total average radiated power, P rad , equations 
3.8 and 3.9 are used to integrate equations 10.1 and 10.2 over 
the hemispherical Gaussian surface above the ground plane at 
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a fixed radius (R) from the antenna. The sum of the integrals 
is the antenna's P rad . With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation parameters using 
the equations from Chapter 3: 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

A max maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

P m numerical distance (horizontal polarization, 6 = 90 °) 



Angle max .elevation angle of maximum directive gain (degrees) 
The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A^J , maximum 
theoretical effective length (l max ) , and the elevation angle of 
maximum directive gain (Angle max ) are all determined for both 
the 0=0 and </>= 7 r /2 vertical planes. 
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XI. THE TERMINATED SLOPING V ANTENNA 



The orientation of the terminated sloping V antenna is 
depicted on the next page in Figure 11.1, where 1 is length of 
the antenna wires (numbered 1 and 2) , H 0 is the feed height 
above ground, H, is the height of the termination points above 
ground, y is the angle subtended by the wires, ij is the angle 
between the x axis and the projection of the antenna legs in 
the x-y plane, a is the angle between the plane which contains 
the antenna and the x-y plane, R is the radial coordinate, 6 
is the elevation coordinate, and 0 is the azimuth coordinate. 

The Mathcad application for the sloping V antenna requires 
the following inputs: 



1 length of the antenna legs (meters) 

H 0 feed height above ground (meters) 

H, height of termination points above ground (meters) 

y angle subtended by the antenna legs (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e r relative dielectric constant of ground plane 

a conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The user inputs the frequency (f 5 ) for which the radiation 
parameters discussed in Chapter 3 are computed. The feed 
height (H 0 ) and termination height (H,) can be any value 
greater than zero. The angle a is found from H 0 , H,, wire 

length (1) , and the half -angle between the wires (y/2) by 



a=sin‘ 1 



■ 

J ' cos (i). 



( 11 . 1 ) 
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application. 



The angle tj is then found from a and 7/2 by 



r| = tan ' 1 



tanl 



(i) 



( 11 . 2 ) 



cos (a) 

Distance from the antenna (R) must meet the far- field 
requirements of Chapter 3.C. The elevation angle index (d) 
sets the 0 coordinate for which a horizontal radiation pattern 
is determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
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increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the sloping V antenna is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The sloping V electric field 
equations are given by [Ref 6: p. 185-186] 



e -j p* 

£' 01 = 3OT o — — — \ cos (a) cos (0) cos (^-ri ) 



l-e- j * lu 1 
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1-e 
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e -j'2PH 0 cos(6) 



where U x = 1-cos (T x ) t/ 3 =l-cos ( Y 3 ) 

cos (T 3 ) =cos (0) sin (a) +sin (0) cos (a) cos ( 4>— ti ) 
cos (Y 3 ) =sin (0) cos (a) cos (<f>— T) ) -cos (0) sin (a) 
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where U 2 = 1-cos (T 2 ) C7 4 =l-cos (Y 4 ) 

cos (Y^ =cos (0) sin (a) +sin(0)cos (a) cos (4>-n) 
cos (T 4 ) =sin(0) cos (a) cos (<J>+r|) -cos (0) sin (a) 



The sloping V' s radiated electric field is a varying 
combination of 6 and $ components, and its radiation pattern 
varies with changes in <p. The total radiated electric field 
is the vector sum of the 6 and $ components. The fields of 
each antenna leg are considered separately, and superposition 
is used to determine the total electric field. Equations 11.3 
and 11.4 are the radiated 6 and $ components, respectively for 
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number one, and equations 11.5 and 11.6 are the 6 and $ 
components, respectively for wire number two. The first two 
terms inside each set of brackets represent the space wave 
components and any additional terms represent the surface wave 
components. I 0 is taken to be one, since a sinusoidal current 
input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant <p of 4 >= 0 and #=7r/2 for 632 discrete 
values of 6 which are equally incremented from -ir/2 to tt/2 : 



X 5 wavelength of the operational frequency (meters) 

/3 free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F c vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
6 for 4 >= 0 and 0=7r/2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding d and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the <p = 0 and 0=7r/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 6 to depict the radiation patterns. The sloping V 
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Mathcad application computes the space wave, surface wave, and 
total electric field radiation patterns and radiation 
parameters in the 0=0 and 0=tt/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 d ) as 0 varies from 0 to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 

components to the total average radiated power, P rad , equations 

3.8 and 3.9 are used to integrate equations 11.3, 11.4, 11.5, 

and 11.6 over the hemispherical Gaussian surface above the 

ground plane at a fixed radius (R) from the antenna. The sum 

of the integrals is the antenna's P rad . With the discrete 

values of the electric field and total average radiated power 

determined, the Mathcad application predicts the following 

radiation parameters using the equations from Chapter 3: 

R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

A^ maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 6= 90°) 

P m numerical distance (horizontal polarization, 0=90°) 

Angle max .elevation angle of maximum directive gain (degrees) 

The directivity (D 0 ) , effective isotropic radiated power 

(EIRP) , maximum theoretical effective area (A^J , maximum 
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theoretical effective length (l max ) , and the elevation angle of 
maximum directive gain (Angle max ) are all determined for both 
the 0=0 and 0=7r/2 vertical planes. 
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XII. THE SIDE-LOADED VERTICAL HALF RHOMBIC ANTENNA 



The orientation of the side -loaded vertical half rhombic 
antenna is depicted in Figure 12.1, where 1 is length of the 
antenna wires (numbered 1 and 2) , the feed is at ground level, 
oi is the angle between the antenna wires and the horizontal 
plane, R is the radial coordinate, 6 is the elevation 
coordinate, and <f> is the azimuth coordinate. 




Side-Loaded Vertical Half -Rhombic Orientation 



FIGURE 12.1: Spatial orientation of the side- loaded 

vertical half rhombic antenna for its Mathcad application. 



79 



The Mathcad application for the vertical half rhombic 



antenna requires the following inputs: 



1 length of the antenna legs (meters) 

a angle of antenna wires with horizontal (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e, relative dielectric constant of ground plane 

a conductivity of ground plane 

d elevation angle index (from Table 3.1) 

The user inputs the frequency (f 5 ) for which the radiation 



parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far- field requirements of 
Chapter 3.C. The elevation angle index (d) sets the 8 
coordinate for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The vertical half rhombic's radiated electric field is a 
combination of 8 and $ components, and its radiation pattern 
varies with changes in <p. The total radiated electric field 
is the vector sum of the 8 and $ components . The radiated 
electric field for the vertical half rhombic antenna is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The radiated fields of each 
antenna leg are considered separately, and superposition is 
used to determine the total electric field. The electric 
field equations for the half rhombic are [Ref 6: p. 198-200] 
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where U, and U 2 are given by 



C/ 1 =l-cos (T x ) ) t/ 2 =l-cos (T 2 ) 

cos (7 X ) =cos (0) sin (a) +sin (0) cos (a) cos (<t>) 

cos (T 2 ) =sin (0) cos (a) cos (<|>) -cos (0) sin (a) 

Equations 12.1 and 12.2 are the d (vertically polarized) and 
$ (horizontally polarized) components, respectively for leg 
number one, and equations 12.3 and 12.4 are the 8 and $ 
components, respectively for leg number two. The first two 
terms inside each set of brackets are the space wave 
components and any additional terms are the surface wave 
components. I 0 is taken to be one, since a sinusoidal current 
input with a maximum of unity is assumed. The requested 

inputs are used to calculate the following variables using a 
constant 0 of 0=0 and 0=7r/2 for 632 discrete values of 8 which 
are equally incremented from -tt/2 to 7 r/ 2 : 



X 5 wavelength of the operational frequency (meters) 

/J free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
8 for 0=0 and 0 =tt/ 2 . The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 8 and $ components . The space wave and surface 
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wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the 0=0 and 0=7r/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 6 to depict the radiation patterns. The vertical 
half rhombic Mathcad application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and 0= 7 r /2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 d ) as 0 varies from 0 to 27r in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 
components to the total average radiated power, P rad , equations 
3.8 and 3.9 are used to integrate equations 12.1, 12.2, 12.3, 
and 12.4 over the hemispherical Gaussian surface above the 
ground plane at a fixed radius (R) from the antenna. The sum 
of the integrals is the antenna 's P rad . With the discrete 
values of the electric field and total average radiated power 
determined, the Mathcad application predicts the following 
radiation parameters using the equations from Chapter 3: 
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R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0 = 90°) 

P m numerical distance (horizontal polarization, 0=90°) 



Angle max .elevation angle of maximum directive gain (degrees) 
The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A„ ax ) , maximum 
theoretical effective length (1^) , and the elevation angle of 
maximum directive gain (Angle max ) are all determined for both 
the </>=0 and <t>=ir/2 vertical planes. 
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XIII. THE TERMINATED SLOPING OR HORIZONTAL RHOMBIC ANTENNA 



The orientation of the terminated sloping or horizontal 
rhombic antenna is depicted on the next page in Figure 13.1, 
where 1 is the length of the antenna wires (numbered 1 to 4) , 
H 0 is the feed height above ground, H 2 is the height of the 
termination point above ground, y is the angle subtended by 
the feed wires (1 and 2) , rj is the angle between the x axis 
and the projection of the feed wires in the x-y plane, a is 
the angle between the plane which contains the antenna and the 
x-y plane, R is the radial coordinate, 6 is the elevation 
coordinate, and <p is the azimuth coordinate. 

The Mathcad application for the rhombic antenna requires 
the following inputs: 



1 length of the antenna legs (meters) 

H 0 feed height above ground (meters) 

H 2 height of termination point above ground (meters) 

7 angle subtended by the antenna legs (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e r relative dielectric constant of ground plane 

o conductivity of ground plane 

d elevation angle index (from Table 3.1) 



The feed height (H 0 ) and termination height (I^) can be any 
value greater than zero. The angle a is found from H 0 , H 2 , 
wire length (1) , and the half-angle between the wires (y/2) by 



a = sin 1 



^2—^o 

2l-cos(-^ 



(13.1) 
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or horizontal rhombic antenna for its Mathcad application. 



The angle rj is then found from a and y/2 by 



■H = tan 



-i 



tan (l 



(13.2) 



cos (a) 

The user inputs the frequency (f 5 ) for which the radiation 
parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far- field requirements of 
Chapter 3.C. The elevation angle index (d) sets the 0 
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coordinate for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the rhombic antenna is 
obtained in a manner analogous to that for the vertical dipole 
discussed in Chapter 4. The radiated field of each wire is 
considered separately, and superposition is used to determine 
the total electric field. The rhombic's electric field 
equations are [Ref 6: p. 183-191] 



E ei = 30I 0 \ cos (a) cos (0) cos ) 

R 



U 1 



(13.3) 



-iy i-*~ 3P2n * e -j**co»w> + ( 1 -r y) Fe lz±^l e - 

U, U-> 



j2$H 0 cos ( 6 ) 



|sin 2 (0) _fn 2 -sin 2 (8) cos (0) W n 2 -sin 2 (8) 

J n 2 cos(0) 



n‘ 



-sin (a) sin (0) 



l-e' inu ' r. x-e- mu > 



-j2pH 0 cos (8) 



a, 



U-, 



* d-rj f. 



l-e 






-j2$H 0 cos (6) 









sin 2 (0) 



\/n 2 - sin 2 (0) 



cos (0 






n‘ 



p-jpR 

= 3 0 J Q — cos (a) sin(<t>-Ti) 



l-e 



-jP-w, 



u, 



(13.4) 



r, l-e" jPi£/3 ^-J 2 PH 0 COS< 6 ) . r< \ c l-e- J * 2U2 ^-j 2 PH 0 cos(e) 
^ e U ^ C/ 3 
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e -jpi? 

Eq 2 =30 1 0 — - — \ -cos (a) cos (0) cos (<J)+t| ) 



R 



_~-iP-^ 2 



1-e 



l h 



(13.5) 



-r 



1-e c -j 2 pw 0 cos(e) + ,, p 1-e jf>Jg4 -j^ph 0 cos ( 0 ) 
U 4 v e U 4 

(sin 2 ( 6 ) -is i-3in 2 (8) cos (6) ) ^-sln 2 (9 ). 

I n 2 ) n 2 cos(6) 



+sin (a) sin (0) 



l-e~ jpjc/2 l-e~ jPJg4 c - j2ptf 0 cos (0) 



U 0 



U A 



d-r v )F e 



_~-iP 10 r t 



1-e 



U A 



e -j2pW o cos(0)| gin 2 (0) 



cos (8 

n 2 



~-j P* 

Eq i2 =-30 I Q — — — cos (a) sin(4>+ri) 



-ipju 4 



1-e 



-jp iu 2 



tt, 



+r h 1 ~ e „ — ! e -i2p« 0 cos(6) + 






£/, 



(13.6) 



-j2PH 0 cos<e) 






^e,3 = 30J 0 - 



>-iP R .* 



i? 



e’ jPiU2 \ -cos (a) cos (0) cos (4>-r|) 






1-e 



a. 



(13 



- -7p.Z£7, 1 -i&lU-, 

•p 1-e K ^-j^pHjcos (6) 4 _p N c, 1-e K ^-j^pHjcos ( 0 ) 

£/, U v> e £/, 



/sin 2 (6) -J^EUZH cos (6) )j^EH 
V n 2 ) n 2 c os 



-sin 2 (0) 
cos (0) 



+sin (a) sin (0) 



1 -e ^ pigl [ p 1-e Jjp ~ zg3 ^-j^p^cosfe) 






U, 



+ a-r v )F e - 



^-j3Pgicos(6)L in 2 (Q) . v/n 2 -sin 2 (0) 



£/, 



■cos (0 






>-j'P* 



g <| , ) 3 = -30 J 0 er ^' e' jpjg; cos (a) sin (4>-r|) 






1-e 



£/, 



(13.8) 



- -i&lU, „ -jpltf, 

+ P 1-e K 3 -j^p^costo) t /,_p s F 1-e _ -j 2 $H x cosm 

h n K h> m n 
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£e .4 = 30I 0 — 



rM* -i 



e -* p cos (a) cos ( 0 ) cos ( c|)+ti ) 



1-e J 



jpiy 2 



a, 



(13.9) 



_P 1-e jPit74 -^p^costei ^ M r \ ^ l-e" jPic/4 

-Iv F ® + (i-r v )F e — - — 



-j 2 $H 1 c os ( 0 ) 



2T 



fsin 2 (9)-i£l- si " 2(9) cos(e)j' / " 2 .~ sip2 . l9 I 



2 COS (0) 



-sin (a) sin (0) 



1-e J ' pJt7a .p, l-e~ jPic/4 ^ -j^p^cos (6) 
H, v U A 



+ d-rjF, 



1-e jPit74 - j^pHjcos (6) 



C4 



sin 2 (6) -^munilcos <e) 






F^ >4 = 3 0J 0 e e" jPil7l cos (a) sin(4>+ri) 



l-e’ jPiU2 



K, 



(13 . 10) 



^■p 1-e jPiU4 rt -j 2 Pflicos ( 0 ) / T _p N c, 1-e JpJC7<1 ^-J' 2 p« 1 cos( 0 ) 

+1 A FT e + FT e 









where U,, U 2 , U 3 , and U 4 are given by 



U ± = 1-cos <Y X ) t/ 2 = l-cos ( Y 2 ) 

£/ 3 =1-cos (Y 3 ) J7 4 = 1-cos (Y 4 ) 

where 

cos (Y^ =cos (0) sin (a ) +sin (0) cos (a) cos <<M) 
cos ( Y 2 ) =cos (0) sin (a ) +sin (0) cos (a) cos (<t>-Ti ) 
cos (Y 3 ) =sin(0) cos (a) cos (<j)-r|) -cos (0) sin (a) 
cos (Y 4 ) =sin(0) cos (a) cos ( 4>-*-ri ) -cos (0) sin (a) 

The rhombic's radiated electric field is a varying combination 
of $ and $ components, and its radiation pattern varies with 
changes in 0. The total radiated electric field is the vector 
sum of the d and $ components. Equations 13.3 and 13.4 are 
the d (vertically polarized) and $ (horizontally polarized) 
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components, respectively for wire one, and equations 13.5 and 
13.6 are the d and $ components, respectively for wire number 
two. Equations 13.7 and 13.8 are the analogous equations for 
wire three, and 13.9 and 13.10 are those for wire four. The 
first two terms inside each set of brackets represent are the 
space wave components, and any additional terms represent the 
surface wave components. I 0 is taken to be one, since a 
sinusoidal current input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant 0 of 0=0 and 0=7r/2 for 632 discrete 
values of 6 which are equally incremented from -7r/2 to tt/2 : 

X 5 wavelength of the operational frequency (meters) 

/3 free space wavenumber for operational frequency 

n index of refraction 

P c complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

T v vertical reflection coefficient 

r h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 

The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
6 for 0=0 and 0=7r/2 . The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 0 and $ components . The space wave and surface 
wave results are then combined for corresponding discrete 
values of 0 to obtain the total radiated electric field 
distribution for the 0=0 and 0=7r/2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
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normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 0 to depict the radiation patterns. The rhombic 
antenna Mathcad application computes the space wave, surface 
wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and <p=ir/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
selected elevation angle (0 d ) as 0 varies from 0 to 2ir in 632 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the 6 and $ electric field 

components to the total average radiated power, P rad , equations 

3.8 and 3.9 are used to integrate equations 13.3 to 13.10 over 

the hemispherical Gaussian surface above the ground plane at 

a fixed radius (R) from the antenna. The sum of the integrals 

is the antenna's P rad . With the discrete values of the electric 

field and total average radiated power determined, the Mathcad 

application predicts the following radiation parameters using 

the equations from Chapter 3: 

R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

A max maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 6= 90°) 

P m numerical distance (horizontal polarization, 0=90°) 

Angle max .elevation angle of maximum directive gain (degrees) 
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The directivity (D 0 ) , effective isotropic radiated power 



(EIRP) , maximum theoretical effective area (A^) , 
theoretical effective length (l max ) , and the elevation 
maximum directive gain (Angle,^) are all determined 
the <f>=0 and $=7r/2 vertical planes. 



maximum 
angle of 
for both 
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XIV. THE TERMINATED SLOPING DOUBLE RHOMBOID ANTENNA 



The orientation of the terminated sloping double rhomboid 
antenna is depicted on the next page in Figure 14.1, where 1, 
is the length of the four short antenna wires (A, D, E, and 
H) , 1 2 is the length of the four long antenna wires (B, C, F, 
and G) , 1 3 is the length of the radials which connect the feed 
to the two termination points, H 0 is the feed height above 
ground, H 3 is the height of the termination points above 
ground, y is the angle between by the radials which connect 
the feed to the termination points, 6 is the angle between the 
termination radials and the short feed wires, 77 is the angle 
between the termination radials and the long feed wires, a is 
the angle between the plane which contains the antenna and the 
x-y plane, R is the radial coordinate, 6 is the elevation 
coordinate, and <p is the azimuth coordinate. 

The Mathcad application for the double rhomboid antenna 
requires the following inputs: 



1 , length of the short antenna wires (meters) 

1 2 length of the long antenna wires (meters) 

1 3 length of the termination radials (meters) 

H 0 feed height above ground (meters) 

H 3 height of termination points above ground (meters) 

7 angle subtended by termination radials (degrees) 

f 5 operational frequency (Hertz) 

R distance from antenna (meters) 

e r relative dielectric constant of ground plane 

a conductivity of ground plane 

d elevation angle index (from Table 3.1) 

The feed height (H 0 ) and termination height (H 3 ) can be any 



value greater than zero. The length of the radials to the 
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TOP VIEW 




Sloping Double Rhomboid Orientation 



FIGURE 14.1: Spatial orientation of the terminated sloping 

double rhomboid antenna for its Mathcad application. 
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termination points must be less than the combined lengths of 
a long and short antenna wire (1 3 < 1, + 1 2 ) • The angle 77 is 
found from the law of sines and cosines by 



t| = s m 



-i< 



11 . 

ySin 

1 2 



cos 



J (i 1 ) 2 + (i 3 ) 2 -(i 2 ) 



(14.1) 



2l 1 l :} 

The angle 6 is also found from the law of cosines by 



6=1 



cos' 



(I 1 ) 2 +(l 3 ) 2 -(l 2 ) 

21 , 1 , 



- Y 



The angle ce is then found from H 0> 
angle between the wires (y/2) by 



(14.2) 

H 3 , 1,, 1 2 , and the half- 



^ 3—^0 

+6|+1 2 *cos 

The user inputs the frequency (f 5 ) for which the radiation 
parameters discussed in Chapter 3 are computed. Distance from 
the antenna (R) must meet the far- field requirements of 
Chapter 3.C, and the elevation angle index (d) sets the 6 
coordinate for which a horizontal radiation pattern is 
determined. Table 3.1 lists possible indices and their 
corresponding elevation angles from .285° to 88.857° in 
increments of about 2°. Indices between those listed can be 
used to interpolate a better estimate of a desired elevation. 

The radiated electric field for the double rhomboid is 
obtained in a manner analogous to that used for the vertical 
dipole discussed in Chapter 4. The electric field equations 
for the double rhomboid are given by [Ref 6 : p. 238-244] 



(H 



(14.3) 



a=sin 



-i 



1,' cos|-| 
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p I y 
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F e B = 30 J 0 -^-^- e ' jM2t/2 ( cos (a) cos ( 9 ) sin ( 0 +^ + 6 ) ( 14 . 12 ) 
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where U 1# U 2 , U 3 , U 4( U 5 , U 6 , U 7/ and U 8 are given by 



C7 1 =l-cos (7 X ) 
C7 3 = 1-cos (Y 3 ) 
C/ 5 = l-cos (T 5 ) 
C7 7 = 1-cos (Y 7 ) 



U 2 = 1-cos (Y 2 ) 

£7 4 = 1-cos(Y 4 ) 

C7 6 = 1-cos (Y 6 ) ' 

C7 8 = 1 -cos (T 8 ) where 



cos (T x ) =cos (0) sin (a) +sin(0) cos (a) sin(<J>-a-6) 
cos (T 2 ) =cos (0) sin (a) +sin (0) cos (a) sin (4>-a-r| ) 
cos (Y 3 ) =cos (0) sin (a) +sin(0) cos (a) sin(4>+a+r|) 
coc(T 4 ) =cos (0) sin(a) +sin(0) cos (a) sin(<j>+a+6) 
cos (Y 5 ) = sin (0) cos (a) sin(<J)-a-6) -cos (0) sin (a ) 
cos (T 6 ) =sin(0) cos (a ) sin (4»-a-r| ) -cos (0) sin (a) 
cos ( Y 7 ) =sin (0) cos (a ) sin (<|)+a+T) ) -cos (0) sin (a ) 
cos (T 8 ) = sin (0) cos (a) sin (<j>+a+6) -cos (0) sin (a) 



The double rhomboid's radiated electric field is a combination 
of 6 and $ components, and its radiation pattern varies with 
changes in <£. The total electric field is the vector sum of 
the d and $ components. The radiated fields of each antenna 
wire are considered separately, and superposition is used to 
determine the total electric field. Equations 14.4 and 14.5 
are the 6 (vertically polarized) and $ (horizontally 
polarized) components, respectively for wire A, and equations 
14.6 and 14.7 are the d and $ components, respectively for 
wire B. Equations 14.8 and 14.9 are the analogous equations 
for wire C, 14.10 and 14.11 those for wire D, 14.12 and 14.13 
those for wire E, 14.14 and 14.15 those for wire F, 14.16 and 
14.17 those for wire G, and 14.18 and 14.19 those for wire H. 
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The first two terms inside each set of brackets represent the 
space wave components and any additional terms . represent the 
surface wave components. I 0 is taken to be one since a 
sinusoidal current input with a maximum of unity is assumed. 

The requested inputs are used to calculate the following 
variables using a constant 0 of 0=0 and 0 = 7 r /2 for 632 discrete 
values of 6 which are equally incremented from -7r/2 to tt / 2 : 



X 5 wavelength of the operational frequency (meters) 

free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numeri cal distance (horizontal polarization) 

T v vertical reflection coefficient 

T h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
6 with 0=0 and 0 =tt/ 2. The total space wave and surface wave 
distributions are determined from vector addition of 
corresponding 6 and $ components. The space wave and surface 
wave results are then combined for corresponding discrete 
values of 6 to obtain the total radiated electric field 
distribution for the 0=0 and 0 = 7 r /2 vertical planes. The space 
wave, surface wave, and total electric field results are then 
normalized with respect to the maximum field intensity of 
each, and the normalized magnitudes are plotted for each 
discrete 6 to depict the radiation patterns. The Mathcad 
double rhomboid antenna application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and 0 = 7 r /2 vertical planes. 
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The variables corresponding to the selected elevation 
angle index (d) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 d ) as <f) varies from 0 to 2n in 632 
equal increments . The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

To find the contributions of the d and $ electric field 
components to the total average radiated power, P rad , equations 
3.8 and 3.9 are used to integrate equations 14.4 to 14.19 over 
the hemispherical Gaussian surface above the ground plane at 
a fixed radius (R) from the antenna. The sum of the integrals 
is the antenna's P rad . With the discrete values of the electric 
field and total average radiated power determined, the Mathcad 
application predicts the following radiation parameters using 
the equations from Chapter 3 : 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective isotropic radiated power (Watts) 

A^, maximum theoretical effective area (square meters) 

l mi maximum theoretical effective length (meters) 

P c numerical distance (vertical polarization, 0 = 90°) 

P m numerical distance (horizontal polarization, 0=90°) 



Angle^ .elevation angle of maximum directive gain (degrees) 
The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A^J , maximum 
theoretical effective length (l^J , and the elevation angle of 
maximum directive gain (Angle^J are all determined for both 
the <p=0 and <p=ir/2 vertical planes. 
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XV. THE VERTICALLY POLARIZED LOG- PERIODIC DIPOLE ARRAY 



The orientation of the vertically polarized log-periodic 
dipole array is depicted in Figure 15.1, where N is the number 
of dipoles, n is an index equal to 0,1,2..N-1, l n is the 
length of the n^+l dipole, d„ is the distance between the n^+l 
and n*+2 dipole elements, H 0 is the height of the antenna feed 
above ground, H„ is the height of the n^+l dipole center, SF is 
the angle between the vertical axis and the array's center 
axis, a, is the angle between the center axis and the line 
which connects the upper tips of the elements, a 2 is the angle 
between the center axis and the line which connects the lower 
tips of the elements, a 3 is the angle between the horizontal 
and the line connecting the lower tips of the elements, R is 
the radial coordinate, 0 is the elevation coordinate, and 0 is 
the azimuth coordinate. 

The Mathcad application for the vertically polarized log- 
periodic array requires the following inputs: 



N number of dipole elements 

1 0 length of the shortest dipole (meters) 

1, length of the second shortest dipole (meters) 

d 0 distance between first two elements (meters) 

rad 0 radius of the shortest dipole (meters) 

H 0 feed height above ground plane (meters) 

SF angle from vertical axis to array axis (degrees) 

f 5 operational frequency (Hertz) 

R distance from array (meters) 

e r relative dielectric constant of ground plane 

a conductivity of ground plane 

w elevation angle index (from Table 3.2) 

ADM transmission line characteristic admittance 

TIMP termination impedance connected to longest dipole 
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FIGURE 15.1s Spatial orientation of the vertically polarized 
log-periodic dipole array for its corresponding Mathcad 
application. 
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The user inputs the frequency for which the radiation 
parameters discussed in Chapter 3 are computed: The lengths 
of the first two elements (1 0 and 1,) are used to calculate the 



log-periodic relationship (r) given by [Ref 10: p. 317] 



The remaining element lengths (1 2 to 1 N .,) are then found from 



As noted in equation 15.1, indices start at zero instead of 1, 
so an index simply refers to the index-plus-one position in 
terms of successive values for a given parameter. This is to 
maintain continuity with Mathcad which also indexes from zero. 
The radius of the first element (rad 0 ) is used to find the 
remaining radii by substituting rad„ and rad^, for l n and l n+ , 
in equation 15.2. The angle between the vertical axis and the 
array axis OiO can be between 0° and 90°. The angles a 2 and 
a 3 are found by substituting 




(15.1) 



1 



1 



n 



(15.2) 






(15.3) 



into 



1-T 



(15.4) 



1 



n+1 



sin (a 2 + a 3 ) - tan (a 3 ) cos (a 2 + a 3 ) 
which yields 
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The angle a, is then found from [Ref 10: pp . 351-354] 

sin (o^) _ sin(a 2 ) (15 5) 

cos (a 1 +a 2 + a 3 ) cos(a 3 ) 

or 

J sin (a,) cos (a, + a,) 

a 1 =tan' 1 — 1 — — - — %— 3 

cos (a 3 ) +sm (a 2 ) sin (a 2 +a 3 ) 

The remaining separations (d, to d N . 2 ) are then found from 
equation 15.4. The dipole half-lengths (hj are found from 
l n /2 , the feed height (H 0 ) must be at least equal to the half- 
length of the longest dipole ( h N _, ) , and distance from the 
array (R) must meet the far-field requirements of Chapter 3.C. 
The elevation angle index (w) sets the d coordinate for which 
a horizontal radiation pattern is determined. Table 3.2 lists 
possible indices and their corresponding elevation angles from 
.57° to 89.1° in increments of about 2.3°. Indices between 
those listed can be used to interpolate a better approximation 
of a desired elevation. The termination impedance (TIMP) is 
assumed to be connected to the array's center axis opposite 
the feed at a distance of h N .,/2 from the longest element. 
Typical termination impedances are the transmission line's 
matched impedance (complex conjugate) or a short circuit 
(zero) . The characteristic admittance (ADM) is simply the 
inverse of the transmission line's characteristic impedance. 

If computed lengths, separations and radii are not 
desired, they can be entered manually as described on the page 
two of the application. Consult a Mathcad manual to be 
certain the entries are made correctly. 
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The open circuit impedance matrices must be calculated to 
determine the current distribution among the dipole elements. 
The matrix of open circuit mutual impedances between the 
dipole elements is represented as 





*o,o 


*0.1 


*0.2 


• * * *0, (N-l) 




*1.0 


*1.1 


*1.2 • 


’ ' • *1, (N-l) 


[Z] = 


*2.0 

• 


*2, 1 


*2, 2 


• ‘ * *2, (N-l) 




■ 

2 (N-l) , 0 


2 (N-l) ,1 


2 (N-l), 2 ' 


• ■ ■ • 

* ‘ ’ 2 (N-l) , (N-l). 



( 15 . 6 ) 



Each main diagonal term (Zj,) is the self - impedance of the i^+l 
dipole. The off-diagonal terms (Z ik ) are the mutual impedances 
between the i^+l and k^+l elements. Subscripts i and k are 
matrix indices from 0 , 1, 2 . . . .N-l . 

The matrix of open circuit mutual impedances between the 
actual dipoles and their images is represented as 





**0,0 


**o,i 


**0,2 • • 


• * **1(N-1) 




**1,0 


**1,1 


**1.2 • • 


• ' **1, (N-l) 


[ ZI ] = 


**2,0 


**2.1) 


**2 ,2 • • 


* * **2 , (N-l) 




22 (N-l) ,0 


22 (N-l) ,1 


22 (N-l), 2 * * 


• * 22 (N-l) , (N-l), 



( 15 . 7 ) 



The main diagonal terms (Zljj) are the mutual impedances 
between the i^+l element and its image, and the off-diagonal 
terms (ZI ik ) are the mutual impedances between the i^+l element 
and the image of the k^+l element. 
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The main diagonal self - impedance terms of equation 15.6 



( Z ; j) are found from [Ref 6: pp. 205-206] 



Z ^- l-costV) IOW-20W] 



(15.8) 



+ e J ' 2phi [i?( V 0 ) -2Q{V 1 ) ]+2[C ) ( U f 0 ) -Q( Uj -Q( Vj ]+2Q( U' 0 ) [1+cos (2 P^) ] ] 

X. X 

where Q(x) =Ci (x) -jSi {x) =f cos (y) dy-j-f ^ - in( ^ ’- ) -dy (15.9) 

« y o y 

U 0 = P [^2rad|+4A|-2/i i ] C7 / 0 = v /2prad i V 0 = p [^2rad|+4A|+2hj 

Ui=P [^radf+hf-ijJ V 1= p [ v /2rad?+i2?+i2 i ] 



For ^=90°, the off-diagonal mutual impedance terms of 
equation 15.6 (Z^) are given by [Ref 11] 



cos [p ] -cos [p (h^h^ ] 

[e" jP [Q ( C7 0 ) - Q ( ^ ) -f? ( U 2 ) ]+ e Jp {h ** hk) [Q ( V 0 ) - Q ( V, ) -Q ( V 2 ) ] 

+ e -jp <*i-A*> [p ( U' 0 -Q ( V x ) -Q ( U 2 ) ] + e Jp (hrhk) [QiV'^-QlUj-QiVj] 
+2Q{$D ik ) (cos [p ( h r h k ) ] +cos [p (h^h^ ] )] 



where Q{x) =Ci (x) -j'Si ( x ) = f 

oo 

C7 0 =P [jDl k+ (h i+ h k )*-(h i+ h k )] 

^ 0 =P [^i.* + (A i -A*) 2 -(A i -A*)] 
t^i = P 
d 2 =P 



cos IZl dy _ j.f sin (y) d 

y { y 

v 0 = p [^? Jt +(h i +A Jt ) 2 +(/i i +/i *) 

<=P [V^i.* + (Ai-A*) 2 + (Ai--ft*)' 
^i = P 

^ 2 = P 
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For ^ less than 90°, the dipoles are in echelon. The 

expression for the mutual impedance terms of equation 15.6 
(Z ik ) for an echelon configuration is given by [Ref 12] 

Z i ( J .=15{e - j * (Ai_SDi '* ) [Q ( U 0 ) -Q ( U x ) ]+ e jP (h ^ SD ^ [Q {V 0 ) - Q (VJ ] (15.11) 

+ e " Jp {h ^^[Q (u' 0 )-Q(U 2 )} + e jp {h ^ SD ^[Q (V / q )-Q{V 2 )] 

+ e " Jp (uj-QiUj ] + e ip (h r 2h *- SD i,J [£> ( V 3 ) -Q ( V 2 ) ] 

+ e ' ip (A i +2 ** +5D i.k> [p ( C/ 4 ) -£> ( t/ 2 ) ] + e jp (iJ i +2i3 * +SD i.* ) [j? ( V 4 ) -Q ( V 2 ) ] 

+2cos ( ph 2 .) e~ j * SDi -\Q (W 2 )-0(W 1 ) ]+2cos ( p^) e j * SDi -\Q (Y 2 )-Q(Y X )] 

+2cos(ph i )[e" J ' p<2 ^ +SDi -* ) [£)(^ 2 ) -Q(W 2 )]+e j * Uh ** SDi - k) [Q(Y 2 ) -Q(Y 3 ) ] ] } 






where Q(x)=Ci(x) 

oo ~ 

t/ 0 = P [v/z?? *+ iSD~^h~) 2 + (SDi't-h,) ] 



V, 



U‘ 






v. 



r ■ I y ± , x -L 

u = P [^i.* + 2 - ( SD i.k~ h i'> 

T ' o=P 

r/ o = P <^D~^hT) 2 + (5^,*+^)] 

i = P [V^iV TsD~^h^) 2 + (SD ik -h i +h k ) } 

-hi+h k ) 2 - (SD^-h^h*)] 



/i=p [^? *+ (SD^-^+A*) 2 - ( SD itk -h^h k ) ] 

u 2 = p <‘SD JiJe +A J +A*) 2 - 

v 2 " r ^ 7:13 T ^ '1 

^3 



v; 



L V — J- # ^ ' JL , A J. a - U / A. JL 

=p [^i iJt + (SD iik +h £ +h k ) 2 + (5D i(jt +h i+ ^)] 

= P <SD i( *-h i+ 2A*) 2 + (SD iik -h i+ 2h k ) ] 

= p [^1.*+ (SZ> iiJt -A i+ 2A*) 2 - (SD iik -h i +2h k ) ] 
r 4 = p \jDl k + (SD~^h~l 2h k ) 2 - (SD lik +h i + 2h k ] 
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V 4 = p [pi ~+ {SD~^ 2 ^) 2 + (SD irk+ h i+ 2 h^) ] 



^=P [plk+SDlt-SD^ 

^i=P [V^i.*+5D? >Jt +SD if J 
^ 2 = P [V^i.* + (SD ik +h k ) 2 - (SD i# *+**) ] 
y 2 =P [\/Di , ( SD i ik +h k ) 2 + (SD i>k+ h k )} 

^ 3 = P [Plk+ (SD i k+ 2h k ) 2 - (5D i>Jt +2A Jt )] 

F 3 = p [^(Dl,^(5D i(J , + 2i2 Jc ) 2 + (SD. i ^2V] 

The D ik terms in equations 15.10 and 15.11 are the perpen- 
dicular distances between the i^+l and k^+l elements given by 

^i.jc=^.i=E d n ‘ sin(T) (15.12) 

n=i 

The SD ik terms for equation 15.11 are the vertical distance 
from the center of the i^+l dipole to the z coordinate of the 
nearest end of the k^+l dipole and are given by 

k-1 

SD i k =SD kii =Y^ d n ' coa(V) -hj (15.13) 

n=i 

For ground planes with high conductivities relative to 
the operating frequency, it is necessary to compute equation 
15.7, [ZI] , for the mutual impedances between the elements and 
their images. The main diagonal terms < ZIj j) are the mutual 
impedances between the i^+l dipole and its own image. These 
Z j; terms are calculated from the mutual impedance equation for 
a collinear configuration given by [Ref 12] 
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(15.14) 



ZJ ifi =15<e' Jp(l,i_l,ii ' i) [C(C7 0 ) -Q(u 1 )]+e j * lhi ~ hii - i) In 



^i.i 



+ e 



^o) -j?( l^ 2 ) ]+e" jP(iJi+ijii ‘ i) ln — 



hii.i-hi i 



hi j + 2 A i 



hi i.i +h i 



+ e 



j'P (h]*hi J ' i ) 



[QiU^-QiU^] 



+ e 



-j'P (hj+hi If f) 



i#i In 






i 



,jK 3 ii 1 +u 1 , 1 )p| Vr ^ -j?(V r 2 )] +e ’ Jp(3ili+j5ij,i> ln 



_ hi i,i +h i 

hi iii ^2h i 



hl i.i +3h i j 



+2cos(Pi2 i ) e jPijii - i [C>(y 2 ) -0(l r 1 )]+2cos(pA i ) e'^^ln 



lil 1 



+ 2COS ( p/2 i ) e JP<2V*ii.i ) ^( yj -p(r 3 )j 
+2cos (P^) e' jP(2Ai+hJi - i> ln 



Ai i,i +2A i 

hii^+h^ 



where Q(x) =Ci (x) -j -Si ( x) =j cos izl cfy- j dy 

C7 0 =p [v/-Di,ic+ Ihi~]ph~) 2 + (hi ik -h i ) ] 

V 0 =P 

^o=P [^i,jc + (Ai^+Zij) 2 - (hi i k +hi) 

^o=P 

tfi-P [^i.* + Thi^P^h ) 2+ Mii.k-hi+h k ) } 

Vi=P 2 - (hi itk -h^h k ) } 

u 2 = P [^i,jc + {hi ~^^) 2 - {hi iik +h^h k ) ] 
v 2 =P [A\* + ihi~^h~^) 2 + (hi i>k +h i+ h k ) ] 

U 3 = P <Ai i<Jt -Ai + 2A*> 2 + ( hi itk -h i + 2h k ) ] 
v 3 =p [^i.*+ (hi^-hs 2 hj 2 - {hi ik -h i + 2 h k ) } 



ill 



(7 4 = p [/d?, *+ {hi~~^h~+ 2hJ 2 - (hi ik +h i +2h k ) } 



V 4 = P []fnTk+ t + *i +2 ^) 2+ ( h ii,k +h i +2h k) ] 

^i = P T-^p J 

w 2 = P [v^i,* + <^i,* + V 2 ~ (^ii.jt+A*)] 

Z 2 = P [jDl k *(hi itk *h k )^ (Ai i(Jt+ ^)] 

P/ 3 = p [^1,*+ UiT~^+2h^) 2 - (Aij^+2/ip] 

Z 3 = P [^ (£>?,*+ U2i~^2h^+ (hi itk +2h k )} 

The symbol hi i; in equation 15.14 is the height of the i^+l 
dipole above the nearest end of its own image, and with i=k is 
obtained from 



hi iik =H i +H k -h k (15.15) 

The off-diagonal mutual impedances between the dipoles and the 
images (ZI ik ) are found using equation 15.13, except the Q(x) 
arguments are those given for equation 15.14, but with SD ik 
replaced by hi ik from equation 15.15. 

The function Q(x) given by equation 15.9 is defined in 
terms of sine [Si(x)] and cosine [Ci (x) ] integrals. Mathcad 
is incapable of evaluating Si (x) and Ci(x) directly, so a 
series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q(x). The arguments for Q(x) are 
always real, so it is unnecessary to find Si(x) and Ci (x) 
expressions valid for complex numbers. For arguments less 
than one, Si(x) and Ci (x) are evaluated by series expansions 
given by [Ref 9: p. 232] 
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(15 . 16) 



si (x) =£ 

73=0 



(_1) Jl x 2ml 

(2n+l) (2n+l) ! 



Ci ( x ) =Y + ln (x) + J2 

13=1 



(-l) n x 2n 
( 2 n) ( 2 n) ! 



where y=0 . 5772156649 



(15.17) 



For arguments greater than one, Si(x) and Ci (x) are 
defined in terms of two auxiliary functions, f (x) and g (x) 
given by [Ref 9: p. 233] 



x° + a,x D +a,x*+a,x“+a, 



\ 



^ x Q +b 1 x 6 +b 2 x 4 +b 2 x 2 +b 4 

|e (x) | < 5-10- 7 



+e (x) 



fix) =A 

x 

where 

a 1 =38 . 027264 
a 2 = 265. 187033 
a 3 =335 .677320 
a 4 =38. 102495 



2^=40 . 021433 
jb 2 = 322. 624911 
jb 3 = 57 0 . 236280 
jb 4 = 157 . 105423 



9ix)=^~ 



x" 



x 8 +a 1 x 6 +a 2 x i +a 3 x 2 +a 4 

x 8 +b 1 x 6 +b 2 x i +b 3 x 2 +b 4 J 



+e (x) 



where 



a x = 42.242855 



|e (x) | < 3 ‘1 0 ' 7 
2^=48 . 196927 



a 2 = 302 .757865 



jb 2 =482 . 4859 84 



a 3 =352 . 018498 



jb 3 = 1114 .97 8885 



a 4 =21 . 821899 



2? 4 =449 .690326 



(15.18) 



(15.19) 



The sine and cosine integral functions are then defined in 
terms of these auxiliary functions by [Ref 9: p. 232] 



Si ix) = — -f(x) cos (x) -gix) sin(x) (15 . 20) 

2 

Ci {x) =f (x) sin (x) -gix) cos (x) (15 . 21) 
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In addition to the mutual impedance matrices, it is 
necessary to calculate the short circuit admittance matrix to 
determine the elemental current distribution. The short 
circuit admittance matrix is given by [Ref 10: p. 321] 



[Y] = 



Y 0 ,o Y 0 .i 



0 



*i.o Yi.i Y 1i2 



0 

0 

0 

0 



Y Y Y 

1 2 , 1 •‘ 2,2 *‘ 2,3 



0 0 
0 0 
0 



0 

0 



Y Y Y 
" I 3,2 J 3,3 -* 3,4 



0 

0 



0 

0 



Y(N-2),(N-2) Y( N _2),(n-2) Y(n-2) , (N-l) 

® Y (N-l) , (N-2) Y( N _ !) (n-x) 



(15.22) 



where the diagonal terms are given by 
F o, o=-i F o c ot (pcf 0 ) and (n-d = Y i t~ JY 0 cot 

and for 0 < n < N-l : 

Y„,„ Y {n _ 2)i w . 2 ) =-ir 0 (COt(pd J 1 .i) + COS(pd n ) 

and the off diagonal terms are given by 

Y (n-l) , ji = Y n ' (n-l)—~3 YqCSC ( $d n _i) 

The term Y 0 is the characteristic admittance of the trans- 
mission line given by the input (ADM), and Y' T is given by 



cos 



y r'V 






YqZjC OS 



(^1 



+jsin 



\ 



(15.23) 



‘w- 1 



where Z T is the termination impedance input (TIMP) . 

The ground conductivity relative to the operating 
frequency (a rel ) is then computed by [Ref 6: p. 640] 



°rej = 



we r e 0 



18000 q 

f-Mhz ®r 



(15.24) 
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to determine which equation to use for the array's current 
distribution. A surface with a relative conductivity greater 
than 20 is considered to be highly conductive relative to the 
transmitted frequency and the current distribution is given by 
[Ref 10: p. 342] 



where [U] is the N X N identity matrix. The entries of the 
[IB] matrix are the base currents of the dipole elements. The 
matrix [IREF] represents the input current to which the [IB] 
distribution is referenced. Since the current input is only 
to the feed at the center of the first element, and since the 
input current is sinusoidal with a maximum of unity, [IREF] is 
a 1 X N matrix with a one as the first entry and N-l zeroes in 
the remaining positions. 

For surfaces with a relative conductivity of less than 
20, the mutual impedances due to the images are neglected and 
the current distribution is given by 



The electric field equations are referenced to the 
current maxima values for the dipole elements, [I] , which are 
obtained from the [IB] matrix using [Ref 6: p. 209] 



[IB] ={[U]+[Y] [Z] + [Y] [ZI])~ 1 [IREF] 



(15.25) 



[IB] ={ [ U] + [ Y] [ Z ] T 1 [IREF] 



(15.26) 



Ji Sin(ph i ) 



IB, 



(15.27) 
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The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a singularity error in the mutual 
impedance calculations, and it is necessary to vary the 
frequency such that no element is an exact integer multiple of 
the wavelength. The change required is only one or two 
percent, and the predicted radiation parameters are still a 
good estimate of those for the original frequency. 

The electric field for the vertical log-periodic array is 
obtained for each individual dipole in a manner analogous to 
that used for the vertical dipole discussed in Chapter 4. The 
equations for the individual elements are combined into a 
single expression for the array's total electric field by the 
array factor. The array factor for the vertical log-periodic 
array is [Ref 6 p. 222] 

N- 1 

A+jS=^ J i e ipHiCOS(0) [cos [p^cos (9) ] -cos (p^)] • (15.28) 
i=0 



1 + e -j2PH/Cos(6) 


/ 

rvu-rv)F e 


[sin 2 (0) - v^ 2 -sin 2 (0) cos(0) Y| 






n 2 ))\ 



^ jpYj-sec (a 2 +cc 3 ) [cos (0) sin (a 2 +a 3 ) +sin (6) cos (a 2 +a 3 ) sin (4>) ] 



where Y ( is the y coordinate of the i^+1 dipole given by 




sin (Y) 



(15.29) 
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The first two terms inside the brackets of equation 15.28 
represent the space wave, the third term represents the 
surface wave, and I ; is the current term for each element based 
on a sinusoidal current input with a maximum of unity. The 
array factor is combined with the element factor to find the 
expression for the total radiated electric field given by 



E _ j6 0e'^ fi A+fB 
0 R sin(0) 



(15.30) 



The requested inputs are used to calculate the following 
variables using a constant 0 of 0= 0 and 0= tt/2 for 312 discrete 
values of 6 which are equally incremented from -tt/2 to tt/2: 



h ; half-length of the i^+l dipole (meters) 

Hj height of the i^+l dipole (meters) 

Y; y coordinate of the i^+l dipole 

Ii current for the i^+l dipole 

\ 5 wavelength of the operational frequency (meters) 

/? free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance for vertical polarization 

T v vertical reflection coefficient 

F e vertical surface wave attenuation factor 



The calculated variables are used to evaluate the far- 
field space wave and surface wave for the discrete values of 
0 with 0= 0 and 0=7r/2 . The space wave and surface wave results 
are combined for corresponding values of 6 to obtain the total 
radiated electric field distribution for the 0=0 and 0=7r/2 
vertical planes. The space wave, surface wave, and total 
electric field results are then normalized with respect to the 
maximum field intensity of each, and the normalized magnitudes 
are plotted for each discrete 6 to depict the radiation 
patterns. The Mathcad vertical log-periodic dipole array 
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application computes the space wave, surface wave, and total 
electric field radiation patterns and radiation. parameters in 
the 0=0 and 0 = tt/ 2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 W ) as 0 varies from 0 to 27r in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

Equations 3.8 and 3.9 are used to integrate equation 
15.30 over the hemispherical Gaussian surface above the ground 
plane at a fixed radius (R) from the array to find total 
average radiated power (P rad ) . With the discrete values of the 
electric field and total average radiated power determined, 
the Mathcad application predicts the following radiation 
characteristics from the equations in Chapter 3 . 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

A^ maximum theoretical effective area (square meters) 

l^ maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

P m numerical distance (horizontal polarization, 0=90°) 



Angle^ .elevation angle of maximum directive gain (degrees) 
The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A^) , maximum 
theoretical effective length (1^) , and the elevation angle of 
maximum directive gain (Angle max ) are all determined for both 
the 0=0 and 0= 7 r /2 vertical planes. 
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As an example, the Mathcad vertical log-periodic array 
application was executed with the following inputs: 



number of elements 
length of first dipole 
length of second dipole 
first separation distance 
radius of first element 
height of antenna feed 
frequency 

distance from the antenna 
relative dielectric constant 
ground conductivity 
characteristic admittance 
termination impedance 
angle from vertical to array 
elevation angle index 



12 

3.246 meters 
3.732 meters 
1.096 meters 
0.00325 meters 
8 . 0 meters 
10’10 6 Hertz 
3000 meters 
4 

1 - 10' 3 

1/450 Mhos 
0 Ohms 
axis 90° 

285 (-17°) 



Figures 15.2 through 15.7 are the space wave and surface wave 
radiation patterns in the 0=0 and 0=7r/2 vertical planes, and 
the designated horizontal plane for this example. 




pattern for 0=0 vertical 
plane . 




FIGURE 15.3: Vertical log- 
periodic surface wave 
radiation pattern for 0=0 
vertical plane. 
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pattern for <p=ir/2 vertical 
plane . 




vertical plane. 




pattern for horizontal plane. 




radiation pattern for 
horizontal plane. 
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The following radiation parameters were predicted by the 
Mathcad application for a sinusoidal current input of one Amp. 



Total power radiated (Watts) 




37.967 


Radiation resistance (Ohms) 




75.934 


Numerical distance (vertical) 




57.126 




<t >= 0 


d>= 7r/2 


Directivity 


3.415 


12.176 


EIRP (Watts) 


129.656 


462.272 


Max eff area (sq meters) 


244.578 


872.012 


Max eff length (meters) 


14.038 


26.506 


Angle mM (degrees) 


21.71 


21.14 



These results are consistent with expectations for this 
particular configuration. Appendix D contains computer 
hardcopies of additional example calculations for the vertical 
log-periodic dipole array and compares predicted radiation 
parameters to those expected based on previous calculations. 
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XVI. THE HORIZONTALLY POLARIZED LOG-PERIODIC DIPOLE ARRAY 



The orientation of the horizontally polarized log- 
periodic dipole array is depicted in Figure 16.1, where N is 
the number of dipoles, n is an index equal to 0,1,2..N-1, l n 
is the length of the n^+l dipole, d„ is the distance between 
the n^+l and ^+2 dipole elements, H 0 is the height of the 
antenna feed above ground, H„ is the height of the n^+l dipole, 
^ is the angle between the vertical axis and the array's 
center axis, a is the angle between the array's center axis 
and the lines which connect the tips of the elements, R is the 
radial coordinate, 6 is the elevation coordinate, and (p is the 
azimuth coordinate. 

The Mathcad application for the horizontally polarized 
log-periodic array requires the following inputs: 



N number of dipole elements 

1 0 length of the shortest dipole (meters) 

1, length of the second shortest dipole (meters) 

d 0 distance between first two elements (meters) 

rad 0 radius of the shortest dipole (meters) 

H 0 feed height above ground plane (meters) 

i' angle from vertical axis to array axis (degrees) 

f 5 operational frequency (Hertz) 

R distance from array (meters) 

e r relative dielectric constant of ground plane 

a conductivity of ground plane 

w elevation angle index (from Table 3.2) 

ADM transmission line characteristic admittance 

TIMP termination impedance connected to longest dipole 



The user inputs the frequency for which the radiation 
parameters discussed in Chapter 3 are computed. The lengths 



122 



Top View 





FIGURE 16.1: Spatial orientation of the horizontally 
polarized log-periodic dipole array for its corresponding 
Mathcad application. 
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of the first two elements (1 0 and lj) are used to calculate the 
log-periodic relationship (r) given by [Ref 10: p. 317] 

T = -r— 2- (73=0,1,2. . . .N-2,N-1) (16.1) 

■‘■n* 1 



The remaining element lengths (1 2 to 1 N .,) are then found from 



1 = — 
ntl x 



(16.2) 



As noted in equation 16.1, indices start at zero instead of 1, 
so the index simply refers to the index- plus -one position in 
terms of successive values for a given parameter. This is to 
maintain continuity with Mathcad which also indexes from zero. 
The dipole half-lengths (hj are found from l n /2 . The radius 
of the first element (rad 0 ) is used to find remaining radii by 
substituting rad„ and rad^, for l n and l n+I into equation 16.2. 
The angle between the vertical axis and array axis (^) can be 
between 0° and 90°, and a is found from [Ref 10: p. 319] 



— ~ = (1-t) cot (a) 

h n + 1 

which can be written as 



a = tan _1 



(1 -t )h 



n + 1 



(16.3) 



(16.4) 



The remaining separations ( d, to d N . 2 ) are found from equation 
16.3. Feed height (H 0 ) can be greater than or equal to zero. 
Distance from the array (R) must meet the far- field require- 
ments of Chapter 3.C. The elevation angle index (w) sets the 
0 coordinate for which a horizontal radiation pattern is 
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determined. Table 3.2 lists possible indices and their cor- 
responding elevation angles from .57° to 89.1°. in increments 
of about 2.3°. Indices between those listed can be used to 
interpolate a better estimate of a desired elevation. The 
termination impedance (TIMP) is assumed to be connected to the 
array's center axis opposite the feed at a distance of h^/2 
from the last element. Typical termination impedances are the 
transmission line's matched impedance (complex conjugate) or 
a short circuit (zero) . The characteristic admittance (ADM) 
is just the inverse of the transmission line's characteristic 
impedance. If computed lengths, separations and radii are not 
desired, they can be entered manually on page two of the 
application. Consult a Mathcad manual to be certain the 
entries are made correctly. 

The open circuit impedance matrices are calculated to 
determine the current distribution among the dipole elements. 
The matrix of mutual impedances between dipole elements is 
represented as 



Each main diagonal term (Z M ) is the self- impedance of the 
i^+l dipole. The off-diagonal terms (Z ik ) are the mutual 
impedances between the i^+l and k^+l elements. Subscripts i 
and k are matrix indices (0, 1,2 . . . .N-l) . 





[Z] 



(16.5) 



Z (N-1),0 Z {N- 1),1 Z (N- 1) ,2 



Z {N- 1) , (W-l), 
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The matrix of open circuit mutual impedances between the 



actual dipoles and their images is represented as 



I ZZ 0,0 


ZZ 0,1 


ZZ 0 , 2 ■ 


• • 1 ZI 0 . ( N-l ) 


ZZ 1,0 


ZZ 1,1 


ZZ 1 , 2 ' 


• 1 ‘ ZI 1, (.N-l) 


ZI (N- 1) ,0 


ZI (N- 1) ,1 


ZZ (N- 1) ,2 1 


• • 1 ZZ (N-1) , (W-l) 



The main diagonal terms (Zljj) are the mutual impedances 
between the i^+l element and its image, and the off-diagonal 
terms (ZI ik ) are the mutual impedances between the i^+l element 
and the image of the k^+l element. 

The main diagonal self -impedance terms of equation 16.5 
are found from [Ref 6: pp. 205-206] 



z. 60 

1,1 1-cos (20^) 



[« 






[£>(C/ 0 )-2£)(C7 1 )] 



(16.7) 



+e j2 * hi [Q{V 0 ) -2C)(U 1 )]+2[C)(C7 / 0 ) -£>(£7 1 )-£>(V 1 )]+2£>(tf / 0 )[ 1+cos (2ph i )] ] 



where Q(x) =Ci ( x ) -j-Si (x) = J £2^1x1 dy-j-J s ^ n ^ dy (16.8) 

oo y o y 

U 0 = p ^f2ra^f+4hj-2h^ U 1 0 =\/2§rad i V 0 = p [y2rad|+4h| + 2h i ] 

dj = p yjradf+hf-h^ V^-p [yiradj+hf+h^ 



The off-diagonal mutual impedance terms of equation 16.5 
are given by [Ref 11] 

z i.k z k.i cos [p (hj—hfc) ] —cos [p (hj+hj'.) ] (16.9) 

[e - J ’ P itls * hk) [QW^-QtUj-QUJjye™ (h ^ hk) [Q ( V 0 ) -Q ( V x ) -Q ( V 2 ) ] 

+ e ~ j * {hi ~ hk) [Q( U ' 0 -Q(V 1 )-Q( U2 )]+e j V (hi ~ h * ) [Q( V ' 0 ) -Qiuy -Q(V 2 )] 
+2£)(P D ik ) (cos [p (h-h*) ] +cos [p (h i +h^) ] )] 
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wheie Q(x) =Ci (x) -j-Si (x) 


Jcos(y) sin (y) 

l y { y 


^o=P [\/Dl^(h i+ h k ) 2 -(h 1 + h k )} 


^o = P [jDl k+ (h i +h k ) 2 + (h 1 + h k )] 


0=P (*!-*»)] 


v'o-P <*,-**> i+ (A,-** 


^i=P [iDlk+hl-h, 


Va=P [/Df.^hl+hi] 


^ 2 =P [jDik+hl-h^ 


^ 2 =P {plk+h 2 k+ h k } 



The D i>k terms are the distances between the i^+l and k^+l 
dipole elements and are given by 

*■-1 

D i.l c=%i = E d n (16.10) 

n~i 



For ground planes with high conductivities relative to 
the operating frequency, it is necessary to compute equation 
16.6, [ZI] , for the mutual impedances between the elements and 
their images. The mutual impedance terms (ZI i>k ) are given by 
equation 16.9 with the Q(x) arguments 



tf 0 =P 
U' 0 =P 



\Jdii,k + (hi + h k ) 2 - (h i +h k ) ] 
Jdih'+UirhJ'-Uii-hj] 
U x = P [)/dii, k +hj-h^ 

U 2 = P [jdif ik+ hZ-h k ] 



^o = P [/dill c + 2+ ] 

v'o-P [/dil k +(h i -h k ) 2 +(h i -h k ) 

Vi=P [/dil k +hf+h^\ 

^2 = P [jdii'k+hk+h)^ 



The di ik entries are the distance from the i^+l dipole to the 
image of the k*+l dipole and are given by 



di iik =/TH^Hj^ (D iik sin ( J ¥) ) 2 



(16 . 11) 
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The function Q(x) given by equation 16.8 is defined in 
terms of sine [Si(x)] and cosine [Ci(x)] integrals. Mathcad 
is incapable of evaluating Si(x) and Ci (x) directly, so a 
series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q (x) . The arguments for Q (x) are 
always real, so it is unnecessary to find Si (x) and Ci (x) 
expressions valid for complex numbers. 

For arguments greater than one, two auxiliary functions, 
f (x) and g(x) are determined by [Ref 9: p. 233] 




(16.12) 



where 



|e(x) | < 5-lCT 7 



a 1 = 38 . 027264 



£>3=40 .021433 



a 2 = 265 . 187033 



b 2 =322.624911 



a 3 = 335 .677320 



£>3 = 570.236280 



a 4 = 38 . 102495 



£> 4 =157 . 105423 




(16 . 13) 



where 



|e (x) | < 3-10' 7 



a 1 =42 . 242855 



£>3=48 . 196927 



a 2 = 302 .757865 



£> 2 =482.485984 



a 3 =352 . 018498 



£> 3 = 1114 .978885 



a 4 =21 . 821899 



£> 4 =449 .690326 
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and Si (x) and Ci (x) are defined in terms of these auxiliary 
functions by [Ref 9: p. 232] 



Si (x) =-2. -f(x) cos (x) -g(x) sin (x) 
Ci (x) =f(x) sin (x) -g(x) cos (x) 



(16 .14) 

(16.15) 



For arguments less than one. Si (x) 
evaluated by series expansions given by [Ref 



and Ci (x) 
9: p. 232] 



are 



Si (x) =£ 

il— 0 



(-1)” x 2ml 
(2n+l) (2n+l) ! 



Ci(x) =y + ln(x) +£ 

73=1 



(- 1) 71 x 2n 
(2 n) (2 n) ! 



where 7=0.5772156649 



(16 . 16) 
(16.17) 



It is also necessary to calculate the short circuit 
admittance matrix to find the elemental current distribution. 
The admittance matrix is given by [Ref 10: p. 321] 



[Y] = 



Y o.o Y o.i 0 

*1.0 * 1.1 *1.2 



0 

0 

0 

0 



Y Y Y 

1 2,1 1 2,2 1 2,2 



0 0 
0 0 
0 



* 3,2 * 3 , 3 * 3 , 4 



0 

0 



Y (N-2) , (N-2) Y (N-2) . (N-2) Y (N-2),(N-1) 



Y (N- 1) , (N-2) Y (N- 1) (W-l) 



00 

0 0 .. 0 

where the diagonal terms are given by 

- v/ 



(16 . 18) 



*o,o~ J 'i r o c ot ( Pc? 0 ) and Y (n-d , (n-d~ y t .7*o co *' ( P^n-i^ 
and for 0 < n < N-l : 

Y n.n Y (N-2 ) . ( AT-2) = - j *o ( ' co t ( P d n _i ) +CO ; S ( p d„) 

and the off diagonal terms are given by 

Y (n-l) ,n =Y n, (n-l) = ~J Y 0 CSC ^ P^n-l^ 
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The term Y 0 is the transmission line's characteristic admit- 
tance obtained from the input (ADM) , and Y' T is given by 



Y 1 =y ' ■ 
* t 1 o 



cos P 



h N - 1 



+j'F 0 Z T sin p 






YqZ ^ os p^lj+jsin^p^ij 



The ground conductivity relative 
frequency (a rel ) is then computed by [Ref 6 



(16.19) 

to the operating 
p. 640] 



q _ 18000 q 

coe r e 0 f Mhz € r 



(16 .20) 



A surface with a relative conductivity greater than 20 is 
considered to be highly conductive relative to the frequency 
and the current distribution is given by [Ref 10: p. 342] 

[IB] ={ [U] + [Y] [Z] + [Y] [ZI] } ’ [IREF] (16.21) 

where [U] is the N X N identity matrix. The entries of the 
[IB] matrix are the base currents of the dipole elements. The 
matrix [IREF] represents the input current to which the [IB] 
distribution is referenced. Since the current is input only 
to the feed at the center of the first element, and since the 
input current is sinusoidal with a maximum of unity, [IREF] is 
a 1 X N matrix with a one as the first entry and N-l zeroes in 
the remaining positions. 

For surfaces with a relative conductivity of less than 
20, the mutual impedances due to the image elements are 
neglected and the current distribution is given by 

[IB] ={ [U] + [Y] [Z] } ‘ [IREF] (16.22) 
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The electric field equations are referenced to the 
current maxima values for the dipole elements, [I] , which are 
obtained from the [IB] matrix using [Ref 6 p.209] 



1 sin ((Jhj) 



(16 . 23) 



The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a singularity error in the mutual 
impedance calculations, and it is necessary to vary the 
frequency such that no element is exactly an integer multiple 
of the wavelength. The change required is only a percent or 
two, and the predicted radiation parameters are still a good 
estimate of those for the original frequency. 

The electric field for the horizontal log-periodic array 
is obtained for each individual dipole in a manner analogous 
to that for the vertical dipole discussed in Chapter 4. The 
equations for the individual elements are combined into a 
single expression for the array's total electric field by the 
array factors. The array factors for the horizontal log- 
periodic array are [Ref 6: p. 214] 



N-l 



i= 0 



S 0 =^ J i [cos[Ph i (cos (0) cos (Y) +sin (0) sin (Y) sin (4>) )]-cos ( ph i ) ] 

(16 . 24) 

n*c os (0) 



1 + p e -j 2 $H i cos < 0 ) + d_r ) p e - j2 $ H i COB (6) /n 2 ~sin 2 (0) 

v v e n 2 c os(0) 



sm* 



n- 



jpV’yCse(T) [cos (0) cos (T) +sin (6) sin (T) sin (<J>) ] 
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N-l 

■T;[cos[P.h i (cos (0) cos (T) +sin (0) sin (Y) sin (<(>) )]-cos ( fihj) ] 

i = 0 

[l+r h e" i2pWiCOS(e) + (l-r^j) F m e -J' 2 P w i CO3 ( 0 )]. (16.25) 

e jpy i c , sc(T) [cos (0) cos (Y) +sin (0) sin (Y) sin (<t>) 3 



where Y; is the y coordinate of the i^+l dipole given by 



Yi = 



1-1 



E d r 



\n=0 



sin (Y) 



(16 . 26) 



The first two terms inside the brackets of equations 
16.24 and 16.25 represent the space wave, the third term 
represents the surface wave, and I- t is the current term for 
each element based on a sinusoidal current input with a 
maximum of unity. The array factors are combined with the 
element factors to find the total radiated electric field 
equation, and the expression for the radiated electric field 
distribution of the horizontal log-periodic dipole array is 
given by [Ref 6 : p. 213] 



E _ -j'60e cos ( 4 >) cos ( 0 ) s 

R 1 -sin 2 (0) cos 2 (<J>) e 

E _ j60e- J ’ p * sin (<j>) ~ 

* R 1 -sin 2 ( 0 ) cos 2 ( 0 ) * 



(16.27) 



(16.28) 



The requested inputs are used to calculate the following 
variables using a constant 0 of 0=0 and 0 = 7 r /2 for 312 discrete 
values of 6 which are equally incremented from - 7 r /2 to ir/2: 
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hj half-length of the i^+l dipole (meters) 

Hj height of the i^+l dipole (meters) 

Yj y coordinate of the i^+l dipole 

Ij current for the i^+l dipole 

X 5 wavelength of the operational frequency (meters) 

/3 free space wavenumber for operational frequency 

n index of refraction 

P e complex numerical distance (vertical polarization) 

P m complex numerical distance (horizontal polarization) 

r v vertical reflection coefficient 

r h horizontal reflection coefficient 

F e vertical surface wave attenuation factor 

F m horizontal surface wave attenuation factor 



The 


calculated 


variables 


are 


used to 


evaluate 


the 


radiated 


far- field 


space wave 


and 


surface 


wave for 


the 



discrete values of 6 with 0=0 and 0=7r/2. The space wave and 
surface wave results are combined for corresponding values of 
8 to obtain the total radiated electric field distribution for 
the 0=0 and 0=7r/2 vertical planes. The space wave, surface 
wave, and total electric field results are then normalized 
with respect to the maximum field intensity of each, and the 
normalized magnitudes are plotted for each discrete 6 to 
depict the radiation patterns. The Mathcad horizontal log- 
periodic dipole array application computes the space wave, 
surface wave, and total electric field radiation patterns and 
radiation parameters in the 0=0 and 0=7r/2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 W ) as 0 varies from 0 to 2tt in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 
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Equations 3.8 and 3.9 are used to integrate equations 
16.27 and 16.28 over the hemispherical Gaussian' surface above 
the ground plane at a fixed radius (R) from the array to find 
total average radiated power (P rad ) . With the discrete values 
of the electric field and total average radiated power 
determined, the Mathcad application predicts the following 
radiation characteristics from the equations in Chapter 3: 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

A^ maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0=90°) 

P m numerical distance (horizontal polarization, 6 = 90°) 



Ang 1 ©max .elevation angle of maximum directive gain (degrees) 
The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A max ) , maximum 
theoretical effective length (1^) , and elevation angle of the 
directivity (Angle max ) are determined for both vertical planes. 

As an example, the Mathcad horizontal log-periodic array 
application was executed with the following inputs: 



number of elements 
length of first dipole 
length of second dipole 
first separation distance 
radius of first element 
height of antenna feed 
frequency 

distance from the antenna 
relative dielectric constant 
ground conductivity 
characteristic admittance 
termination impedance 
angle from vertical to array axis 
elevation angle index 



12 

3.246 meters 
3.732 meters 
1.096 meters 
0.00325 meters 
8 . 0 meters 
20 * 10 6 Hertz 
3000 meters 
4 

5 • 10 3 
1/450 

0 Ohms 
90° 

285 (-17°) 
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Figures 16.2 through 16.7 are the space wave and surface wave 
radiation patterns in the 0=0 and 0=ir/2 vertical planes, and 
the designated horizontal plane for this example. 




pattern for 0=0 vertical 
plane. 



FIGURE 16.3: Horizontal log- 
periodic surface wave 
radiation pattern for 0=0 
vertical plane. 






pattern for 0= 7r/2 vertical radiation pattern for 0=7r/2 
plane. vertical plane. 
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pattern for horizontal plane. 




radiation pattern 
horizontal plane. 



for 



The following radiation parameters were predicted by 
the Mathcad application for a sinusoidal current input of one 
Amp. The predicted radiation parameters are: 



Total power radiated (Watts) 
Radiation resistance (Ohms) 
Numerical distance (vertical) 
Numerical distance (horizontal) 



< 6»0 



Directivity 
EIRP (Watts) 

Max eff area (sq meters) 
Max eff length (meters) 
A 11 ?!®®** (degrees) 



0.838 
63.301 
15 . 003 
4.904 
89 . 7 



75.546 

151.092 

93.742 

3398.15 



<6-tt/2 

18.384 

1388.82 

329.159 

22.971 

24.0 



These results are consistent with expectations for this 
particular configuration. Appendix E contains computer 
hardcopies of additional example calculations for the 
horizontal log -periodic dipole array and compares predicted 
radiation parameters to those expected based on previous 
calculations . 
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XVII. THE HORIZONTAL YAGI-UDA ARRAY 



The orientation of the horizontal Yagi-Uda array is 
depicted in Figure 17.1, where 1 NR is the length of the 
current driven element, NR is the number of parasitic 
reflector elements, ND is the number of parasitic director 
elements, n is an index equal to 0, 1, 2 . .NR+ND, and l n is the 
length of the n^+l element starting with the outermost 
reflector and counting consecutively toward the opposite end. 
The d„ terms are the separation distances between the n^+l and 
^+2 elements, H 0 is the height of all antenna elements above 
ground, R is the radial coordinate, 6 is the elevation 
coordinate, and 4> is the azimuth coordinate. 

The Mathcad application for the horizontal Yagi-Uda array 

requires the following inputs: 

NR number of reflector elements 

ND number of director elements 

l n length of the n^+l element (meters) 

d,, distance between n^+l and ^+2 elements (meters) 

rad„ radius of the n^+l element (meters) 

H 0 antenna height above ground plane (meters) 

f s operational frequency (Hertz) 

R distance from array (meters) 

e T relative dielectric constant of ground plane 

a conductivity of ground plane 

w elevation angle index (from Table 3.2) 

Unlike the log-periodic applications, all element 
lengths, radii, and separation distances must be input 
manually on page two of the application. Consult a Mathcad 
manual to be certain the entries are made correctly. The user 
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array for its corresponding Mathcad application. 
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inputs the frequency for which the radiation parameters 
discussed in Chapter 3 are calculated. As noted in the first 
paragraph, indices start at zero instead of 1, so the index 
simply refers to the index-plus -one position in terms of 
successive values for a given parameter. This is to maintain 
continuity with Mathcad which also indexes from zero. The 
dipole half-lengths (hj are found from l n /2. Antenna height 
(H 0 ) can be any value greater than or equal to zero, and the 
distance from the array (R) must meet the far- field 
requirements of Chapter 3.C. The elevation angle index (w) 
sets the 6 coordinate for which a horizontal radiation pattern 
is determined. Table 3.2 lists possible indices and their 
corresponding elevation angles from .57° to 89.1° in 
increments of about 2.3°. Indices between those listed can be 
used to interpolate a better approximation of a desired 
elevation. 

The open circuit impedance matrices are calculated to 
determine the current distribution among the dipole elements. 
The matrix of mutual impedances between the dipole elements is 
represented as 





N 

o 

o 


* 0.1 


^ 0,2 




* 0 . 


(JV-1) 




* 1 , 0 


* 1.1 


* 1.2 






(N-l) 


[Z] = 


^ 2,0 


^ 2,1 


to 




2 2 . 


(N-l) 




• 


. 


. 






. 



(17 .1) 



'( W - 1),0 ^( N - 1),1 ^( N - 1 ) ,2 



z 



(N-l) , (N- 1) 
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Each main diagonal term ( Z ; ; ) is the self impedance of the i^+l 
element. The off-diagonal terms (Z ijk ) are the mutual 
impedances between the i^+l and k^+l elements. Subscripts i 
and k are matrix indices (0 , 1 , 2 . . . .NR+ND) . 

The matrix of open circuit mutual impedances between the 
actual elements and their images is represented as 





' ZIo.o 


ZI o.i 


ZIo.2 ■ • 


• • 2 J 0 . (#_!) 




ZZ i,o 


ZI i.i 


ZI ll2 . . 


• • ZZ 1, (N-l) 


[ZI] = 


ZI 2 , o 


ZI 2,1 


ZI 2 ,2 ■ • 


■ • ZI 2, (N- 1) 




ZI(N-l) , 0 


ZZ (N- 1),1 ZZ (N- 1),2 1 • 


• * ZZ (N-1) , (N-l). 



The main diagonal terms (ZI^) are the mutual impedances 
between the i^+l element and its image, and the off-diagonal 
terms (ZI ik ) are the mutual impedances between the i^+l element 
and the image of the k^+l element. 

The main diagonal self - impedance terms of equation 17.1 
are found from [Ref 6: pp. 205-206] 



Z ± N \e~ j2?hi [J?(t/ 0 )-2!?(y 1 )l (17.3) 

1,1 1-COS (2ph i ) L L 0 1 J 

+ e j2 * b ‘ [Q{V 0 ) -2Q(V 1 )]+2[Q(U , 0 ) -Q(U 1 ) -Q{ V 1 )]+2£)(C/ / 0 )[1-hcos (2ph i ) ] ] 

X X 

where Q(x) =Ci (x) -j-Si ( x ) = j dy-j-f dy (17.4) 



£7 0 = p [ v /2rad J 2 +4h 1 2 -2h i ] U , 0 =s/2^rad i V 0 = p [^2rad 2 +4h?+2h i 

C/^P [^2 Tadf+hj-h i ] V^-p y2radJ7h?+h^ 
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The off-diagonal mutual impedance terms of equation 17.1 
are given by [Ref 11] 

Z ilk =Z kii = — [p ihi _ hk) ] 6 -° C os [p (h i+ h k ) ) ' (17 - 5) 

[ e -jP(V**)[ 1? ( i7o ) -q(u i ) -Q(U 2 )] + e jfiihi * hk) [Q(V 0 ) -Q(V X ) -Q(V 2 )] 

+ e -n [ q(u , 0 -q(v 1 )-q(u 2 ) ]+ e {hi ' hk) [Q ( v' 0 ) -Q ( U x ) -Q ( v 2 ) ] 
+2£)(pD i #Jt )(cos [p (h r h k ) } +cos [ P (2^+A*) ])] 



X X 



where Q(x) =Ci (x) -j-Si (x) 


-f c°s <y) d y- j-f sin ( y'dy 
« y o y 


U 0 = P yDl k +(hsh k )*-(h i+ h k )] 


V 0 = P [^i,Jc + (*i + ^) 2+ (*i-*A) ] 


U' 0 =P [ v /l3? Jt +(h i -h Jt ) 2 -(h i -h Jc ) 


v'„=|) '♦(V-ht)] 


^i=P [jDlk+hf-hi] 


y i=P [jVi.k+hj+h^ 


U 2 = P [jDl k ,ht-h k 


V 2 — p [v/d? k +h£+h k ] 



The D ik terms are the distances between the i^+l and k^+1 
dipole elements and are given by 

Jc-l 

D i. k =D kii =Y, d n (17.6) 

n=i 

The mutual impedances between the actual elements and 
their images must also be computed to determine the elemental 
current distribution. The mutual impedance terms of equation 
17.2 ( Zlj k ) are given by equation 17.5, except the D i>k terms in 
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the Q(x) arguments are replaced by di ik terms which represent 
the distance from the i^+l element to the image of the k^+l 
element and are given by [Ref 6: pp. 205-208] 

di itk =^H^D 2 iik (17.7) 



The function Q(x) given by equations 17.4 is defined in 
terms of sine [Si(x)] and cosine [Ci(x)] integrals. Mathcad 
is incapable of evaluating Si (x) and Ci (x) directly, so a 
series expansion and polynomial approximation are written into 
the Mathcad code for both the sine and cosine integrals to 
evaluate the functions Q (x) . The arguments for Q (x) are 
always real, so it is unnecessary to find Si (x) and Ci (x) 
expressions valid for complex numbers. 

For arguments less than one, Si (x) and Ci(x) are 
evaluated by series expansions given by [Ref 9: p. 232] 



Si (x) =£ 



(-l) n x 2n+1 



n=0 (2n+l) (2/1+1) ! 

(-l) n x 2n 



Ci ( x ) =y + ln (x) +£ 



(2/3) (2 n) ! 



(17.8) 

(17.9) 



where y=0 . 5772156649 



For arguments greater than one, Si (x) and Ci (x) are 
expressed in terms of auxiliary functions, f (x) and g (x) , 
given by Ref 9: p.232] 



Si (x) = — -f(x) cos (x) -g(x) sin(x) 
2 



(17 . 10) 



Ci (x) -f (x) sin (x) -g(x) cos (x) 



(17 .11) 
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The auxiliary functions f (x) and g(x) are evaluated by the 
polynomial approximations [Ref 9: p. 233] 



f U)=- 

x 



x 8 +a 1 x 6 +a 2 x 4 +a 3 x 2 +a 4 N 



k x 6 +b 1 x 6 +b 2 x A +b 2 x 2 +b 4 j 
where |e(x) | < 5-10' 7 



+e (x) 



a x =38. 027264 
a 2 =265 . 187033 
a 3 =335. 677320 



a 4 =38. 102495 



^=40 . 021433 
b 2 = 3 22 .624911 
jb,= 570 . 236280 



jb 4 = 157 . 105423 



g{x) =— 






x°+a 1 x D +a,x*+a,x‘ ! +a J , 



k x 8 +jb 1 x 6 +b 2 x 4 +b 2 x 2 +b 4 



+e (x) 



where 



a x = 42 . 242855 
a 2 -302 .757865 
a 3 = 352 . 018498 



a 4 =21 . 821899 



|e (x) | < 3 *10' 7 
^=48 . 196927 
i? 2 =482 . 4859 84 



Jb 3 = 1114. 978885 
jb 4 =449 .690326 



(17 . 12) 



(17 . 13) 



When the mutual impedances are found, the array's base 
current distribution is computed from [Ref 10 p. 258-259] 

[JB] ={[Z] +T h ' 0 [ZI]y 1 [VREF ’] (17.14) 

where 



r =r 

1 h,0 a A 



6 = 0 ' 



1 — n 

1+n 



and n is the complex index of refraction. The entries of the 
[IB] matrix are the base currents of the dipole elements, and 
[VREF] is the voltage matrix to which the [IB] distribution is 
referenced. Since the Yagi-Uda has only one driven element 
(corresponding to the NR* index) , and since a sinusoidal 
voltage response with a maximum of unity is assumed across the 
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input terminals, [VREF] is a 1 X NR+ND matrix with a 1 as the 
NR* entry and NR+ND- 1 zeroes in the remaining positions. The 
electric field equations are referenced to the current maxima 
values for the dipole elements, [I] , which are determined from 
the [IB] matrix using [Ref 6 p. 209] 



j.= -Bi 

J sin (pAj) 



(17.15) 



The mutual impedance calculations are not valid when any 
element length is an exact integer multiple of the wavelength. 
When this occurs, there is a singularity error in the mutual 
impedance calculations, and it is necessary to vary the 
frequency such that no element is exactly an integer multiple 
of the wavelength. The change required is only a percent or 
two, and the predicted radiation parameters are still a good 
estimate of those for the original frequency. 

The electric field for the horizontal Yagi-Uda array is 
obtained for each element in a manner analogous to that used 
for the vertical dipole in Chapter 4. The equations for the 
individual elements are simplified into two expressions for 
the array's electric field components by the array factor. 
The Yagi-Uda' s array factor is given by [Ref 6: p. 204] 

NR+ND- 1 

S = £ J i [cos[ph i (sin(0) sin(<f>))]-cos (p^)] e J ' pyiSin<0)sin<4,) (17.16) 

2 = 0 

where Y ± is the y coordinate for the i th + 1 element given by 

^i = E d n (17.17) 

n=0 
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The array factor is multiplied by the element factors to 
find the equations for the total radiated electric field 
distribution of the horizontal Yagi-Uda array, 



E e =-j6 0 



e cos (<t>) 



jPH 0 gos (0) 



(S)- (17.18) 



R l-sin 2 (0) cos 2 (<J>) 

[cos (0) (l-r v e- J ' 2PH « cos(e) ) + (l-r v ) F e e - j2 » H ° cos(e) • 



y(n 2 -sin 2 (8) 

,2 



E^=j60 



1 V 

e -m 



sin ^ e )-^ 2 - 31 ” 2191 

J2 2 



sin (<|>) 



R 1-sin 2 (0) cos 2 (0) 



cos (0) 
iPHoCOS (e> (S') • (17.19) 



[l+r A e~ j2pH ° cos(0> + (l-r ft ) F 6 _ ^' 2 P w o cos (0) j 



The first two terms inside the brackets of each equation are 
the space wave, the third terms are the surface wave, and I, 
is the current term for each element based on a sinusoidal 
voltage across the input terminals with a maximum of unity. 

The requested inputs are used to calculate the following 
variables using a constant <j) of <£= 0 and <£=tt/ 2 for 312 discrete 
values of 6 which are equally incremented from -7r/2 to tt/2: 



h 

Y 

I 



X 5 




half-length of the i^+l element (meters) 
y coordinate of the i^+l element 
current for the i^+l element 

wavelength of the operational frequency (meters) 
free space wavenumber for operational frequency 
index of refraction 

complex numerical distance (vertical polarization) 
complex numerical distance (horizontal polarization) 
vertical reflection coefficient 
horizontal reflection coefficient 
vertical surface wave attenuation factor 
horizontal surface wave attenuation factor 
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The calculated variables are used to evaluate the far- 



field space wave and surface wave for the discrete values of 
6 with 0=0 and 0=7t/2. The space wave and surface wave results 
are combined for corresponding values of 6 to obtain the total 
radiated electric field distribution for the 0=0 and 0=7r/2 
vertical planes. The space wave, surface wave, and total 
electric field results are then normalized with respect to the 
maximum field intensity of each, and the normalized magnitudes 
are plotted for each discrete 0 to depict the radiation 
patterns. The Mathcad horizontal Yagi-Uda array application 
computes the space wave, surface wave, and total electric 
field radiation patterns and radiation parameters in the 0=0 
and 0=7T / 2 vertical planes. 

The variables corresponding to the selected elevation 
angle index (w) are used to evaluate the radiated electric 
field components for the space wave and surface wave at the 
fixed elevation angle (0 W ) as 0 varies from 0 to 2ir in 312 
equal increments. The horizontal radiation patterns are then 
plotted for the space wave, surface wave, and total radiated 
electric field just as those for the vertical planes. 

Equations 3.8 and 3.9 are used to integrate equations 
17.18 and 17.19 over the hemispherical Gaussian surface above 
the ground plane at a fixed radius (R) from the array to find 
total average radiated power (P rad ) . Since the Yagi-Uda current 
distribution is referenced to a sinusoidal voltage input with 
a maximum of unity, the radiation resistance can not be 
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calculated from equation 3.10 as in previous applications. 
Instead, radiation resistance must be found with respect to 
the unity voltage input using 



zad 



L5dE 



2 R 



or 



R 



zad 



=_L5i 



zad 



2P 



zad 



(17 . 20) 



With the discrete values of the electric field and total 
average radiated power, the Mathcad application predicts the 
following radiation characteristics from Chapter 3: 



R rad radiation resistance (Ohms) 

D 0 directivity 

EIRP effective radiated isotropic power (Watts) 

maximum theoretical effective area (square meters) 

l max maximum theoretical effective length (meters) 

P e numerical distance (vertical polarization, 0 = 90°) 

P m numerical distance (horizontal polarization, 0 = 90°) 



Angle max .elevation angle of maximum directive gain (degrees) 



The directivity (D 0 ) , effective isotropic radiated power 
(EIRP) , maximum theoretical effective area (A^) , maximum 
theoretical effective length (1^) , and elevation angle of the 
directivity (Angle max ) are determined for both vertical planes. 

As an example, the Mathcad horizontal Yagi-Uda array 
application was executed with the following inputs: 



number of reflectors 
number of directors 
element lengths 
element separations 
element radii 0.001, 0 

height of antenna 
frequency 

distance from the antenna 
relative dielectric constant 
ground conductivity 
elevation angle index 



1 

1 

15.6, 15, and 14 
7 . 5 and 6 . 5 
0.001, and 0.001 
8.0 

10 • 10 6 
3000 
72 
4 

285 



meters 

meters 

meters 

meters 

Hertz 

meters 



(- 17 °) 
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Figures 17.2 through 17.7 are the space wave and surface wave 
radiation patterns in the 0=0 and 0=7r/2 vertical planes, and 
the designated horizontal plane for this example. 




pattern for 0=0 vertical 
plane. 




pattern for 0=7r/2 vertical 
plane . 




radiation pattern for 0=0 
vertical plane. 




FIGURE 17.5: Horizontal 
Yagi-Uda surface wave 
radiation pattern for 0=7r/2 
vertical plane. 
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FIGURE 17.6: Horizontal 
Yagi-Uda space wave radiation 
pattern for horizontal plane. 




FIGURE 17.7: Horizontal 
Yagi-Uda surface wave 
radiation pattern for 
horizontal plane. 



The predicted radiation parameters based on a one volt 
response across the input terminals are: 



Total power radiated (Watts) 
Radiation resistance (Ohms) 
Numerical distance (vertical) 
Numerical distance (horizontal) 



0.00482 

103.760 

0.04363 

2.26-10 6 



< 6-0 



(6=7t/ 2 



Directivity 1.679 

EIRP (Watts) 0.0081 

Max eff area (sq meters) 120.26 

Max eff length (meters) 11.506 

Angl^ (degrees) 89.7 



14.631 

0.0705 

1047.9 

33.966 

41.714 



These results are consistent with expectations for this 
particular configuration. Appendix F contains computer 
hardcopies of additional example calculations for the 



horizontal Yagi-Uda array and compares predicted radiation 
parameters to those expected based on previous calculations. 
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XVIII. REMARKS AND CONCLUSION 



This thesis is the culmination of nine months of research 
and computer programming. No new electric field equations 
were derived exclusively for this project. It was directed by 
NAVMARINTCEN prior to commencing that existing electric field 
equations would provide the basis for the radiation parameter 
predictions. The expressions used in this thesis are all 
previously derived by pioneers in antenna radiation theory 



such as Sommerfeld, Norton, Cox, King, 


Ma, 


and Walters. 


The 


equations 


are 


compiled here solely 


for 


the purpose 


of 


crediting 


the 


source references and 


to 


describe how 


the 



Mathcad code was assembled. 

The accuracy of the predicted radiation characteristics 
is totally dependent upon the extent to which the equations 
used realistically model the actual radiated electric fields 
of the antennas. It is near impossible to obtain analytic 
results which accurately model an antenna's radiation 
characteristics in all cases. There are too many operational 
and environmental variables to obtain one general, all-purpose 
expression. Even if all possible variables could be accounted 
for in a single expression, computer processing times would be 
unacceptably long. The equations used herein are simple 
enough to be evaluated by Mathcad in a reasonable length of 
time, but the radiation parameter predictions are accurate 
enough to be useful for antenna analysis. 
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As addressed in appendices A-F, the equations presented 
in previous chapters and evaluated by the Mathcad applications 
provide radiation parameters consistent with expectations for 
the inputs which have been executed to date. All antenna 
configurations which have been computed exhibit radiation 
characteristics consistent with other computational programs 
and empirically obtained patterns and parameters. There is no 
way to test the accuracy of the Mathcad predictions other than 
empirical measurements for each antenna configuration, clearly 
a task beyond the scope and purpose of this report. However, 
adequate analysis of the Mathcad results has been provided to 
demonstrate that the applications provide very good estimates 
of antenna radiation parameters as a base level analysis tool. 

Most of the applications have computations times under 
ten minutes. However, even with the simplified equations used 
by Mathcad, some of the applications can take upwards of two 
hours to compute on a 33 MHz 80386 PC. The applications for 
the rhombic, double rhomboid, and array (with more than four 
to five elements) antennas can take an exorbitantly long time 
for power calculations. For this reason, the applications 
should be run on the fastest PC available; a 50 MHz 80486 is 
preferable. Computation times can also be reduced by 
computing only the radiation patterns until it is determined 
which frequencies are of greatest interest. Then the average 
power calculations for those frequencies can be executed. 
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The Mathcad applications were written for NAVMARINTCEN to 
provide antenna analysts with the capability to predict 
radiation parameters based solely on antenna physical 
dimensions and ground properties. The Mathcad results 
computed thus far are consistent with expectations and are 
most likely providing accurate predictions of actual antenna 
radiation characteristics. Only extended use of the 
applications and empirical confirmation of results will prove 
out the accuracy of the Mathcad code . Results obtained thus 
far certainly justify continued use of the Mathcad code as an 
antenna analysis tool . 
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APPENDIX A; 

VERTICAL DIPOLE ARRAY COMPUTER OUTPUT 

This appendix contains computer hardcopies from the Mathcad 
vertical dipole application which show the input values and 
predicted radiation characteristics for two sample calculations. 
The configuration is a half-wave dipole at one quarter wavelength 
above soil (e r =lO and cr=10' 2 ) for the first example and the same 
configuration above seawater (e r =72 and cr=4) for the second. 
Reference 6 [p. 91] provides the radiation patterns and gain 
predictions from several sources for the configuration in the 
first example. 

The elevation of maximum directive gain is slightly higher 
for the Mathcad output than for those given in reference 6, but 
the overall radiation patterns are very similar. The 4.17 (6.2 
dB) value of directivity for Mathcad is quite a bit higher than 
the maximum gain values of about zero dB in reference 6 . 
However, since a half-wave dipole has a free- space directivity of 
1.64 (2.1 dB) , one would expect the actual directivity to be 
closer to the Mathcad prediction because of the effect of the 
reflected wave (constructive and destructive interference) on the 
space wave, and the lower total average radiated power resulting 
from ground plane losses. The seawater example yields results 
consistent with expectations. With respect to the soil example, 
the seawater example's directivity is slightly higher due to a 
stronger reflected wave, and the surface wave is stronger at 
grazing angles (0*=9O°) due to higher conductivity. 
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VERTICAL DIPOLE 



This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted vertically along the z-axis in a rectangular coordinate system with the feed 
at the antenna's center at a set height above the surface. Required inputs are the 
antenna length, feed height above the surface, transmitted frequency, distance from 
the antenna, the conductivity and dielectric constant of the surface below the 
antenna. The planar earth model is assumed in predicting radiation patterns. 
Predicted operating frequencies assume that the antenna is a quarter- wavelength, 
half-wavelength, three-quarter- wavelength, or full-wavelength dipole. A sinusoidal 
current input with a maximum of unity is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna. Plotted radiation patterns are valid for any vertical plane passing through 
the antenna axis. The electric field magnitudes to which the Radiation Patterns are 
normalized are displayed below their respective plots. Polarization is vertical for all 
vertical dipoles. 

Input the Dipole j ^ 

length in meters 






Wavelengths and Frequencies 



= 41 


*2 = 21 


x, , ±1 

3 3 


X 4 := 1 


'« = k 


f 2 := — 

X 2 


f 3 = “ 

h 




Possible Operating Frequencies (Hertz) 






^ Dipole: 
4 


f t = 1.5- 10 7 


^ Dipole: 


f 2 = 3* 10 7 


Dipole: 

4 


f 3 = 4.5- 10 7 


A Dipole: 


f 4 = 6* 10 7 


Input the operating 
Frequency (Hertz) 


f 5 := 30 10 S 


Input the Distance 
from the Antenna 
(meters) 


R 3000 


Input the ground 
Dielectric Constant 


er = 10 


Input the ground 
Conductivity 


o •= 10" 2 




Input the Height of 
the Antenna Feed 


Hq •= 2.5 






(meters) 
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Radiation Patterns valid for any Vertical Plane passing through the Antenna 



Space Wave Radiation Pattern 




xref^ , xlj 



Maximum Space Wave Electric Field 

Intensity (Volts per meter) max(MagEl) = 0.01733 
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Surface Wave Radiation Pattern 




-1 - 0.5 0 0.5 1 

Hret^ . h2j 



Maximum Surface Wave Electric Field 

Intensity (Volts per meter) max( MagE2) = 2.69275* 10 
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Combined Space and Surface Wave Radiation Pattern 




. k3. 



Maximum Radiated Electric Field 

Intensity (Volts per meter) max(MagE3) = 0.01732 
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X 

2 



power = 30- 



cos(/3'h'COs( f )) - cos(/3-h) 



’ - H 0 -co*(0) -J P (» + 



sin(£) 



-0 



R - Hq-cos( O 



R t H 0 



sqE3j := 



(l E 3 >l ) 2 

2 - ( 120 *) 



Directivity 



4'X'R^' max( sqE3) 
power 



Radres ;= 2- power 



Total Power Radiated (Watts) 



Radiation Resistance (Ohms) 



Directivity (maximum Power Gain 
assuming 100 % Antenna efficiency) 



Effective Isotropic Radiated 
Power (EIRP) (Watts) 



Maximum Effective Area 
(Along Radial of Directivity) 
(square meters) 



Maximum Effective Length 
(Along Radial of Directivity) 2* 
(meters) 



power = 10.79029 



Radres = 21.58059 



Directivity = 4.16784 



Directivity power = 44.97225 



A*') ^ Directivity 

= 33.16664 

4'X 



Radres- ( Ac') Directivity . 

^ = 2.75579 

480- tl 2 
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Numerical Distance for 
Vertical Polarization 



| Pe Q | = 74.95924 



Elevation Angle of Maximum 
Power Gain (Degrees) 



Angle3 



19.39778 \ 
19.39778/ 
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VERTICAL DIPOLE 



This application calculates far field radiation patterns and parameters associated 
with vertical thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted vertically along the z-axis in a rectangular coordinate system with the feed 
at the antenna’s center at a set height above the surface. Required inputs are the 
antenna length, feed height above the surface, transmitted frequency, distance from 
the antenna, the conductivity and dielectric constant of the surface below the 
antenna. The planar earth model is assumed in predicting radiation patterns. 
Predicted operating frequencies assume that the antenna is a quarter- wavelength, 
half- wavelength, three-quarter- wavelength, or full- wavelength dipole. A sinusoidal 
current input with a maximum of unity is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna. Plotted radiatioh patterns are valid for any vertical plane passing through 
the antenna axis. The electric field magnitudes to which the Radiation Patterns are 
normalized are displayed below their respective plots. Polarization is vertical for all 
vertical dipoles. 

Input the Dipole 1-5 h - - 

length in meters 2 

Wavelengths and Frequencies 



A t = 41 


X 2 21 


X 3 := 

3 3 


A 4 := 1 




f 2 = — 
X 2 


f 3 = — 
X 3 


li 

tM 


Possible Operating Frequencies (Hertz) 






^ Dipole: 
4 


f t = 1.5- 10 7 


A , 

- Dipole: 


f 2 = 3- 10 7 


3A IV 1 
— — Dipole: 

4 


f 3 = 4.5- 10 7 


A Dipole: 


f 4 = 6-10 7 


Input the operating 
Frequency (Hertz) 


f 5 30- 10 6 


Input the Distance 
from the Antenna 
(meters) 


R := 3000 


Input the ground 
Dielectric Constant 


CT := 72 


Input the ground 
Conductivity 


o=4 




Input the Height of 
the Antenna Feed 


Hq := 2.5 





(meters) 



160 



Radiation Patterns valid for any Vertical Plane passing through the Antenna 



Space Wave Radiation Pattern 




, « 1 | 



Maximum Space Wave Electric Field 

Intensity (Volts per meter) max(MagEl) = 0.03394 
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Surface Wave Radiation Pattern 




-1 -0.5 0 0.5 1 



Hret^ . h2j 

Maximum Surface Wave Electric Field 

Intensity (Volts per meter) max(MagE2) = 0.03331 
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Combined Space and Surface Wave Radiation Pattern 




, h3j 



Maximum Radiated Electric Field 

Intensity (Volts per meter) max(MagE3) = 0.03398 
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power ■= 30- 



2 



-0 



cos(/9-h-cos( f ) ) - cos(/3-h) 
sin( ?) 



-J p (R - H q oor(a) 

e v / 

R - H 0 -cos( ?) 




sqE3j ;= 



(l E 3 .l ) 2 

2 ( 120 *) 



Directivity 



4-x-R 2 -max( sqE3) 
power 



Radres := 2- power 



Total Power Radiated (W atts) 

Radiation Resistance (Ohms) 

Directivity (maximum Power Gain 
assuming 100% Antenna efficiency) 



power = 37.64371 
Radres = 75.28741 

Directivity = 4.59967 



Effective Isotropic Radiated 
Power (E1RP) (Watts) 



Maximum Effective Area 
(Along Radial of Directivity) 
(square meters) 



Maximum Effective Length 
(Along Radial of Directivity) 2- 
(meters) 



Directivity power = 173.14876 



A 

• Directivity 

^ = 36.60304 

4tc 



Radres- ( Ae') 2 Directivity 

^ = 5.40735 

480- x 2 
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Numerical Distance for 
Vertical Polarization 



| Pe 0 1 = 0.39252 



Elevation Angle of Maximum 
Power Gain (Degrees) 



Ang!e3 



8.55784 \ 
8.55784/ 
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APPENDIX B: 

VERTICAL MONOPOLE ARRAY COMPUTER OUTPUT 

This appendix contains computer hardcopies from the Mathcad 
vertical monopole application which show the input values and 
predicted radiation characteristics for two sample calculations. 
The configuration is a quarter-wave monopole above soil (e r =10 and 
<r=10‘ 2 ) for the first example and the same configuration above 
seawater (e r =72 and <7=4 ) for the second. Reference 6 [p. 89] 
provides the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The elevation of maximum directive gain is slightly higher 
for the Mathcad output than for those given in reference 6, but 
the overall radiation patterns are very similar. The 3.24 (5.1 
dB) value of directivity for Mathcad is quite a bit higher than 
the maximum gain values of about zero dB in reference 6. 
However, since a half-wave dipole has a free- space directivity of 
1.64 (2.1 dB) , the actual directivity of a quarter-wave monopole 
above a ground plane should be closer to the Mathcad prediction 
because of the effect of the reflected wave (constructive and 
destructive interference) on the space wave and the lower total 
average radiated power resulting from ground plane losses. The 
seawater example yields results consistent with expectations. 
With respect to the soil example, the seawater example's 
directivity is slightly higher due to a stronger reflected wave, 
and the surface wave is stronger at grazing angles (0«9O°) due to 
higher conductivity. 
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VERTICAL MONOPOLE 



This application calculates far field radiation patterns and parameters associated 
with vertical thin- wire monopole antennas (diameter « wavelength). The antenna 
is mounted vertically along the z-axis in a rectangular coordinate system with the 
feed at the origin. Required inputs are the antenna length, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. Predicted operating frequencies assume that the antenna is an eighth- 
wavelength, quarter-wavelength, three-eighths-wavelength, or half-wavelength 
monopole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electk field intensity 
transmitted by the antenna. The electric field magnitudes to which the patterns are 
normalized are displayed below their respective plots. Radiation patterns are valid 
for any vertical plane passing through the z-axis, because the radiation pattern is 
symmetrical with respect to phi. Polarization is vertical for all vertical monopoles. 



Input the Monopole 
length in meters 

Wavelengths and Frequencies 



1 2- h 



X, = 


41 


X 2 


= 21 


X, , ii 
3 3 


h 


= 1 


f, = 


c 

h 


f 2 


II 


f 3 •■= — 

X 3 


( 4 


c 



Possible Operating Frequencies (Hertz) 



^ Monopole: fj = 5*10* 



^ Monopole: f 3 = 1.5- 10 7 



^ Monopole: fg = 1 • 10^ 



^ Monopole: f| = 2* 10^ 



Input the operating 
Frequency (Hertz) 



f 5 := 10- 10 6 



Input the Distance from ^ 
the Antenna (meters) 



Input the ground 
Dielectric Constant 



Input the ground 
Conductivity: 



o := 



10 



-2 
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Radiation Patterns valid for any Vertical Plane passing through the Antenna 



Space Wave Radiation Pattern 




, Hi. 



Maximum Space Wave Electric Field 

Intensity (Volts per meter) max ( MagEl ) = 0.01 1 1 5 
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Surface Wave Pattern 




Kref^ . h2j 



Maximum Surface Wave Field Intensity 
(Volts per meter) 



max(MagE2) = 8.23677-10 
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Combined Space and Surface Wave Pattern 




-1 -0.5 0 0.5 1 



wefy , h3j 

Maximum Radiated Electric Field 

Intensity (Vohs per meter) max( MagE3) = 0.01 1 14 
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power ;= 



7.5 

7 * 



X 

2 



JO 



cost 


0hcos(f)) - cos(/3 h) ... 




sin(j8-h-cos(?)) ... 




i 


+ (-(cos(?)sin(0h))) 





sin(?) 



n -cos(f) - <> 



n 2 - sin(f) 2 



n -cos(f) 



*£» T 



sin(0' 



cos(/3*h*cos(f)) - cos(/3*h) 



V 



sm(j8-h cos(f)) •••! | 

+ (-(cos(Qsin(/3h)) 

stn(n 



sqE3j ;= 



(l E 3 ii ) 2 

2 ( 120 *) 



Directivity = 



4*x-R 2 -max( sqE3) 
power 



Radres ;= 2* power 

Total Power Radiated (Watts) power = 5.74454 

Radiation Resistance (Ohms) Radres = 11.48908 

Directivity (maximum Power Gain 

assuming 100% Antenna Efficiency) Directivity = 3.23898 



Effective Isotropic Radiated 

Power (EIRP) (Watts) Directivity • power = 18.60648 



Maximum Effective Area 
(Along Radial of Directivity) 
(square meters) 



(Xc') 2 - Directivity 

= 231.97519 

4-x 
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Maximum Effective Length 
(Along Radial of Directivity) 
(meters) 



Radres- (Xg) Directivity 



480*7r‘ 



5.31775 



Numerical Distance for 
Vertical Polarization 



|Pe Q | = 14.91116 



Elevation Angle of Maximum 
Power Gain (Degrees) 



Angle3 = 



25.67353 

25.67353 
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VERTICAL MONOPOLE 



This application calculates far field radiation patterns and parameters associated 
with vertical thin- wire monopole antennas (diameter « wavelength). The antenna 
is mounted vertically along the z-axis in a rectangular coordinate system with the 
feed at the origin. Required inputs are the antenna length, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. Predicted operating frequencies assume that the antenna is an eighth- 
wavelength, quarter-wavelength, three-eighths-wavelength, or half-wavelength 
monopole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electic field intensity 
transmitted by the antenna. The electric field magnitudes to which the patterns are 
normalized are displayed below their respective plots. Radiation patterns are valid 
for any vertical plane passing through the z-axis, because the radiation pattern is 
symmetrical with respect to phi Polarization is vertical for all vertical monopoles. 



Input the Monopole 
length in meters 



Wavelengths and Frequencies 



X t :• 41 


X 2 21 


x 3 ,±! 

3 3 


X 4 := 1 




f 2 := — 

X 2 


f 3 := — 

X 3 


'• ■ k 


Possible Operating Frequencies (Hertz) 






^ Monopole: 


f t = 5-10 6 


^ Monopole: 
4 


f 2 = 1- 10 7 


Monopole: 


f 3 = 1.5* 10 7 


^ Monopole: 


f 4 = 2I0 7 



Input the operating ^ jq jqG Input the Distance from ^ ^qqo 
F requency (Hertz) ® the Antenna (meters) 



Input the ground a - 12 Input the ground 
Dielectric Constant Conductivity: 
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Radiation Patterns valid for any Vertical Plane passing through the Antenna 



Space Wave Radiation Pattern 




, hIj 



Maximum Space Wave Electric Field 

Intensity (Volts per meter) max(MagEl) =0.01855 
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Surface Wave Pattern 




Href^ , h2j 



Maximum Surface Wave Field Intensity 

(Volts per meter) max(MagE2) = 0.01976 
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Combined Space and Surface Wave Pattern 




-1 -0.5 0 0.5 1 

Href^ , «3j 



Maximum Radiated Electric Field 

Intensity (Vohs per meter) max( MagE3) = 0.01976 
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** 

2 



power ;= 



7.5 

R 2 



.0 



cos(0hcos(f)) - cos(jfMi) 



+ )' 



sin(0*h*cos( f)) 

+ (-(cos(n-s«n()8h))) 



sin(H 



n -cos(f) 
n 2 *cos( f) 



JJ2 

+ Ja 2 - 



sm(0' 

sin( 0 ‘ 



cos(/3-h 

K 



cos(f)) - cos(/3-h) ... 
sin (/S*h*cos(f )) •••! I 

+ (-(cos(f)sin(fln)) 
sin(?) 



sqE3j ••= 



(l E3 il) 2 

2 -( 120 -*) 



Directivity = 



4-xR 2 max(sqE3) 

power 



Radres ;= 2- power 



Total Power Radiated (Watts) 



Radiation Resistance (Ohms) 



Directivity (maximum Power Gain 
assuming 100% Antenna Efficiency) 



Effective Isotropic Radiated 
Power (EIRP) (Watts) 



Maximum Effective Area 
(Along Radial of Directivity) 
(square meters) 



power - 16.19454 
Radres = 32.38909 

Directivity = 3.61503 
Directivity *power = 58.54372 

fXe) 2 ' Directivity 

v < ■-- = 258.9073 

4** 
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Maximum Effective Length 
(Along Radial of Directivity) 
(meters) 



2 - 



i ; 

Radres* c - Directivity 

480-x 2 



Numerical Distance for 
Vertical Polarization 



|Pe Q | = 0.04363 



Elevation Angle of Maximum 
Power Gain (Degrees) 



Angle 3 = 




9.4327 
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APPENDIX C: 

HORIZONTAL DIPOLE ARRAY COMPUTER OUTPUT 

This appendix contains computer hardcopies from the Mathcad 
horizontal dipole application which show the input values and 
predicted radiation characteristics for four sample calculations. 
The configuration is a half-wave dipole at one- quarter, one -half, 
and three-quarters wavelengths above soil (e r =lO and cr=l0' 2 ) for the 
first three examples, respectively, and a half-wave dipole at one- 
quarter wavelength above seawater (e r =72 and c= 4) for the fourth. 
Reference 6 [pp. 86-88] provides the radiation patterns and gain 
predictions from several sources for the configurations in the 
first three examples. 

The radiation patterns and maximum directive gain computed by 
the first three Mathcad examples are almost identical to the those 
given in reference 6. The directivity predictions are slightly 
higher for Mathcad, but the overall similarity between the 
predicted radiation characteristics is noteworthy. The seawater 
example also yields results consistent with expectations. With 
respect to the soil example, the seawater example's directivity is 
slightly higher due to a stronger reflected wave, and the surface 
wave is slightly stronger at grazing angles (0<= 90°) due to higher 
conductivity. For horizontal polarization, the higher conductivity 
surface below the antenna does not result in a greatly enhanced 
surface wave and increased directivity as it does for vertical 
polarization. 
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HORIZONTAL DIPOLE 



This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 

Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter-wavelength, half-wavelength, three-quarter-wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 

Input the dipole I - 15 h = - 

length in meters 2 



Input the index of the elevation angle for 
which to calculate the horizontal radiation 
pattern (from the angle index table) 



Wavelengths and Frequencies 



:■ 41 


*2 * 21 


x - 41 

X 3 - t 




f ' ' *i 


II 

CJ 

hmrn 


f 3 := — 

X 3 


f 4 • 


Possible Operating Frequencies (Hertz) 






^ Dipole: 
4 


f t = 5- 10 S 


^ Dipole: 


<2 = 


3X TV . 
Dipole: 

4 


f 3 = 1.5- 10 7 


X Dipole: 


'4 = 



180 



Input the operating 
Frequency (Hertz) 



f 5 := 10- 10 6 




0 * 



2- x 



X 



5 



Input the Height of 
Antenna Feed (meters) 



H 0 := 7.6 



Input the ground 
Dielectric Constant 



ct = 30 



Input the Distance 
from Antenna (meters) 



R := 3000 



Input the ground 
Conductivity: 



O := 3- 10" 2 



Index of Refraction 



Rdj := R - Hq-cos^0|) 
RTj := R + Hq'COS^0j^ 



n = 



CT - 



N 



i- 



18000o 

('s 10 ' 6 ) 



Complex Numerical Distance 
for Vertical Polarization 




Complex Numerical Distance 
for Horizontal Polarization 



Pmj := 



Vertical Reflection Coefficient 



Horizontal Reflection Coefficient 
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Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 




Max E-Field Intensity (Volts per meter) max( MagElP) = 0.03632 
Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 0.5 1 



Max E-Field Intensity (Volts per meter) max(MagE2P) = 1.03856- 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




Hret^ . h3Pj 

Max E-Fieldlntensity (Vohs per meter) max(MagE3P) = 0.03632 
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Radiation Patterns in Phi=0 Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

xref^ , k!0| 



Max E-Field Intensity (Volts per meter) max(MagElO) = 0.03594 
Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

Krel^ , x20j 



Max E-Field Intensity (Volts per meter) max(MagE20) = 3.84453- 10 * 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




xref^ , h30 



Max E-Field Intensity (Volts per meter) max(MagE30) = 0.03594 
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Radiation Patterns in Horizontal Plane 
Space Wave Radiation Pattern 




Phi=0 



-1 - 0.5 0 0.5 1 

Kfefj, . hIHj 

Max E-Field Intensity (Volts per meter) max(MagElH) = 0.02634 
Surface Wave Radiation Pattern (Horizontal Plane) 




Phi=0 



Kret^ , h2H| 



Max E-Field Intensity (Volts per meter) max(MagE2H) = 2.62147-10 



r 6 
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Combined Space and Surface Wave Radiation Pattern (Horizontal Plane) 




Phi=0 



Max E-Field Intensity (Volts per meter) max(MagE3H) — 0.02634 
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Total Power Radiated (W alts) 



power - 37.80286 



Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 



Radres = 75.60572 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 

Phi = pi/2 Wane Phi = 0 Plane 

DirectivityP = 5.23417 DirectivityO = 5.12535 

Effective Isotropic Radiated Power (EIRP) (W atts) 

Phi = pi/2 Plane Phi - 0 Plane 

DirectivityP- power = 197.86674 DirectivityO- power = 193.75275 

Maximum Effective Area (Along Radial of Directivity) 

(square meters) 

Phi = pi/2 Plane Phi = 0 Plane 

( Ac) DirectivityP ( Ac) DirectivityO 

= 374.87006 = 367.07586 

4-x 4-x 
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Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres- Directivity? 



= 17.34132 2 



480-x‘ 



Radres* (A«j) 2 * DirectivityO 



480* ic 4 



Numerical Distances 



Vertical Polarization 



Horizontal Polarization 



PeJ = 5.04617 



|Pm Q | = 1.92562- 10 4 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 



65.8954 \ 
65.8954/ 



AngleO = 



89.85737 

89.85737 



17.1601 
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HORIZONTAL DIPOLE 



This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 

Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter- wavelength, half- wavelength, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 

Input the dipole 1 = 75 h = - 

length in meters 2 

Input the index of the elevation angle for 

which to calculate the horizontal radiation d := 535 

pattern (from the angle index table) 

Wavelengths and Frequencies 



X. := 41 










Possible Operating Frequencies (Hertz) 



^ Dipole: f« = 1* 10 7 

4 ' 



^ Dipole: f 2 = 2-10 7 



3-X 7 

~ Dipole: f 3 = 3- 1(V 



X Dipole: f^ = 4*10 



7 
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Input the operating 
Frequency (Hertz) 



f 5 := 20- 10* 




0 = 



2n 



X 



5 



Input the Height of 
Antenna Feed (meters) 




Input the ground 
Dielectric Constant 



a 30 



Input the Distance 
from Antenna (meters) 



R := 3000 



Input the ground 
Conductivity: 



o := 3- 10" 2 



Index of Refraction 



Rdj := R - Hfl’cos(0j) 
RTj = R + Hq'COs(0j^ 



n 




j- 



18000 - o 

M-‘) 



Complex Numerical Distance 
for Vertical Polarization 




Complex Numerical Distance 
for Horizontal Polarization 



Pm, 



-j/3-Rr, 

2-sin^0,) 2 



cos 



W - J " 2 - ( s “( 9 i) 2 ) 



Vertical Reflection Coefficient 




Horizontal Reflection Coefficient 
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Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 



KlPj 

Max E-Field Intensity (Volts per meter) max(MagElP) = 0.03737 
Surface Wave Radiation Pattern (Phi=pi/2) 




. h2Pj 

Max E-Field Intensity (Volts per meter) max(MagE2P) = 8.03309’ 10 ^ 
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Radiation Patterns in Phi=0 Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 




-1 - 0.5 0 0.5 1 

xret^, xtOj 

Max E-Field Intensity (Volts per meter) max(MagElO) = 0.01561 
Surface Wave Radiation Pattern (Phi=0) 




xret^ , x20j 

Max E-Field Intensity (Volts per meter) max(MagE20) = 1.54902* 10 * 
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Radiation Patterns in Horizontal Plane 
Space Wave Radiation Pattern 




-1 -0.5 0 0.5 1 



wef^ . kIH. 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagElH) = 0.037B4 
Surface Wave Radiation Pattern (Horizontal Plane) 




-1 -0.5 0 0.5 1 

Kret^ , k2Hj 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagE2H) = 1.62046- 10 ^ 
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Total Power Radiated (Watts) 



power = 27.68975 



Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 



Radres = 55.3795 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 

Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP = 7.56408 DirectivityO = 1.32021 

Effective Isotropic Radiated Power (EIRP) (W atts) 

Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP- power = 209.44744 DirectivityO- power = 36.55637 

Maximum Effective Area (Along Radial of Directivity) 

(square meters) 

Phi = pi/2 Plane Phi = 0 Plane 

( At) DirectivityP ( At) DirectivityO 

V } A = 135.43431 ±-^—2 = 23.63833 

4-x 4 -x 
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Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres* ( Xc) • DirectivityP 
^ = 8.92079 2* 



480 *' 



Radres* 2 - DirectivityO 



480*x‘ 



Numerical Distances 



Vertical Polarization 



Horizontal Polarization 



PeJ = 15.283 



Pm 0 | = 2.48%* 10 4 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 



28.24089 \ 
28.24089 / 



AngleO = 



/ 41.64818\ 
\ 41.64818/ 



3.7269 
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HORIZONTAL DIPOLE 



This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 

Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter- wavelength, half- wave length, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. AU 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 

Input the dipole 1^5 h = - 

length in meters 2 

Input the index of the elevation angle for 

which to calculate the horizontal radiation d = 535 

pattern (from the angle index table) 



Wavelengths and Frequencies 



ii 


x 2 = 2 1 


x - 4,1 
X 3 - t 


\ t . i 


" - k 


A 2 


f 3 = ^ 
A 3 




Possible Operating Frequencies (Hertz) 






- Dipole: 
4 


iy = 1.5* 10 7 


^ Dipole: 


h = 3- 10 7 


3X TV . 
Dipole: 


f 3 = 4.5- 10 7 


X Dipole: 


f 4 = 610 7 
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Inpul (he operating 
Frequency (Hertz) 



f, 30 10* X 5 := f 

*3 



& -- 



2 * 



\ 



5 



Input the Height of 
Antenna Feed (meters) 



H 0 := 7.6 



Input the ground 
Dielectric Constant 



ex- 30 



Input the Distance 
from Antenna (meters) 



R := 3000 



Input the ground 
Conductivity: 



o = 3 10“ 2 



Index of Refraction 



Rdj = R - Hq-cos^Aj) 
Rfj := R + Hq COS^j^ 



n := 



CT - 



j- 



18000-0 

( f 5- *) 



Complex Numerical Distance 
for Vertical Polarization 






n 



/ 



Complex Numerical Distance 
for Horizontal Polarization 



PflQ| 



Vertical Reflection Coefficient 



Horizontal Reflection Coefficient 
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Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 

kIP. 

Max E-Field Intensity (Volts per meter) max(MagElP) = 0.03799 
Surface Wave Radiation Pattern (Phi=pi/2) 




, «2Pj 

Max E-Field Intensity (Volts per meter) max( MagE2P) = 6.21706* 10 7 
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Radiation Patterns in Phi=0 Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

hIOj 

Max E-Field Intensity (Volts per meter) max( MagElO ) = 0.03422 
Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

Mref^. h20| 



Max E-Field Intensity (Volts per meter) max(MagE20) = 9.38091*10 




Radiation Patterns in Horizontal Plane 
Space Wave Radiation Pattern 




-1 -0.5 0 0.5 1 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagElH) = 0.02972 
Surface Wave Radiation Pattern (Horizontal Plane) 




xrel^ , k2H| 

Max E-Field Intensity (Volts per meter) max(MagE2H) = 1.14409* 10 ^ 
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Total Power Radiated (Watts) 


power = 31.69054 


Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 


Radres = 63.38107 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 



Phi = pi/2 Plane 


Phi = 0 Plane 


DirectivityP = 6.8304 


DirectivilyO = 5.54301 



Effective Isotropic Radiated Power (EIRP) (W atts) 

Phi - pi/2 Plane Phi = 0 Plane 

Directivity P- power = 216.45893 Directi vityO- power = 175.66087 

Maximum Effective Area (Along Radial of Directivity) 

(square meters) 



Phi = pi/2 Plane 


Phi = 0 Plane 


( Xc) DirectivityP 

V J - 54.35457 

4-x 


( Xe) 2 * DirectivilyO 

v ' - 44.10985 

4-* 
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Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres* (X^) 2 * DirectivityP 






6.04592 2- 



480*x‘ 



Radres* ^* DirectivityO 



480*x‘ 



Numerical Distances 



Vertical Polarization 



Horizontal Polarization 



Pe Q | = 26.28165 



Pm 0 1 = 3.21687* 10’ 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 



18.542 \ 
18.542/ 



AngleO = 



89.85737 \ 
89.85737/ 



5.44643 
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HORIZONTAL DIPOLE 



This application calculates far field radiation patterns and parameters associated with 
horizontal thin-wire dipole antennas (diameter « wavelength). The antenna is 
mounted above and parallel to the x-axis in a rectangular coordinate system. The 
feed is at the center of the antenna at a set height directly above the origin. 

Required inputs are the antenna length, feed height above the surface, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. Predicted operating frequencies assume that the antenna is a 
quarter- wavelength, half- wavelength, three-quarter- wavelength, or full wave- 
length dipole. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. The radiation patterns are plotted for the phi=0 and phi=pi/2 
vertical planes parallel and perpendicular to the x-axis respectively. A horizontal 
radiation pattern is plotted at an elevation selected by the index from the elevation 
angle index table. Polarization is a combination of vertical and horizontal 
depending upon spatial orientation relative to the dipole. 

Input the dipole 1-3 h = - 

length in meters 2 

Input the index of the elevation angle for 

which to calculate the horizontal radiation d = 535 

pattern (from the angle index table) 

Wavelengths and Frequencies 



A, = 41 










Possible Operating Frequencies (Hertz) 



Dipole: f t = 2.5* 10 7 




Dipole: f 3 = 7.5- 10 7 



A Dipole: f 4 = 1*10 
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Input the operating j = 50> jqS 
Frequency (Hertz) 5 



X 5 r 



J8 = 



2-x 



Input the Height of H = 1 5 * n P ul g rount ^ 

Antenna Feed (meters) ® Dielectric Constant 



CT := 72 



Input the Distance ^ _ 3 qqq Input the ground o - A 

from Antenna (meters) Conductivity: 



Rdj = R - Hq'Cos(0j) 
Rr, = R + Hq- COS^0jj 



Index of Refraction 



n := 



cr - )• 



18000 - o 

( V 10- 6 ) 



Complex Numerical Distance 
for Vertical Polarization 



Pe, 



- j *0*Rrj 
2-sin(0j ) 1 



cos 



« * 1 



n 2 - sin(0|)' 



omplex Numerical Distance p . ' cos ( e ) * J„ 2 - (sio(0,) 2 ) 

)r Horizontal Polarization ' -> ' ^ '/ y \ \ V / 

Z * S 10 / C 7 |j 



Vertical Reflection Coefficient 



Horizontal Reflection Coefficient 



,. T .. (" 2 - cos ( 9 i)) - ( 


J - sln M 2 ) 


(n 2 . cos (<!,)) 4- ( 

, L .. cos 0’i) - (i 


Jn 2 - s!n(«,) 2 ) 


I" 2 - sin («l) 2 ) 


* "1 

cos(0|) 


+ G 


fn 2 - sin(0,) 2 ) 
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Radiation Patterns in Phi=Pi/2 Plane (Perpendicular to Dipole) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 

wet^ . hIPj 



Max E-Field Intensity (Volts per meter) max(MagElP) = 0.03926 
Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 0.5 1 

, «2Pj 



Max E-Field Intensity (Volts per meter) max(MagE2P) = 8.83121-10^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




-t -0.5 0 0.5 1 

. k3Pj 



Max E-Field Intensity (Vohs per meter) raax(MagE3P) = 0.03926 
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Radiation Patterns in Phi=0 Plane (Parallel to Dipole) 
Space Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

xref^. utOj 

Max E-Field Intensity (Volts per meter) max(MagElO) = 0.03925 
Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

«ret^ , h20j 



Max E-Field Intensity (Volts per meter) max(MagE20) = 3.2224*10 



Combined Space and Surface Wave Radiation Pattern (Phi=0) 




, h30| 

Max E-Field Intensity (Volts per meter) max ( MagE30 ) = 0.03925 
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Radiation Patterns in Horizontal Plane 
Space Wave Radiation Pattern 




-1 -0.5 0 0.5 1 




Max E-Field Intensity (Volts per meter) max( MagElH ) = 0.02647 
Surface Wave Radiation Pattern (Horizontal Plane) 




-1 -0.5 0 0.5 1 

xref^ , x2H| 



Phi=0 



Max E-Field Intensity (Volts per meter) 



max( MagE2H) = 4.07737- 10~ 8 
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Combined Space and Surface Wave Radiation Pattern (Horizontal Plane) 




xret^ , x3Hj 

Max E-Field Intensity (Volts per meter) max(MagE3H) = 0.02647 



Phi=0 
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Total Power Radiated (Watts) 



power = 41.98012 



Radiation Resistance (versus 
Maximum Antenna Current) 
(Ohms) 



Radres = 83.96025 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 

Phi = pi/2 Plane Phi = 0 Plane 

DireclivityP = 5.50819 Directi vityO = 5.50321 

Effective Isotropic Radiated Power (EIRP) (Watts) 

Phi = pi/2 Plane Phi = 0 Plane 

DireclivityP- power = 231.23457 DirectivityO* power = 231.02563 



Maximum Effective Area (Along Radial of Directivity) 
(square meters) 



Phi = pi/2 Plane 

( A s ) 2 * DireclivityP 

= 15.77981 

4*x 



Phi = 0 Plane 

( Ag) 2 * DirectivityO 
4 * 



15.76555 
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Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres* ( Xk) • Directivity? 

^ = 3.74932 2* 



480-x 4 



Radres* • DirectivityO 



480* x 4 



Numerical Distances 



Vertical Polarization Horizontal Polarization 

Pe 0 | = 1.08943 |Pm 0 | = 2.26469- 10 6 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 

Phi = pi/2 Plane Phi = 0 Plane 



AngleP = 



78.73217 \ 
78.73217/ 



AngleO = 



89.85737 \ 
89.85737 / 



3.74762 
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APPENDIX D: 

VERTICAL LOG-PERIODIC DIPOLE ARRAY COMPUTER OUTPUT 



This appendix contains computer hardcopies from the 
Mathcad vertical log-periodic dipole array application which 
show the input values and predicted radiation characteristics 
for two sample calculations. The configuration is a given by 
the inputs on the first two pages of each printout. The first 
antenna is mounted above soil (e r =10 and <7=10' 3 ) and the second 
above seawater (e r =72 and cr= 4). Reference 6 [p. 114] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, and the elevation of the maximum 
directive gain is also slightly higher, but the overall 
similarity between the predicted radiation characteristics is 
still quite good. The seawater example also yields results 
consistent with expectations. With respect to the soil 
example, the seawater example's directivity is higher due to 
a stronger reflected wave, and the surface wave is stronger at 
grazing angles (0«= 90°) due to higher conductivity. 
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VERTICAL LOG PERIODIC DIPOLE ARRAY 



This application calculates far field radiation patterns and parameters associated with 
vertical log periodic dipole arrays. The antenna is oriented such that the projection 
of its center axis lies on the positive y-axis of a rectangular coordinate system. The 
antenna axis can be oriented with respect to the vertical at any angle between zero 
and ninety degrees. The feed is at the center of the shortest element directly above 
the origin. The dipole elements are bisected by the antenna axis and are parallel to 
the z-axis. Required inputs are the number of dipole elements, shortest and second 
shortest element lengths, separation between the shortest and second shortest 
elements, radius of the shortest element, height of the feed above the surface, 
angle of the antenna axis with the vertical, characteristic admittance of the line 
feeding the antenna, termination impedance of the antenna, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=0 and phi=pi/2 vertical 
planes perpendicular and parallel to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is vertical for all vertical log-periodic antennas. 



Input the number 
of elements 



Input the length 
of the shortest 
element (meters) 

Input the distance 
from shortest to 
second shortest 
element (meters) 



N 12 



Input the shortest 
element's Radius 
(meters) 



radg ■•= .00918 



Input the length 

Iq := 3.673 of the second ^ = 4.373 

shortest element 
(meters) 

Input the height of 

dg = 1 .0766 (he shortest element Hg := 6.87 
above the surface 
(meters) 



Input the operating 
Frequency (Hertz) 



18 10 6 



Input the Distance 
from the Antenna 
(meters) 



Input the ground 
Conductivity 



j IQ- 3 Input the ground 
Dielectric Constant 



R := 3000 
cr - 10 



Input the Characteristic 1 

Admittance of the Line ADM = — 
Feeding the Antenna 
(Ohms) 



Input the Termination 

Impedance connected TIMP = 450 

to longest element 

(Ohms) 
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Input the angle between Input the elevation index 

the vertical (z axis) and ¥ = 78 for which to calculate the w : 285 

the antenna axis (degrees) horizontal radiation pattern 

(from angle index table) 

NOTE: The Vertical Log-Periodic Dipole mutual impedance calculations are not 
valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 



If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables 1 and d, using the define key (shift, colon), and entering the lengths and 
spacings by separating each successive entry by a comma. 



Calculated length 
of elements from 
shortest to longest 



Calculated distance 
between successive 
elements-shortest 
to longest 



Calculated radii 
of elements from 
shorlest to longest 



1 



d 



rad 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X -Axis) 

Space Wave Radiation Pattern (Phi=Pi/2) 




-t -0.5 0 0.5 1 



. MtPj 

Max E-Field Intensity (Volts per meter) max(MagElP) = 0.05599 
Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 05 1 



wel^ , k2P { 

Max E-Field Intensity (Volts per meter) max(MagE2P) = 0.00124 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




-t -0.5 0 0.5 1 

. h3Pj 



Max E-Field Intensity (Volts per meter) max(MagE3P) = 0.05583 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 




kIO. 

Max E-Field Intensity (Volts per meter) max( MagElO) = 0.03221 
Surface Wave Radiation Pattern (Phi=0) 




wef^ . k20j 

Max E-Field Intensity (Volts per meter) max(MagE20) = 7.64667 1 10 * 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

Hfet^ , k 30 j 

Max E-Field Intensity (Volts per meter) max(MagE30) = 0.0321 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 



wef^ , hIHj 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagElH) = 0.05438 
Surface Wave Radiation Pattern (Horizontal) 




. u2Hj 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagE2H) - 1.09861*10 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 




, h3H. 

Max E-Field Intensity (Volts per meter) max(MagE3H) = 0.05427 



Phi=0 
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wer0 := 



30 



*-R‘ 



x 



2 







j-p-Y -s* c(a2 + a3)-( cos(i; )-sfn( a2 + a.3) + sfn( 

1 e Y 
T 



power ;= powers 



Radres - 2 power 



sqE3Pj = 



(| MagE3P[|) 2 
2- ( 120- x) 



DirectivityP > 



4* x- R^- max( sqE3P) 
power 



sqE30| = 



(|MagE30,|) 2 
2( 120 -x) 



DirectivilyO 



4- x- R^- max( sqE30 ) 
power 



Total Power Radiated (Watts) 



power = 34.29572 



Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 



Radres = 68.59144 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 

Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP = 13.63346 DirectivityO = 4.50794 

Effective Isotropic Radiated Power (EIRP) (W atts) 

Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP- power = 467.56917 DirectivityO- power = 154.6031 
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Maximum Effective Area (Along Radial of Directivity) 
(square meters) 



Phi = pi/2 Plane 



Phi = 0 Plane 



(X 5 ) 2 -DirectivityP 



4-x 



= 301.36553 



2 - DireclivilyO 



4-x 



= 99.64739 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 



Phi = pi/2 Plane 



Phi = 0 Plane 



Radres- ( Ac) • DirectivityP 
^ = 14.8097 2- 



480- ic* 



Radres- (Ag) 2 - DirectivityO 



480- r' 



Numerical Distance for Vertical Polarization 
Pe 0 1 = 50.7 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 14.28571 



AngleO = 



1 3.71 429 \ 
13.71429/ 



8.51594 
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VERTICAL LOG PERIODIC DIPOLE ARRAY 



This application calculates far field radiation patterns and parameters associated with 
vertical log periodic dipole arrays. The antenna is oriented such that the projection 
of its center axis lies on the positive y-axis of a rectangular coordinate system. The 
antenna axis can be oriented with respect to the vertical at any angle between zero 
and ninety degrees. The feed is at the center of the shortest element directly above 
(he origin. The dipole elements are bisected by the antenna axis and are parallel to 
the z-axis. Required inputs are the number of dipole elements, shortest and second 
shortest element lengths, separation between the shortest and second shortest 
elements, radius of the shortest element, height of the feed above the surface, 
angle of the antenna axis with the vertical, characteristic admittance of the line 
feeding the antenna, termination impedance of the antenna, transmitted frequency, 
distance from the antenna, the conductivity and dielectric constant of the surface 
below the antenna. The planar earth model is assumed in predicting radiation 
patterns. A sinusoidal current input with a maximum of unity is assumed. All 
radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=0 and phi=pi/2 vertical 
planes perpendicular and parallel to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is vertical for all vertical log-periodic antennas. 



Input the number 


N : = 


12 


Input the shortest 


radg := .00918 


of elements 


element's Radius 
(meters) 








Input the length 


*0 = 


3.673 


Input the length 


l t := 4.373 


of the shortest 


of the second 


element (meters) 






shortest element 
(meters) 




Input the distance 




1.0766 


Input the height of 


H 0 := 6.87 


from shortest to 


d Q •= 


the shortest element 


second shortest 






above the surface 




element (meters) 






(meters) 




Input the operating 


f 5 = 


18- 10 6, 


Input the Distance 


R := 3000 


Frequency (Hertz) 


from the Antenna 
(meters) 








Input the ground 


0 •= 


4 


Input the ground 


CT := 72 


Conductivity 






Dielectric Constant 


Input the Characteristic 




1 


Input the Termination 


TIMP ■•= 450 


Admittance of the Line 


ADM 


450 


Impedance connected 


Feeding the Antenna 




to longest element 




(Ohms) 






(Ohms) 
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Input the angle between Input the elevation index 

the vertical (z axis) and ^ •- 75 for which to calculate the w = 285 

the antenna axis (degrees) horizontal radiation pattern 

(from angle index table) 

NOTE: The Vertical Log-Periodic Dipole mutual impedance calculations are not 
valid if (here is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 



If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables 1 and d, using the define key (shift, colon), and entering the lengths and 
spacings by separating each successive entry by a comma. 



Calculated length 
of elements from 
shortest to longest 



Calculated distance 
between successive 
elements-shortest 
to longest 



Calculated radii 
of elements from 
shortest to longest 



I 



d 



rad 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X -Axis) 

Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 

, hIPj 



Max E-Field Intensity (Volts per meter) max(MagElP) = 0.10756 
Surface Wave Radiation Pattern (Phi=pi/2) 




Kref^ . k2P. 

Max E-Field Intensity (Volts per meter) max(MagE2P) = 0.12392 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 0.5 1 

w . k3Pj 



Max E-Field Intensity (Volts per meter) max( MagE3P) = 0.12392 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 




-1 -0.5 0 0.5 1 

xr*^. xIOj 

Max E-Field Intensity (Volts per meter) max(MagElO) = 0.06687 
Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 



xret^ . x20j 

Max E-Field Intensity (Volts per meter) max ( MagE20 ) = 0.07827 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

nret^ , x30j 

Max E-Field Intensity (Volts per meter) max(MagE30) = 0.07827 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 



Phi=0 



, kIHj 

Max E-Field Intensity (Volts per meter) max(MagElH) = 0.06821 
Surface Wave Radiation Pattern (Horizontal) 




uref^ , «2H. 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagE2H) - 4.93203-10 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 




. h3H. 

Max E-Field Intensity (Volts per meter) max(MagE3H) = 0.0682 



Phi=0 
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ve T0 : = 



30 

ic-R‘ 





r 'w 


r nr 






2 


2 




Z 




1C 


0 




y 










2 







I -e 
Y 



j-p-Y -$ec(a2 + a3) ( cos(^ ) sm(a2 + a3) i 



power := powers 



Radres = 2 power 



sqE3P| := 



( | MagE3P j | ) 

2 -( 120 - it ) 



DireclivityP = 



4-x-R -max( sqE3P) 
power 



sqE30j 



( | MagE30j | ) ' 



2-( 120- w ) 
Total Power Radiated (Watts) 



Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 



DirectivityO = 



4-x R -max( sqE30) 



power 
power = 99.48703 

Radres = 198.97405 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 

Phi = pi/2 Plane Phi = 0 Plane 

Directivity P = 23.15454 DirectivityO = 9.2377 

Effective Isotropic Radiated Power (EIRP) (Walts) . 

Phi = pi/2 Plane Phi = 0 Plane 

DireclivityP- power = 2.30358*10^ DirectivityO- power = 919.030% 



sin(C) 
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Maximum Effective Area (Along Radial of Directivity) 
(square meters) 



Phi = pi/2 Plane 

/ X*?') DirectivhyP 

= 511, 

4- x 



82763 



Phi = 0 Plane 

( Xe) DireclivityO 

= 204.19793 

4-x 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres- ( Xe) • DirectivityP 
= 32.87188 2 



480- x‘ 



Radres- ( Xe\ • DirectivityO 
= 20.76291 



480* x* 



Numerical Distance for Vertical Polarization 



Pe 0 | = 0.14135 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 0 



AngleO = 



0 

0 , 
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APPENDIX E: 

HORIZONTAL LOG- PERIODIC DIPOLE ARRAY COMPUTER OUTPUT 

This appendix contains computer hardcopies from the 
Mathcad horizontal log-periodic dipole array application which 
show the input values and predicted radiation characteristics 
for two sample calculations. The configuration is a given by 
the inputs on the first two pages of each printout. The first 
antenna is mounted above soil (e r =4 and a=10' 3 ) and the second 
above seawater (e r =72 and a= 4). Reference 6 [p. 108] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, but the overall similarity among 
the predicted radiation patterns is remarkably good. The 
seawater example has results very much like the soil example. 
For horizontal polarization, the higher conductivity surface 
below the antenna does not result in a greatly enhanced 
surface wave and increased directivity as it does for vertical 
polarization. 
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HORIZONTAL LOG PERIODIC DIPOLE ARRAY 



This application calculates far field radiation patterns and parameters associated with 
horizontal log periodic dipole arrays. The antenna is oriented such that the 
projection of its center axis lies on the positive y-axis of a rectangular coordinate 
system. The antenna axis can be oriented with respect to the vertical at any angle 
between zero and ninety degrees. The feed is at the center of the shortest element 
directly above the origin. The dipole elements are bisected by the antenna axis and 
are parallel to the x-axis. Required inputs are the number of dipole elements, 
shortest and second shortest element lengths, separation between the shortest and 
second shortest elements, radius of the shortest element, height of the feed above 
the surface, angle of the antenna axis with the vertical, characteristic admittance of 
the line feeding the antenna, termination impedance of the antenna, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. A sinusoidal current input with a maximum of unity is assumed . 
All radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=pi/2 and phi =0 vertical 
planes parallel and perpendicular to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is a varying combination of horizontal and vertical 
depending on spatial orientation with respect to the antenna. 



Input the number ^ ^ 

of elements 

Input the length 

of the shortest = 7.219 

element (meters) 



Input the distance 

from shortest to ;= 1.659 

second shortest 
element (meters) 



Input the operating 
Frequency (Hertz) 



f 5 := 12 * 10 & 



Input the shortest 

element's Radius rat ^o " 001^4 
(meters) 



Input the length 

of the second ^ = 8.297 

shortest element 

(meters) 

Input the height of 

the shortest element Hq ■= 9.4 

above the surface 

(meters) 

Input the Distance 

from the Antenna R := 3000 

(meters) 



Input the ground _ j , jq- 3 

Conductivity 



Input the ground 
Dielectric Constant 



Input the Characteristic 
Admittance of the Line 1 

Feeding the Antenna " 300 

(Ohms) 



Input the Termination 
Impedance connected 
to longest element 
(Ohms) 



TIMP 



300 
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Input the angle between Elevation index for which 

the vertical (z axis) and ♦ = 90 to calculate horizontal pattern w r 285 

the antenna axis (degrees) (from angle index table) 



NOTE: The Horizontal Log-Periodic Dipole mutual impedance calculations are not 
valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 



If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables 1 and d, using the define key (shift, colon), and entering the successive 
lenglhs and spacings by separating each successive entry by a comma. 



Calculated length 
of elements from 
shortest to longest 



Calculated distance 

between successive Calculated radii 

elements-shortest elements from 

to longest shortest to longest 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 

, kIPj 



Max E-Field Intensity (Volts per meter) max(MagElP) = 0.08512 
Surface Wave Radiation Pattern (Phi=pi/2) 




xre^ . k2P ( 

Max E-Field Intensity (Volts per meter) max(MagE2P) = 4.54242-10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




xref^ . x3P| 

Max E-Field Intensity (Volts per meter) raax( MagE3P ) = 0.0851 1 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 

Space Wave Radiation Pattern (Phi=0) 




Kfet^ , KlO. 

Max E-Field Intensity (Volts per meter) max(MagElO) = 0.03495 
Surface Wave Radiation Pattern (Phi=0) 




-1 - 0.5 0 0.5 1 



Hret^ , k20 ( 

Max E-Field Intensity (Volts per meter) max(MagE20) = 1.45261-10 * 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

, h30j 



Max E-Field Intensity (Volts per meter) max(MagE30) = 0.03495 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 

Kfet^ , hIH. 

Max E-Field Intensity (Volts per meter) max(MagElH) = 0.0676 
Surface Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 

. k2H. 



Max E-Field Intensity (Volts per meter) max(MagE2H) = 3.44012* 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 

, k3Hj 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagE3H) = 0.06758 
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power# ••= 



30 



x-R‘ 





* 1.5707963 


* 1.5707963 




S 










y 




- 1.5707963 - 


0 







cos(()*cos(?)’Iy j P Y csc( 

c 



1 - (sin(f) cos(O ) 2 



powers ;= 



30 



x-R‘ 



x 

2 



-- -0 



x 

2 



E 



sio( $ )I 



j-p-Y •csc(4 / )•( cos(t;) co*( 



1 - (sin(n*cos( 0 )' 



power - power# +- power# 



Radres = 2 power 



sqE3P| = 



( | MagE3P, | ) ' 



2*( 120x) 



Direct! vityP := 



4-i-R -max( sqE3P) 



power 



= (|MagE3°i[) 2 Djreci . = 

' 2(120x) power 



Total Power Radiated (W atts) power = 83.98604 

Radiation Resistance (versus Ra( | res = 167.97209 

Maximum Input Current) 

(Ohms) 

Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 

Directi vityP = 12.93749 DireclivityO = 2.18175 
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Effective Isotropic Radiated Power (EIRP) (Watts) 



Phi = pi/2 Plane Phi = 0 Plane 

DireclivityP- power = 1.08657* 10 3 DirectivityO- power = 183.23657 

Maximum Effective Area (Along Radial of Directivity) 

(square meters) 

Phi = pi/2 Plane Phi = 0 Plane 



( Aj) DireclivityP 
4-x 



= 643.45795 



DirectivityO 

4-ff 



108.51135 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres- ( Ar) • DirectivityP 
^ = 33.86438 2- 



480-x‘ 



Radres- (A«-) DirectivityO 



480- k 2 



13.90659 



Numerical Distances 



Vertical Polarization 



Pe Q I = 69.28584 



Horizontal Polarization 
Pm 0 | = 1.26447* 10 3 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 

Phi = pi/2 Plane Phi = 0 Plane 



AngleP = 31.42857 



AngleO = 



89.71 429 \ 
89.71429/ 
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HORIZONTAL LOG PERIODIC DIPOLE ARRAY 



This application calculates far field radiation patterns and parameters associated with 
horizontal log periodic dipole arrays. The antenna is oriented such that the 
projection of its center axis lies on the positive y-axis of a rectangular coordinate 
system. The antenna axis can be oriented with respect to the vertical at any angle 
between zero and ninety degrees. The feed is at the center of the shortest element 
directly above the origin. The dipole elements are bisected by the antenna axis and 
are parallel to the x-axis. Required inputs are the number of dipole elements, 
shortest and second shortest element lengths, separation between the shortest and 
second shortest elements, radius of the shortest element, height of the feed above 
the surface, angle of the antenna axis with the vertical, characteristic admittance of 
the line feeding the antenna, termination impedance of the antenna, transmitted 
frequency, distance from the antenna, the conductivity and dielectric constant of the 
surface below the antenna. The planar earth model is assumed in predicting 
radiation patterns. A sinusoidal current input with a maximum of unity is assumed . 
All radiation patterns are normalized with respect to the maximum electric field 
intensity transmitted by the antenna in the plane of interest. The electric field 
magnitudes to which the patterns are normalized are displayed below their 
respective plots. Radiation patterns are plotted for the phi=pi/2 and phi=0 vertical 
planes parallel and perpendicular to the y-axis respectively. A horizontal radiation 
pattern is plotted at an elevation selected by the index from the elevation angle 
index table. Polarization is a varying combination of horizontal and vertical 
depending on spatial orientation with respect to the antenna. 



Input the number ^ ^ 

of elements 

Input the length 

of the shortest *o := 7.219 

element (meters) 



Input the distance 

from shortest to ;= 1.659 

second shortest 

element (meters) 



Input the operating 
Frequency (Hertz) 



f 5 := 12- 10 6 



Input the shortest 

element's Radius rat *g = -00*44 
(meters) 



Input the length 

of the second *i := 8-297 

shortest element 

(meters) 

Input the height of 

the shortest element h 0 •= 9.4 

above the surface 

(meters) 



Input the Distance 
from the Antenna 
(meters) 



R := 3000 



Input the ground 
Conductivity 



Input the ground 
Dielectric Constant 



Input the Characteristic 
Admittance of the Line 1 

Feeding the Antenna " to 

(Ohms) 



Input the Termination 
Impedance connected 
to longest element 
(Ohms) 



TIMP 



300 
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Input the angle between Elevation index for which 

the vertical (z axis) and ¥ - 90 to calculate horizontal pattern w = 285 

the antenna axis (degrees) (from angle index table) 



NOTE: The Horizontal Log-Periodic Dipole mutual impedance calculations are not 
valid if there is an antenna element whose length is an exact integer multiple of the 
wavelength. If this occurs, there will be a singularity error in the mutual impedance 
calculations. If this problem arises, it will be necessary to vary the operating frequency 
such that no element is exactly an integer multiple of the wavelength. 



If the log-periodic calculations of element length and spacing do not represent the 
desired antenna configuration, you can enter the values directly by selecting the 
variables 1 and d, using the define key (shift, colon), and entering the successive 
lengths and spacings by separating each successive entry by a comma. 



Calculated length 
of elements from 
shortest to longest 



Calculated distance 

between successive Calculated radii 

elements-shortest °f elements from 

to longest shortest to longest 



1 

jr 



7.219 



8.297 



9.53598 



10.95997 



12.5966 



14.47762 



16.63954 



19.12429 



21.98009 



25.26233 



29.03471 



33.37041 



d 

v 

1.659 

1.90674 

2.19146 

2.51871 

2.89483 

3.32711 

3.82394 

4.394% 

5.05125 

5.80554 

6.67247 



rad 

I 



0.00144 


0.00166 


0.0019 


0.00219 


0.00251 


0.00289 


0.00332 


0.00381 


0.00438 


0.00504 


0.00579 


0.00666 



247 




Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 - 0.5 0 0.5 1 

Hfet^ , hIPj 



Max E-Field Intensity (Volts per meter) max(MagElP) = 0.13724 
Surface Wave Radiation Pattern (Phi=pi/2) 




Max E-Field Intensity (Volts per meter) max(MagE2P) = 3.19958*10 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




, k3Pj 

Max E-Field Intensity (Volts per meter) max ( MagE3P ) = 0. 1 3724 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 

Space Wave Radiation Pattern (Phi=0) 




wel^ , kIOj 

Max E-Field Intensity (Volts per meter) max( MagElO ) = 0.05843 
Surface Wave Radiation Pattern (Phi=0) 




-1 - 0.5 0 0.5 1 

wetj. . x20j 



Max E-Field Intensity (Volts per meter) max(MagE20) = 4.46943- 10 * 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 

. h30j 



Max E-Field Intensity (Volts per meter) max(MagE30) = 0.05843 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




-1 - 0.5 0 0.5 1 

Hfet^ , KlHj 

Max E-Field Intensity (Volts per meter) max(MagElH) = 0.09414 
Surface Wave Radiation Pattern (Horizontal) 




-1 - 0.5 0 0.5 1 

. k2Hj 



Phi=0 



Max E-Field Intensity (Volts per meter) 



max( MagE2H) = 5.572%- 10 




Combined Space and Surface Wave Radiation Pattern (Horizontal) 




-1 - 0.5 0 0.5 1 

, «3Hj 



Max E-Field Intensity (Volts per meter) max(MagE3H) = 0.09414 



Phi=0 
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power# = 



30 








* 1.5707963 


M .5707963 




E 










y 


► 


- 1.5707963 « 


0 







cos ( £ ) ■ cos ( J ) • 1^ j-p-Y csc(1 

. e ^ 



1 - (sin(O cos(|))‘ 



powers = 



30 



x R* 



x 

2 



-- -0 



x 

2 



£ 



sin ( £ ) • I 



jpY Y -csc(^)(cos(C)cos(H 



1 - (sin(f) cos(O)' 



power := power# + powers 

2 



Radres - 2 power 



( | MagE3 p | | ) 1 4*R°max(sqE3P) 

sqE3Pj = ^ _ / DirectivityP = / 



2- ( 120- x) 



power 



c „ ( | MagE30j | ) 2 4 * R 2 max(sqE30) 

sqE30j ;= --- - ---- — '-f— DireclivityO = \ m / 



2( 120x) 

Total Power Radiated (Watts) 



power 



Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 



power = 261.928 
Radres = 523.856 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP = 10.78579 Directi vityO = 1.95488 
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Effective Isotropic Radiated Power (EIRP) (Watts) 



Phi - pi/2 Plane Phi - 0 Plane 

DirectivityP* power = 2.8251*10'* DirectivityO* power = 512.03765 

Maximum Effective Area (Along Radial of Directivity) 

(square meters) 

Phi = pi/2 Plane Phi = 0 Plane 



( Ae) ^* DirectivityP 

V V = 536.4412 

4*x 



(A^)^* DirectivityO 
4*x 



97.22772 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres* ( Ae) • DirectivityP 

V 7 = 54.60488 2- 



480* 



Radres* ^A^ ^* DirectivityO 



480* x‘ 



23.24692 



Numerical Distances 



Vertical Polarization 



Pe Q I = 0.06283 



Horizontal Polarization 
Pm 0 | = 2.26211- 10 6 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 

Phi = pi/2 Plane Phi = 0 Plane 



AngleP = 34.85714 



AngleO = 



/ 56.57143 
\ 56.57143 
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APPENDIX F; 

HORIZONTAL YAGI-UDA ARRAY COMPUTER OUTPUT 



This appendix contains computer hardcopies from the 
Mathcad horizontal Yagi-Uda array application which show the 
input values and predicted radiation characteristics for two 
sample calculations. The configuration is a given by the 
inputs on the first three pages of each printout. The first 
antenna is mounted above soil (e r =4 and c=10' 3 ) and the second 
above seawater (e r =72 and a=4) . Reference 6 [p. 107] provides 
the radiation patterns and gain predictions from several 
sources for the configuration in the first example. 

The radiation patterns and maximum directive gain 
computed by the first Mathcad example are almost identical to 
the those given in reference 6. The directivity prediction is 
slightly higher for Mathcad, but the overall similarity among 
the predicted radiation patterns is very good. The seawater 
example has results very much like the soil example. For 
horizontal polarization, the higher conductivity surface below 
the antenna does not result in a greatly enhanced surface wave 
and increased directivity as it does for vertical 
polarization. 
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HORIZONTAL YAG1-UDA ARRAY 



This application calculates far field radiation patterns and parameters associated with 
horizontal Yagi-Uda arrays. The antenna is oriented such that the projection of its 
center axis lies on the positive y-axis of a rectangular coordinate system. The feed 
is at the center of a single driven element between a set number of longer parasitic 
reflector elements and shorter parasitic director elements. Antenna elements are 
bisected by the antenna axis and are parallel to the x-axis. The center of the first 
reflector is directly above the origin at a set height on the z-axis. Required inputs 
are the number of reflector and director elements, individual element lengths, 
individual element radii, separation between adjacent elements, height of the feed 
above the surface, transmitted frequency, distance from the antenna, conductivity 
and dielectric constant of the surface below the antenna. The planar earth model is 
assumed in predicting radiation patterns. A sinusoidal voltage with a one volt 
maximum across the input terminals is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna in the plane of interest. The electric field magnitudes to which the patterns 
are normalized are displayed below their respective plots. Radiation patterns are 
plotted for the phi=pi/2 and phi=0 vertical planes parallel and perpendicular to the 
y-axis respectively. A horizontal radiation pattern is plotted at an elevation 
selected by the index from the elevation angle index table. Polarization is a 
varying combination of horizontal and vertical depending on spatial orientation with 
respect to the antenna. 



NOTE: The Yagi-Uda mutual impedance calculations are not valid if there is an 
antenna element whose length is an exact integer multiple of the wavelength. If 
this occurs, there will be a singularity error in the mutual impedance calculations. 
This should not be a problem since Yagi-Uda Antennas are designed to operate at 
frequencies such that the driven element is approximately one-half-wavelenglh long, 
the reflectors are slightly longer, and the directors are slightly shorter. If the problem 
does arise, it will be necessary to vary the operating frequency such that no element 
is exactly an integer multiple of the wavelength. 



Input the number of = j 
reflector elements 



y := 0 .. ( NR +- ND) 



Input the number of j 

director elements 

6 = 0 .. ( NR + ND - 1 ) 
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This entire page is ailoted for these entries to allow for as many elements as necessary. 



Input the lengths and radii of the Yagi 
elements, starting with the outermost 
reflector and proceeding inward to the 
driven element. After entering the driven 
element length and radius, enter those of 
the director elements, starting with the 
director next to the driven element and 
proceeding outward to the final element. 



Input the distances between 
successiv elements in the same 
order as the length and radius 
inputs. There should be one less 
entry in the separation distance 
array than there is in the length 
and radius array. 



I = 
Y 



15.6 



15 



14 



rad = d„ - 

Y 5 



.001 




7.5 


.001 




6.0 


.001 
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Input the operating f g = 10- 10 6 
Frequency (Hertz) 

_ 3 

Input the ground o - 110 
Conductivity 



Input the height of H 0 = 7.5 
the array above the 
surface (meters) 

Input the ground «r := 4 
Dielectric Constant 



Input the Distance 
from the Antenna 
(meters) 



R := 3000 



Input the index of the elevation 

angle for which to calculate the w := 285 

horizontal radiation pattern 

(from the angle index table) 



Index of Refraction 



n •= 




j- 



18000 - o 



Complex Numerical Distance for Vertical Polarization 



Pe ; 



-j-3(R + H„.cos(»|)) 
2-sin^^ ^ 





2 



Complex Numerical Distance for Horizontal Polarization 

-j-/3-(R + H o -cos(0 | ')) [ rz 7 77 

Pnl i = 7 ~^2\ C0S W * J" ' ( S> "W ) 

2.(sto(0,) 2 ) 



Vertical Reflection Coefficient 



Horizontal Reflection Coefficient 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 



xref^ , xlP. 

Max E-Field Intensity (Volts per meter) max(MagElP) = 5.13224' 10 
Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 0.5 1 



. k2Pj 

Max E-Field Intensity (Volts per meter) max( MagE2P) = 4.08957* 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




, h3P| 

Max E-Field Intensity (Volts per meter) max( MagE3P) = 5.13132* 10 * 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 

Space Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 t 



x10| 

Max E-Field Intensity (Volts per meter) max(MagElO) = 1.35264-10 * 
Surface Wave Radiation Pattern (Phi=0) 




-1 -0.5 0 0.5 1 



xret^ . x20j 

Max E-Field Intensity (Volts per meter) max(MagE20) = 5.32341-10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




-1 - 0.5 0 0.5 1 

, h30j 



Max E-Field Intensity (Volts per meter) max(MagE30) = 1.35264*10 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 

Kfet^ . kIHj 



Phi=0 



Max E-Field Intensity (Volts per meter) max(MagElH) = 3.79643* 10 * 
Surface Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 



Kret^ , k2H. 

Max E-Field Intensity (Volts per meter) max( MagE2H) = 2.73697* 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 




-t -0.5 0 0.5 t 

, h3H. 



Max E-Field Intensity (Volts per meter) max( MagE3H) = 3.7942* 10 * 
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power# ■- 



30 

icR‘ 





-1. 5707963 


* 1.5707963 




Z 










y 




-1.5707963 - 


0 







cos{$)-cos(S ) I jpD sin(c: 

1 — . e y 



1 - (sin(f)-cos(l))' 



power 0 = 



30 

*R 5 



it 

2 



- ‘0 



x 

2 



E 



sin(01 



J p D -sh(C)-sln(^) 



1 - (sin(0’Cos(())‘ 



(cos^0- 



power ;= power# 4- power0 



Radres 



( vref mr ) 2 

2 • power 



sqE3Pj 



(iMagESPjl ) 2 
2 ( 120 *) 



Directi vityP := 



4*R 2 - max( sqE3P) 
power 



sqE30j 



( | MagE30 ( | ) 2 

2( 120x) 



DirectivityO 



4 • ic ■ R 2 - max ( sqE30 ) 
power 



Total Power Radiated (Watts) power = 0.00266 

Radiation Resistance (versus Radres = 187.94973 

Maximum Input Current) 

(Ohms) 

Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP = 14.84639 DirectivityO = 1.03163 
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Effective Isotropic Radiated Power (EIRP) (Walts) 

Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP- power = 0.0395 DireclivityO- power = 0.00274 

Maximum Effective Area (Along Radial of Directivity) 
(square meters) 

Phi = pi/2 Plane Phi = 0 Plane 



2 * DirectivityP 
4- % 



1.06329* 10 3 



( Aj) 2 - DireclivityO 
4-x 



73.88536 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres- ( Ac) • DirectivityP 
^ = 46.04814 2- 



480- 



Radres- ( Ac) • DireclivityO 
^ = 12.13848 



480- 



Numerical Distances 



Vertical Polarization 



Pe Q = 57.12622 



Horizontal Polarization 
Pidq I = 1.09911* 10 3 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 

Phi = pi/2 Plane Phi = 0 Plane 



AngleO = 



89.71429 

89.71429 



AngleP = 36 
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HORIZONTAL YAG1-UDA ARRAY 



This application calculates far field radiation patterns and parameters associated with 
horizontal Yagi-Uda arrays. The antenna is oriented such that the projection of its 
center axis lies on the positive y-axis of a rectangular coordinate system. The feed 
is at the center of a single driven element between a set number of longer parasitic 
reflector elements and shorter parasitic director elements. Antenna elements are 
bisected by the antenna axis and are parallel to the x-axis. The center of the first 
reflector is directly above the origin at a set height on the z-axis. Required inputs 
are the number of reflector and director elements, individual element lengths, 
individual element radii, separation between adjacent elements, height of the feed 
above the surface, transmitted frequency, distance from the antenna, conductivity 
and dielectric constant of the surface below the antenna. The planar earth model is 
assumed in predicting radiation patterns. A sinusoidal voltage with a one volt 
maximum across the input terminals is assumed. All radiation patterns are 
normalized with respect to the maximum electric field intensity transmitted by the 
antenna in the plane of interest. The electric field magnitudes to which the patterns 
are normalized are displayed below their respective plots. Radiation patterns are 
plotted for the phi=pi/2 and phi=0 vertical planes parallel and perpendicular to the 
y-axis respectively. A horizontal radiation pattern is plotted at an elevation 
selected by the index from the elevation angle index table. Polarization is a 
varying combination of horizontal and vertical depending on spatial orientation with 
respect to the antenna. 



NOTE: The Yagi-Uda mutual impedance calculations are not valid if there is an 
antenna element whose length is an exact integer multiple of the wavelength. If 
this occurs, there will be a singularity error in the mutual impedance calculations. 
This should not be a problem since Yagi-Uda Antennas are designed to operate at 
frequencies such that the driven element is approximately one-half-wavelength long, 
the reflectors are slightly longer, and the directors are slightly shorter. If the problem 
does arise, it will be necessary to vary the operating frequency such that no element 
is exactly an integer multiple of the wavelength. 



Input the number of ^ = j 
reflector elements 



Input the number of = j 

director elements 



y := 0 .. ( NR * ND ) 



6 = 0 .. ( NR ^ ND - 1 ) 
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rhis entire page is alloted for these entries to allow for as many elements as necessary. 



Input the lengths and radii of the Yagi 
elements, starting with the outermost 
reflector and proceeding inward to the 
driven element. After entering the driven 
element length and radius, enter those of 
the director elements, starting with the 
director next to the driven element and 
proceeding outward to the final element. 

1 = rad := 

y y 



15.6 




.001 


15 




.001 


14 




.001 



Input the distances between 
successiv elements in the same 
order as the length and radius 
inputs. There should be one less 
entry in the separation distance 
array than there is in the length 
and radius array. 




7.5 

6.0 
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Input (he operating fg = 10* 10^ 
Frequency (Hertz) 



Input the ground o - 4 
Conductivity 

Input the Distance ^ _ ^qoq 

from the Antenna 

(meters) 

Input the index of the elevation 
angle for which to calculate the w ~ 
horizontal radiation pattern 
(from the angle index table) 



Input the height of Hq = 7.5 
the array above the 
surface (meters) 

Input the ground sr := 72 
Dielectric Constant 

Index of Refraction 



18000 - o 




Complex Numerical Distance for Vertical Polarization 




- j*0' (R + H o -cos(0j)) 

2-sin(0j^ 




n - sm 



W 



Complex Numerical Distance for Horizontal Polarization 

- J-zS-fR + Hq- cos ( 6 ,\ ") [ rZ 7 77 

Pm i ’ iv. VV'F^ • c0, ffl * i" - ( sin ( 9 i) ) 

2-(sm(9 | ) z 



Vertical Reflection Coefficient 



Horizontal Reflection Coefficient 





COS^0|^ - 
COS^0|) +■ 
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Radiation Patterns in Phi=Pi/2 Plane 
(Perpendicular to X-Axis and Dipole Elements) 
Space Wave Radiation Pattern (Phi=Pi/2) 




-1 -0.5 0 0.5 1 

xref^ , xIP. 



Max E-Field Intensity (Volts per meter) max(MagElP) = 6.85595*10 * 

Surface Wave Radiation Pattern (Phi=pi/2) 




-1 -0.5 0 0.5 1 

xref^ . x2Pj 



Max E-Field Intensity (Volts per meter) max(MagE2P) = 2.16533*10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=pi/2) 




, h3Pj 

Max E-Field Intensity (Volts per meter) max(MagE3P) = 6.85595*10 * 
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Radiation Patterns in Phi=0 Plane 
(Perpendicular to Phi=pi/2 Plane) 

Space Wave Radiation Pattern (Phi=0) 




-1 - 0.5 0 0.5 1 



kIOj 

Max E-Field Intensity (Volts per meter) max(MagElO) = 2.32256- 10 * 
Surface Wave Radiation Pattern (Phi=0) 




, k20j 

Max E-Field Intensity (Volts per meter) max( MagE20) = 9.86%5- 10 ^ 
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Combined Space and Surface Wave Radiation Pattern (Phi=0) 




Max E-Field Intensity (Volts per meter) 



max( MagE30) = 2.32256-10 
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Radiation Patterns in Horizontal Plane 
(Parallel to Ground) 

Space Wave Radiation Pattern (Horizontal) 




Phi=0 



nret^ . hIHj 

Max E-Field Intensity (Volts per meter) max(MagElH) = 4.25018*10 
Surface Wave Radiation Pattern (Horizontal) 




Phi=0 



, k2H ( 

Max E-Field Intensity (Volts per meter) max(MagE2H) = 3.10467*10 
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Combined Space and Surface Wave Radiation Pattern (Horizontal) 




-1 -0.5 0 0.5 1 

, h3H. 



Max E-Field Intensity (Volts per meter) max( MagE3H) = 4.25018- 10 



Phi=0 
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power# ;= 



30 

*R J 





• 1.5707963 


' 1.5707963 




E 










y 


- 


- 1.5707963 - 


0 







cos ( £ )cos( f )*I 

1 - (sm(?)-cos(€))‘ 



jp-D sin(c: 

„ y 



power# 



30 

*-R : 



* 

2 



-- >0 



% 

2 






sin( 0*1 



1 - (sin(O cos( O)' 



jp-D jm(;) jin(4) 
e * *^cos^/3- 



power = power# +■ powers 



Radres 



(vm*hb ) 2 

2 power 



sqE3Pj 



(| MagE3Fj|)^ 
2 -( 120 *) 



DirectivityP = 



4*R^-roax( sqE3P) 
power 



sqE30j 



(|MagE30 ( |) 2 
2 ( 120 *) 



DirectivityO 



4*R^max( sqE30) 
power 



Total Power Radiated (Watts) power = 0.00482 



Radiation Resistance (versus 
Maximum Input Current) 
(Ohms) 



Radres = 103.76004 



Directivity (or Maximum Power Gain assuming 100% Antenna Efficiency) 
Phi = pi/2 Plane Phi = 0 Plane 

DirectivityP = 14.63143 DirectivityO = 1.67913 
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Effective Isotropic Radiated Power (EIRP) (Watts) 

Phi = pi/2 Plane Phi = 0 Plane 

DireclivityP- power = 0.07051 DirectivityO- power = 0.00809 

Maximum Effective Area (Along Radial of Directivity) 
(square meters) 



Phi = pi/2 Plane 

2 r»* 



Phi = 0 Plane 



( • DireclivityP ~ • DirectivityO 

v = 1.0479* 10 J — : = 120.25911 



4-x 



4-x 



Maximum Effective Length (Along Radial of Directivity) 

(meters) 

Phi = pi/2 Plane Phi = 0 Plane 



Radres- ( Xr) • DireclivityP 
= 33.96559 2* 



4 



480- x‘ 



Radres- ( X *) DirectivityO 

^ = 11 . 



50637 



480- it 4 



Numerical Distances 



Vertical Polarization 



Pe Q I = 0.04363 



Horizontal Polarization 
Pm Q | = 2.26206* 10 6 



Elevation Angle of Directivity (Maximum Gain) 
above the Horizon (Degrees) 



Phi = pi/2 Plane 



Phi = 0 Plane 



AngleP = 41.71429 



AngleO = 



89.71429' 

89.71429; 
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